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OXIDATION OF INTERMEDIATES OF THE TRICARBOX YLIC 
ACID CYCLE BY EXTRACTS OF AZOTOBACTER AGILE* 


By Roy REPASKET AND P. W. WILSON 


INSTITUTE FOR ENZYME RESEARCH AND DEPARTMENT OF BACTERIOLOGY, UNIVERSITY 
OF WISCONSIN, MADISON 


Communicated by E. B. Fred, February 19, 1948 


Investigating the terminal carbohydrate metabolism of Azotobacter 
vinelandii 0, Stone and Wilson! * demonstrated by manometric and isotopic 
methods that the tricarboxylic acid (TCA) eyele is a primary path of oxida 
tion. Karlsson and Barker*® concluded from the results of simultaneous 
adaptation and isotopic dilution experiments with whole cells of A. agile 
1.4 that acetate is not oxidized by any of the conventional metabolic path 
ways. As two organisms so closely related would not be expected to pos 
sess completely different terminal oxidative mechanisms, the oxidation of 
intermediates of the TCA eycle by A. agile has been reinvestigated using 
cell-free extracts as well as whole cells. 

Methods... Azotobacter agile 4.4 was grown on Burk’s nitrogen-free min 
eral salts medium with 2°; sucrose as the source of carbon. ‘The cultures 
were grown in 500-ml. Erlenmeyer flasks containing 100 ml. of medium 
that had been inoculated with 5°, of an active 24-hr. culture and incubated 
for 15 hrs. at 80°C. on a rotary shaker; both the quantity of inoculum and 
the incubation time were selected to give a large yield of young cells with a 


high protein content about 12 grams of wet cells per liter containing 5-6", 


protein. The stock culture was maintained by 24-hr. serial transfer. 

Details of the methods for harvesting the cells and for preparing the cell 
free extracts are described by Wilson and Knight.! Several methods for 
making the enzyme preparations were tested, including grinding in a 
Booth-Green mill, with glass beads and with alumina, but none was as 
satisfactory as breaking in a Raytheon 10 Ke sonie oscillator. 

Twelve grams of wet cells were suspended in a total volume of 40 ml. 
with cold 0.1 1/7 KHCOs,, pH 8.4, and exposed for 10 minutes in the sonic 
oscillator cooled by circulating ice water. The endogenous respiration of 
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such a suspension is usually 100-200 ul O, per hr. per mg. N, but after 
dialyzing 5 hrs. against two changes of 0.04 \/ phosphate buffer, pH 7 
(two liters), the endogenous activity drops below 50 ul. The preparations 
contained 1.4 to 2.2 mg. biuret protein nitrogen per ml.; extracts with less 
than 1.2 mg. N per ml. are unsatisfactory. In the experiments reported 
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Oxidation of intermediates of tricarboxylic acid cycle by whole cells of 
Azotobacter agile 4.4.. Flask contents: substrate as indicated; 6 uM. Mg 
phosphate; 40 hM. PO, pH 7; 1 ml. of cell suspension (0.179 mg. N/m.) 
and H.O to 3.0 ml; 0.2 ml. 20°, KOH in center well. Temperature of bath 


33°C, 


here the extracts were tested the same day as prepared; they can be kept 
overnight at — 10°C. without appreciable loss in activity. 

All experiments were made by conventional manometric procedures; 
oxygen uptake was measured at 33°C., the initial reading was made after 5 
minutes of equilibration after which the substrate was tipped. Details 
of contents of the fiasks are supplied in the tables and figures. The re- 
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agents were commercial products with the exception of the CoA (20%, pan- 
tothenic acid) kindly supplied by Dr. Helmut Beinert and the TPN (tri- 
phosphopyridinenucleotide) prepared by the method of Kornberg and 
Horecker.° 

Results. Figure 1 illustrates typical results obtained with whole cells 
of Azotobacter agile 4.4. A period for adaptation was required to oxidize 
all members of the TCA cycle except acetate, which was immediately 
oxidized, and citrate, which was not oxidized even after an hour. An 
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FIGURE 2 

Oxidation of intermediates of tricarboxylic acid cycle by cell-free extracts 
of Azotobacter agile 4.4. Flask contents: 2 w6M. MgCh, 500 pg. DPN, 16 ug 
TPN, 19 ug. CoA, 30 uM. glutathione, | ml. enzyme (RDI1; 2.09 mg. N/m), 
3M. substrate and water to 3.0 ml.; 0.2 ml. 20°) KOH in center well. Tem- 
perature of bath 33°C. Corrected for endogenous (16 wl. Ov/hr. ) 


entirely different result is obtained, however, when cell-free enzyme prepa 
rations are used (Fig. 2). Although the whole cells required adaptative 
periods of 15 to 60 minutes before the maximum rate of oxidation of vari 
ous substrates is reached, extracts from these cells immediately oxidized 
all the TCA intermediates tested. This is most strikingly demonstrated 
with citrate; oxidation of isocitrate occurred at an even greater rate (not 
shown). 


The dialyzed extracts must be supplemented with certain coenzymes for 
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maximum rates of oxidation with citrate, a-ketoglutarate and acetate. 
The quantities required under the conditions of our experiments are indi- 
cated in the captions for figures 2 and 3. The oxidation of isocitrate by 
enzymes from animal and yeast cells requires TPN. Extracts of A. agile 
slowly oxidize citrate in the presence of 500 wg. DPN (table i), but when as 
little as 4 ug. TPN is also added the Qo,(N) is increased almost threefold; 
16 ug. TPN per vessel caused a slightly higher rate of oxidation. Although 


SPARKING OF ACETATE 
9M ACETATE 
1 ym SUCCINATE 


va 


BL o, UPTAKE 
S 
° 


° 
1 pM SUCCINATE 


Fi | ——— 
ee 
P i 


V 


° o— 0 


_  O0 
ane 9 JM ACETATE 
—_ 
"ete 
10 20 30 40 60 
MINUTES 
FIGURE 3 








Oxidation of acetate by cell-free extracts of Asotobacter agile 4.4. Flask 
contents: | ml. extract (RD9; 1.78 mg. N/ml.), 2 o6M. MgCh, 500 pz 
DPN, 16 wg. TPN, 19 wg. CoA, 80 uM. glutathione, 4 uM. ATP, substrate as 
indicated, and water to 3.0 ml; 0.2 ml. 20°, KOH in center well. Tempera- 
ture of bath 33°C. Corrected for endogenous (18. wl. Oo/hr.), 


experiments were conducted with citrate as substrate, isocitric dehydro- 
genase activity had been demonstrated, and it is probable that it is the 
isocitric dehydrogenase of A. agi/e that is TPN-linked. Coenzyme A 
also stimulated the oxidation of citrate, but the effects may arise from the 
CoA requirement for oxidation of a-ketoglutarate, presumably the product 
of citrate oxidation. 

DPN, CoA and Mgt* were all found to increase the rate of oxidation of 
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a-ketoglutarate (table 2). Lindstrom® has demonstrated these same re- 


quirements for oxidation of this compound by a soluble a-ketoglutaric 
oxidase from Astobacter vinelandii, and Sanadi and Littlefield’ found CoA 
and DPN are required for the oxidation of a-ketoglutarate by a soluble 
The dialyzed extract from A. agile 


preparation from animal tissue. 
DPN is 


oxidizes a-ketoglutarate slowly without added DPN; if 50 ug. 


TABLE 1 


OXIDATION OF CITRATE BY EXTRACTS OF A sotobacter agile 
wl. Or UPTAKE? . 
Qoo(N) 30 min 60 min 


PREPARATION SYSTEM 


RD Complete" 
$ wg. TPN 191 157 256 
Minus TPN 69 63 133 
156 144 200 
5A 
111 


249 215 340 


RDI1O Complete 
Minus TPN 47 $5 
Minus CoA 56 62 


“ Complete system contains 16 wg. TPN; for other constituents see figure 2. 


» Figures in all tables corrected for endogenous respiration. 


RDO 1.78 mg. N/ml.; 10 uM. citrate per vessel; endogenous, 18 ul. Oo/hr. 


RDIO 1 86mg. N/ml.; 3 uM. citrate per vessel; endogenous, 21 yl. O./hr. 


TABLE 2 


OXIDATION OF @-KETOGLUTARATE BY EXTRACTS OF A. agile 

ul. Oo UPTAKE . 

SYSTEM Qo N) 30 MIN 60 MIN 
RDS Complete“ 199 140 174 
50 wg. DPN 85 58 82 


Minus DPN 42 32 55 


PREPARATION 


RDO Complete 137 117 136 
Minus DPN 51 47 75 


Minus CoA 50 43 75 
Minus Mg** 83 79 127 


* Complete system contains 500 wg. DPN; for other constituents see figure 2, 


RDS 1.70 mg. N/ml.; endogenous, 17 wl. O./hr. 
RDO 1.86 mg. N/ml.; endogenous, 21 wl. O./hr. 


added per vessel, the Qo,(N) is increased slightly but still is only about 40°; 
of that obtained with 500 ug. DPN per vessel. Since the reactive group 
of CoA is believed to be a sulfhydryl group, glutathione was used as a re- 
ducing agent (cysteine is itself oxidized and thus raises the apparent endog- 
In the absence of added DPN or CoA the rate of a- 
CoA levels above 12-15 


enous respiration). 
ketoglutarate oxidation 1s one-third of maximum. 
ug. per vessel did not increase the oxygen uptake. 

Data in table 2 also demonstrate the stimulatory effect of Mg++. The 
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over-all oxygen uptake per hour is about the same as in the complete sys- 
tem, but the rate of oxidation is more rapid with added Mg**. ATP, 
cocarboxylase, and protogen have no effect on a-ketoglutarate oxidation. 
Acetate alone, although immediately and rapidly oxidized by whole 
cells, is not oxidized to any extent by cell-free extracts (table 3). With 
the addition of small quantities of succinate (1 uM.), the net oxygen up- 
take in an hour corrected for oxidation of succinate alone is 5-6 times the 
uptake of acetate alone. As can be seen from figure 3, this “‘sparking’’ of 
acetate still continues at a rapid rate after the oxidation of the individual 
substrates has leveled off. 
TABLE 3 
ACETATE OXIDATION BY EXTRACTS OF A. agile 
ut. O 
UPTAKE 
PREPARATION SYSTEM SUBSTRATE QouwN 60 MIN 
RD12 Complete’ 9 uM. Acetate 11 20 
9 uM. Acetate + 104 158 
1 uM. Succinate 
1 uM. Succinate ¢ 64 
No ATP 9 uM. Acetate 0 
9 uM. Acetate + 
1 uM. Succinate 
1 uM. Succinate 
RDO Complete 9 uM. Acetate 
9 uM. Acetate + 
1 uM. Succinate 
1 uM. Succinate 
No CoA 9 uM. Acetate 
9 uM. Acetate 4 
1 uM. Succinate 
1 uM. Succinate 


“ See figure 3 for complete system. 


RDIZ 2.12 mg. N/ml.; endogenous, 18 ul. O./hr 
RDO 1.78 mg. N/ml.; endogenous, 21 wl. Oo/hr 


Discussion. Two terminal oxidative cycles are generally postulated for 
bacteria, the dicarboxylic and the tricarboxylic,*® although the data of 
Karlsson and Barker with whole cells of Azotobacter agile suggest that 
neither of these functions in this organism. The chief difference between 
the di- and tricarboxylic acid cycles is in the mechanism for formation of 
succinate: whether directly from two moiecules of acetate or via citrate. 


One reason for rejecting the operation of the TCA cycle in A. agile has been 
the apparent inability of cells which immediately oxidize acetate to oxidize 
other members of the cycle without a period of adaptation (principle of 
simultaneous adaptation). However, this particular objection need no 
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longer be considered a serious handicap since such cells do possess the neces- 
sary enzyme systems, 1.e., cell-free preparations readily oxidize the inter- 
mediates. 


Further evidence supporting the existence of the TCA cycle in A. agile 


is the similarity of coenzyme and metai requirements for the oxidation of 
acetate, citrate, and a-ketoglutarate by extracts of these cells and these 
requirements in preparations from animal, yeast, and other bacterial spe- 
cies. 

Finally, the ‘‘sparking’’ of acetate oxidation by succinate is strong evi- 
dence for the participation of a TCA cycle.’ Either oxalacetate or a pre- 
cursor must be available to condense with acetate; thereafter, oxalacetate 
is regenerated at each turn of the cycle. A requirement for a dicarboxylic 
acid would not be expected if two molecules of acetate condensed to form 
succinate. Also, it should be noted that Stern ef a/.'° observed in 1950 
that A. agi/e possesses the enzyme for condensing acetate with oxalacetate. 
Recently the reactions for formation of citrate from acetate have been 
clarified ;'') '* in light of this work it seems significant that extracts of A. 
agile readily oxidize acetate only when ATP and CoA are present (table 3). 

Although it is recognized that manometric experiments alone are not 
sufficient to demonstrate conclusively a metabolic pathway, it is empha- 
sized that all results with cell-free extracts of A. agi/e are compatible with 
the view that the TCA cycle is a terminal oxidation pathway in this or- 
ganism. These results do not agree with conclusions based on experi- 
ments with cell suspensions but do agree with analogous (and more exten- 
sive investigations) made with the closely related organism, A. vinelandit. 

Summary. The ability of unadapted cells of Azotobacter agile 4.4 and 
of cell-free extracts from such cells to oxidize intermediates of the tricar- 
boxylic acid cycle has been compared. Adaptive periods were required 
before whole cells oxidize all substrates except acetate and citrate. Ace- 
tate was immediately oxidized; citrate was not oxidized even alter an 
hour. Cell-free extracts, however, oxidized all these substrates immedi- 
ately. 

Dialyzed extracts of A. agile 4.4 had coenzyme requirements similar to 
those for the corresponding enzymes of the tricarboxylic acid cycle found in 
amimal tissue. These included: TPN for citrate oxidation; CoA, DPN, 
and Mg** for a-ketoglutarate oxidation; and Mg**, CoA, and ATP for 
acetate oxidation. Acetate is not oxidized by these extracts unless it 1s 


+ 


“sparked.” 

All these findings are compatible with the conclusion that a tricarboxylic 
acid cycle functions in A. agile 4.4 as has been previously established for 
A. vinelandii. 


Acknowledgment. We wish to express our appreciation for the many 
helpful suggestions offered by Dr. D. E. Green. 
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A BIOCHEMICALLY SPECIFIC METHOD FOR 
ENZYME ISOLATION 


By Leonarp S. LERMAN* 
INSTITUTE OF RADIOBIOLOGY AND Bropuysics, THE UNIVERSITY OF CHICAGO 
Communicated by Linus Pauling, January 29, 1953 


It is at present accepted that the action of an enzyme requires as a first 
step the temporary combination of the substrate with a specific site on the 
enzyme, and that competitive inhibitors act by reversibly substituting for 
the substrate in combination with the specific site. If then, a specific 
competitive inhibitor or a suitable substrate of an enzyme is coupled to the 
surface of an inert solid, it may be expected that the enzyme will be re- 
versibly bound to the solid through combination with the attached group, 
while enzymes with different specificities or other proteins will not be bound. 
The enzyme could later be released by a change in the medium to condi- 
tions unfavorable for combination, as, for example, change in pH or salt 
concentration. Alternatively, by the use of a solution of a competitive 
inhibitor (not necessarily the same as that coupled on the adsorbent) the 
enzyme will, according to the principle of mass action, be returned into 
solution by an additional specific process. Thus the isolation of the enzyme 
would depend directly on its catalytic specificity rather than, as in conven- 
tional procedures, its properties as a protein. 

Mushroom tyrosinase was chosen for an experimental investigation be- 
cause of convenience and the simplicity of the associated chemical proced- 
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ures. With suitable adsorbents a significant purification of the enzyme 
was achieved. 

Experimental. 7 yrosinase: An aqueous extract of mushrooms minced 
in cold acetone was precipitated with ammonium sulfate, according to the 
procedure of Mallette, ef a/.;'_ the precipitate was redissolved and dialyzed 
thoroughly against distilled water. The resulting enzyme solution, usu- 
ally diluted with malonate buffer, pH 5.8, was used in all experiments 
For measurement of activity, aliquots of 0.05 ml. were incubated at 4°C. 
with 0.50 ml. of a solution containing 0.1 per cent catechol, 0.5 per cent 
gelatin, and 0.5 per cent piperidine in 0.1 1/7 malonate buffer, pH 5.8. 
After twenty-four hours the mixture was diluted with 1.00 ml. of water, 


TABLE 1 
Cellulose < NN 
Cellulose 
Cellulose 
Cellulose 
Cellulose 


NN 


Cellulose NN COOH 


Cellulose NN—«< NN COOH 


OH 


and the optical density of 390 my was measured in a 1.00 cm. cell. Because 
of the well-known reaction inactivation’ of tyrosinase, the optical density 
reaches a stationary value, which remains unchanged for at least another 
day. Its dependence on the amount of enzyme was found to be linear up to 
a density of 0.92, with a standard error of estimate of less than 0.01. One 
unit of enzyme activity was arbitrarily defined as the amount which, 
under the above conditions, would yield a density of 42.4. Ludwig and 
Nelson’ found a similar stationary value and proportionality to the amount 
of enzyme in the total oxygen uptake. Relative measurements of total 
protein concentrations were made by the colorimetric ninhydrin method of 


Kunkel and Ward. ! 
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Adsorbents..- Since it was found that azodyes prepared by coupling an 
aromatic amine with phenol inhibit tyrosinase, various adsorbents, listed 
in table 1, containing p-azophenol groups and related groups were pre- 
pared from aromatic ethers of cellulose. | Powdered cellulose (Solka- 
Floe BW200°) was partially etherified by one of two procedures: (a) reflux- 
ing in benzene with phosphorus tribromide, about two grams per hundred 
grams of cellulose, followed by refluxing in alcohol with an excess of sodium 
m-hydroxyphenoxide, or (b) refluxing with sodium ethylate and p-nitro- 
benzyl bromide, about two grams per hundred of cellulose, followed by 
reduction with sodium dithionite. The original fibrous character was not 


Activity 


Enzymatic 





Effluent fractions 
FIGURE 1 
Frontal chromatograms obtained from the percolation of diluted crude tyrosinase 


through a 0.20-g. column of adsorbent III or a similar column of adsorbent V. The 
fractions were 0.05 ml. each, and the holdup volume was about 0.45 ml. 
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FIGURE 2 
Frontal chromatograms from 0.12-g. columns of adsorbents IV and III. The frae 
tions were 0.05 ml. each, and the holdup volume was about 0.35 ml 


appreciably changed. Coupling of the diazonium salts was indicated by 
the usual color changes; the colors were fast to washing with dilute acids, 
alkalis, ete. Similarly constituted dyes which were not bound through 
the ether linkage were found to be easily washed away. 

enzyme Adsorption. Demonstration of adsorption by the cellulose de- 
rivatives, tests on specificity, and rough estimation of the capacity of the 
adsorbents were carried out by the technique of frontal analysis. As the 
enzyme solution percolated continuously through a column of the adsorb- 
ent, successive fractions of the effluent were collected and tested for ac- 
tivitv. In the absence of adsorption the enzyme appeared immediately 
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after the interstitial fluid in the column (the holdup volume) was dis- 
placed by the entering solution. When adsorption occurred, the enzyme 
appeared later according to the capacity of the adsorbent. For quantita- 
tive experiments an aliquot of crude enzyme, usually just enough to satu- 
rate the adsorbent, was introduced inte a column followed immediately 
by buffer of the same pH and later by the eluent. 

Results.— Frontal chromatograms (Figs. | and 2) at pH 5.8 or 6.3 showed 
cellulose p-azophenol preparations I, II, and III to be effective adsorbents, 
while no retention of the enzyme by IV, V, VI, or VII was found except 
for the slight non-specific adsorption shown also by untreated cellulose. 
With certain of the crude tyrosinase preparations a small, well-defined 
shoulder in the frontal chromatograms previous to the emergence of the 
maximum concentration of the enzyme was noted, indicating the presence 
of at least a second catalytically active substance. Some dependence on 
salt concentration was found, in that dilution of the malonate to 0.001 7 

PABLE 2 
ADSORPTION OF 0.30 ML. CRUDE ENZYME SOLUTION ON FIVE GRAMS OF ADSORBEN1 
III, ru&k CoLUMN THEN WASHED WITH MALONATE AND ELUTED WITH PYROPHOSPHATE 
ACTIVITY PROTEIN 
PER ML PER ML 


(ARBITRARY (ARBITRARY VOLUME FRACTION OF 
UNITS) UNITS) MI TOTAL ACTIVITY 


Crude enzyme 1.00 122.4 0.30 1.00 
Effluent from wash 0.08 
Eluate fractions 

2+3 0.001 0.00 

4 5 0.025 0.264 0.08 
6 0.106 0.214 0.35 
8 ) 0.032 0 007 0.11 
1+ 11 0 007 , 0.02 


+ 
+ 
+ 


¢ 


Total recovery in eluate 0.56 


inhibited adsorption, and addition of sodium chloride to 0.5 7 greatly in 
creased non-specific adsorption. However, in the vicinity of 0.04 to 
0.1 Af near pH 6 malonate, phthalate, and chloride were satisfactory and 
interchangeable. The addition of considerable amounts of extraneous 
protein tended to inhibit adsorption. 

On the passage through the column of a small sample of crude enzyme 
followed by buffer, nearly all of the original protein, now devoid of enzyme 
activity, was found in the effluent. The enzyme had remained on the 
adsorbent, and the enzyme was retained through protracted washing with 
malonate but was rapidly eluted by 0.02 17 pyrophosphate at pH 8.3, 0.1 
M glycine at pH 9.6, or other buffers of similar pH. Elution was also 
effected by a solution of sodium benzoate, a strong competitive inhibitor, 
at pH 6.3. Data on a single process of adsorption, washing and elution 
by pyrophosphate are given in table 2... Approximately 56 per cent of the 
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original enzyme was recovered in all fractions of the eluate; the most active 


fraction contained 62 per cent of the total recovered. The ratio of en- 
zyme activity to protein concentration in that fraction was 61 times as 
great as in the original crude preparation. Since highly purified prepa- 
rations of tyrosinase are more rapidly inactivated than crude preparations 
while acting on catechol,’ a somewhat higher value might have been indi- 
cated by a different method of enzyme assay. 

Continuing study of the method is directed toward determining its ap- 
plicability to the purification of various enzymes, improving the capacity 
of the adsorbents, and achieving a degree of purification closer to that 
theoretically possible. 

* Schenley Fellow, 1949-1951. Present address: Laboratory of Chemical Embry 
ology, University of Colorado Department of Medicine, Denver, Colo. 

' Mallette, M. F., etal., Arch. Biochem., 16, 283 (1948) 

* Nelson, J. M., and Dawson, C. R., Advances in Enzymology, 4, 99 (1944). 

’ Ludwig, B. J., and Nelson, J. M., J. Am. Chem. Soc., 61, 2601 (1939). 

' Kunkel, H.G., and Ward, S. M., J. Biol. Chem., 182, 597 (1950), 

> Provided by the Brown Company, Chicago 


THE SOLUBILITY OF THE DIMETHYLETHERATE OF 
ALUMINUM BROMIDE IN METHYL BROMIDE IN THE 
PRESENCE OF SATURATED ITYDROCA RBONS* 


By LutHerR M. FosteER AND CHARLES A. KRAUS 
METCALF RESEARCH LABORATORY, BROWN UNIVERSITY 
Communicated February 4, 1953 


In the course of an investigation of the polymerization of ethylene in 
methyl bromide in the presence of aluminum bromide at —33°C., it appeared 
that the ethylene, or its polymerization products, did not complex with 
all the aluminum bromide present, although much more than one mole of 
ethylene reacted per mole of aluminum bromide.' At the same time, the 
hydrocarbon recovered on decomposing the complex with water was of 
higher molecular weight than expected, even though the proportion of 
ethylene to aluminum bromide was little greater than unity. 

It was thought that the uncombined aluminum bromide might be sepa- 
rated out by adding dimethyl ether which forms the stable complex 
AlBrs:(CH;)2O, whose solubility in methyl bromide ts relatively low, partic- 
ularly at lower temperatures. However, the solubility of the etherate 
was still rather high, even at —7S°C. Ina later investigation, propane was 
added to the methyl bromide solution of the etherate with a view to reducing 
its solubility. The solubility of the etherate in pure propane is negligible. 
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On carrying out some preliminary experiments, the expected large solubil- 
ity decrease on addition of propane was not found and, moreover, singulari- 
ties appeared in the solubility curve which could hardly be ascribed to 


experimental error. Accordingly, a series of careful measurements was 
carried out with the propane-methyl bromide mixtures at —7S°C. These 
measurements established the correctness of the earlier observations. Meas 
uremi.:'ts were then carried out at the same temperature with n-butane, 
isobutane and n-pentane as diluents. 

All materials employed were carefully tested and, if necessary, purified. 
Aluminum bromide was distilled three times under reduced pressure; the 
product was colorless. 

After purification, the vapor pressure of methyl bromide was measured 
over the temperature range —S9.S° to +40°C. The experimental values 
are given to 0.25% by the empirical equation: 

— 588.517 
logiwp = r — 2.36660 logy 7 + 14.3997. 
The boiling point as interpolated was found to be 3.66 + 0.05°, in good 
agreement with 3.58 + 0.05°, as reported by Egan and Kemp.* The freez 
ing point was found to be —93.65 + 0.05° in agreement with —93.66 + 
0.05° as reported by Egan and Kemp. 

Dimethylether, said to have a purity of 99.9°,, was passed over activated 
alumina. At 25° and 760 mm., the vapor was found to have a density of 
1.918 g. 1. in agreement with the value reported by Kennedy, Segenkahn 
and Aston.* To establish the purity of the propane, its vapor pressure was 
measured between —95.50 and —42.00°. Experimental values are re 


produced within 0.25°, by the empirical equation 


logy p = — rat? neeees logw 7 + 21.2951. 
Tr 

By interpolation, the boiling point was found to be —42.12 + 0.05° in 
agreement with Rossini's! value of —42.1 

All solubility measurements were carried out at —7S8.0S + 0.05°C. in a 
dry ice-acetone bath whose temperature was manually controlled. Tem 
peratures were measured by means of a copper-constantan thermocouple, 
the tip of which was encased in a thin-walled glass tube and immersed in 
the solution. The thermocouple was calibrated by means of standard melt- 
ing point samples and could be read to 0.02 

The solutions were prepared in a tube having a capacity of approximately 
150 ce.; it was provided with a magnetic stirrer which was operated manu 
ally. A weighed sample of aluminum bromide in a fragile bulb was in 
troduced into the solubility tube and about 60 cc. of methyl bromide was 
condensed on the bromide. When the bromide was in solution, dimethyl 
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ether, in approximately 10% excess, was condensed in the solubility tube. 
The desired amount of hydrocarbon was then distilled into the system. 
Thereafter, the thermocouple was introduced. 

The mixture was kept at —78° for 20 hours (3 in the case of propane) 
and stirred from time to time. At the end of that period, a fritted filter 
was introduced. This filter was connected with a receiver of known volume 
and a quantity of solution was passed over into the receiver. Thereafter, 
the solution in the receiver was evaporated and the receiver was detached 
and weighed, thus yielding the weight of compound in a known volume of 
solution. 
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FIGURE 1 


I. Solubility of AIBr,(CH;).0 in propane-methyl 
bromide mixtures. II. Solubility of BuyNBr in pro- 
pane-methyl bromide mixtures 


The product obtained on evaporation of the solute was shown to be the 
mono-etherate of aluminum bromide by analysis for aluminum and bro- 
mine. In the propane series, the residual product was evacuated for two 
minutes only and the results are slightly high by a few per cent in the methyl] 
bromide-rich solutions (10°, in pure methyl bromide). In the case of 
other hydrocarbons, the product was pumped for 30 minutes. With a 
solution containing 50.8 mole per cent propane, the weight obtained analyt- 
ically was 98.60, of that obtained on weighing the residue. With the other 
hydrocarbons, the error was less than one per cent. 
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In the case of the isobutane mixtures, the receiver was provided with a 
pair of electrodes and the conductance of the solutions was measured as 
they were withdrawn. 

The density of the liquids employed was determined at —S1°C. as fol 
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Solubility of AlBr,(CH;)-O0 in mixtures of met iy! bro 
mide and n-butane and isobutane. 


‘ows: methyl bromide, 1.9383; propane, 0.6265; n-butane, 0.68.0; 
isobutane, 0.6681; x-pentane, 0.7148. 

Measured values of the solubility of the dimeth, lether-aluminum bro 
mide complex in propane-methylbromide mixtures are recorded in table | 


and are shown graphically in figure |, curve I. 


TABLE 1 


SOLUBILITY OF AIBrg*(CH,)sO IN n-PROPANE-METHYL BROMIDE MIXTURES AT —78.08 


MOLE ‘; SOLUBILITY, MOLE “; SOLUBILITY, MOLE ‘, SOLUBILITY, 
PROPANE MOLES/L PROPANE MOLES/L PROPANE MOLES/L. 


00.00 0. 2668 9 96 0.1050 51.95 0.1163 
410 65 0.1958 50.76 0.1397 52.33 0.1168 
17 28 0.1396 51.23 0 1528 5411 0 0983 
$8 &5 QO. 1115 } ) 0.1244 58.56 0 0649 
19.31 QO 1021 { i 0 1288 66 68 0 O280 

10000 0 Q000 
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The solubility curve consists of two branches; the first, at lower propane 
concentrations, 1s concave toward the axis of composition, the second is 
convex toward this axis. The two branches are joined by a nearly vertical 
section at a composition in the neighborhood of 50 mole per cent propane. 
At 49.3 mole per cent propane, the solubility reaches a value of 0.105 m./1. 
It then rises sharply to a value of 0.153 m./1. at a composition of 51.2 mole 
per cent propane. The solubility branch at higher propane concentrations 
is of normal form. 

To check the experimental method, the solubility of tetrabutylammonium 
bromide was determined for the propane-methyl bromide mixtures. The 
results are shown graphically on curve II, figure 1. This curve is of the 


TABLE 2 


SoLcuBriLiry OF AlBry:(CH,;)O IN n-BUTANE-METHYL BROMIDE MIXTURES AT —78.08 


MOLE ‘% SOLUBILITY, MOLE °¢ SOLUBILITY, MOLE ° ¢ SOLUBILITY 
BUTANE MOLES/L BUTANE MOLES’ L BUTANE MOLES/1I 


0.00 0.2420 16.37 0.1794 40). 23 0.1192 
& 40 Q 2094 18.97 0.1570 41.25 0 0996 
Ripe 0 15384 22.70 0 1406 42.24 0.0742 
12.70 0.1374 30.58 0. 1323 $4.75 0.0703 
13.28 0). 2279 36.54 1292 16.80 0.0649 
13 71 0.2492 38 60 1250 65.95 0.0194 
100.00 0 0000 
TABLE 3 
SOLUBILITY AND ConpucTANCE OF AIBrs:(CHs)0 IN| ISOBUTANE-METHYL BROMIDE 
MIXTURES AT —78.08 
SPEC MOLEC SPEC MOLEC 


MOLE ‘, SOLUBILITY COND COND MOLE ‘; SOLUBILITY, COND COND 
BUTANE MOLES/L x 10 x 10% BUTANE MOLES/I x 10° x 1 


0.00 2420 28 68 118.51 14.27 0.2272 3.984 17.58 
& 00 1958 8.025 10 9S 16 02 0 2011 3.235 16.08 
10.67 15838 859 30.69 17.59 0.1720 2 269 13.19 
11 49 1255 3.783 30) 14 20 01 0.1562 

12.28 1769 280) 24.19 25. 86 0. 1408 0.578 4.11 
13.22 24S 5.159 20.98 35.82 0.1319 QO 149 1.13 


expected form and the experimental values show no scatter. The observed 
discontinuities in the solubility curve of the aluminum bromide-ether com- 
plex cannot be accounted for on the basis of experimental error. 

To eliminate the possibility of some unlooked-for interaction between 
methyl bromide and propane, the vapor pressure of mixtures of methyl 
bromide and propane was measured at —78.08° from 0.0 to 43.5 mole per 
cent propane. The curve found was a regular one which was concave to- 
ward the axis of composition and which gave no evidence of any special 
type of interaction between the components. 

The solubilities of aluminum bromide-dimethyletherate in n-butane- 
methyl bromide mixtures are presented in table 2 and those for the iso 





Vo. 88, 1955 CHEMISTRY: FOSTER AND KRAUS 241 


butane mixtures in table 3. The solubilities for both systems are shown 
graphically in figure 2, values in n-butane being represented by open circles, 
for isobutane by black circles. The solubility cur.c; for the two systems 
appear to be identical. The slight differences found at higher concentra- 
tions of hydrocarbons are probably due to experimental error. It may be 
noted that the purity of the two butanes was established by means of their 


absorption spectra. 

The solubility curve for the butane mixtures differs markedly from that of 
the propane mixtures. The solubility curve falls more rapidly with addi- 
tion of butane than with that of propane and the first break comes in the 
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FIGURE 3 


Solubility and conductance of AIBrg(CH,).O in iso 
butane-methyl bromide mixtures 


neighborhood of 10 mole per cent of hydrocarbon. The solubility then 
increases to a value of 24.9 as against 24.2 for pure methyl bromide. From 
the upper cusp of the solubility curve, the solubility begins to decrease in 
what might appear to be normal fashion but, on further addition of butane, 
the curve flattens out with a solubility in the neighborhood of 0.14. 

In the neighborhood of 40 mole per cent butane, the solubility falls 
sharply from a value of 0.119 to 0.074. Thereafter, the solubility decreases 
in normal manner. This second discontinuity does not appear in the pro- 
pane or the pentane mixtures. 

In order to obtain a check on the measured solubilities in isobutane mix- 
tures, the receptacle into which the solutions were withdrawn was provided 
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with a pair of electrodes and the resistance of the solutions was measured 
at —78.08° prior to evaporation. This provided a value of the resistance 
from which the specific conductance and the molecular conductance of the 
solutions could be computed. In table 4, the composition of the mixture 
is given in column |, the solubility of the etherate in column 2, the specific 
conductance of the solution in column 3 and the molecular conductance 
in column 4. The specific conductances are shown graphically in figure 33, 
curve I; the solubilities of the etherate are shown ascurve II. The molec- 
ular conductances are shown in figure 4. 





° 
< 
8 
rs 
3 
g 
4 
8 
fe) 
$ 














20 40 60 60 
MOLE PERCENT OF ISO-BUTANE 
FIGURE 4 


Molecular conductance of AIBr,;-(CH,;)0 in isobutane 
methyl bromide mixtures. 


As may be seen from figure 3, the curve of specific conductances closely 
parallels that of solubilities. At the lower cusp of the solubility curve, the 
specific conductance is 3.78 X 10~*, at the upper cusp it is 5.16 XK 107°; 
between the two cusps, the conductance increases in the ratio of 1.36: 1. 
If the molecular conductance of the salt in the two mixtures were the same, 
then the salt concentration would have increased in that ratio. However, 
as many be seen from figure 4, as well as from table 3, the molecular con- 
ductance at the lower cusp is 30.1 X 10~* while at the upper cusp it is 
21.0 X 107-4, a ratio of 1.43:1. If the solution at the upper cusp had the 
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same molecular conductance as at the lower cusp, the conductance ratios 
between the two cusps would be 1.95. In other words, on the basis of con- 
ductance measurements, we conclude that the concentration at the upper 
cusp is 1.95 times that at the lower cusp. Now the measured solubility 
at the lower cusp is 0.1255 and at the upper cups it is 0.2458 and the ratio 
is 1.96. Thus, the ratio of salt con- eee 

centrations at the two cusps as de- [ 
termined by conductance measure- 





iments closely checks the measured P 


solubility values. 

A careful inspection of figure 4 will 
show that two small jogs occur in the 
curve of molecular conductance at 
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Solubilities of the dimethylether- er . 

7 ; iG; Solubility of  AIBry;(CH,;).O in n 
ate of aluminum bromide in n-pen- pentane-methy! bromide mixtures 
tane-methyl bromide mixtures are 
presented in table 4 and are shown graphically in figure 5. 


rABLE 4 
SOLUBILITY OF AIBry-(CH,)c0O IN n-PENTANE-METHYL BROMIDE MIXTURES AT — 78.08 


MOLES ‘; SOLUBILITY MOLE “; SOLUBILITY MOLE ‘{ SOLUBILITY 
PENTANE MOLES (I PENTANE MOLES /1 PENTANE MOLES/I 


32.58 0. 1800 16.99 0.0509 52.06 0.0531 
410.68 0.1400 17 60 0.0651 64.46 0. O160 
43.92 0 1091 $813 0.0733 100.00 0.0000 
15.74 0. O776 1%. 72 0 O754 
16.45 0.0621 99 13 0 0686 


On addition of hydrocarbon, the solubility falls off somewhat more rapidly 
with pentane than with propane, but much less rapidly than with butane. 
The discontinuities in the solubility curve, with pentane, come at a slightly 
lower concentration of hydrocarbon (47 m. ©;) than with propane (49.3 
m.‘;). The ratio of solubilities at the lower and upper cusps is 1.5 for 
both propane and pentane, but the actual solubilities with pentane are 
much lower (one half) than with propane. The solubility curve for the 
two butanes differs strikingly from those for propane and pentane. 
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Whenever a discontinuity occurs in the solubility curve of a solid com- 
pound, a change must be taking place in the state of the saturating phase. 
In the case of butane, for example, a phase change appears to set in at the 
composition of the lower cusp, which change is completed on reaching the 
composition of the upper cusp. With the butanes, there appears to be a 
second phase change at about 41 mole per cent butane in that the solubility 
drops precipitately from 0.169 to 0.074. This second discontinuity is 
absent with the other two hydrocarbons. 

The etherate of aluminum bromide is a very stable compound and is 
always recovered without change on evaporation of its solutions. It 
would seem that, in some way, the etherate interacts with the hydrocarbon 
to form a compound much more soluble than the simple etherate. Presum- 
ably, at the lower cusp, a phase change takes place that leads to a higher 
solubility. In the case of the butanes, there appears to be a second phase 
change. 

The small solubility change on initial addition of hydrocarbon 1s surpris- 
ing and difficult to account for. The solubility of the etherate in pure 
hydrocarbons is negligibly small. On addition of hydrocarbon we should 
expect the solubility to fall off much as it does for tetrabutylammonium 
bromide or for the etherate at higher concentration of hydrocarbon. 

It would seem that the solution of this problem lies in the characteriza- 
tion of the saturating phase which is in equilibrium with the solution along 
the different branches of the solubility curve. Separation of the solid 
phase from the supernatant solution at —78° should not prove too difficult 
an operation, 

* This paper is based on a portion of a thesis presented by Luther M. Foster to the 
Graduate School of Brown University, in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy in June, 1947. Experimental details may be obtained 
by consulting this thesis which is on file in the library of Brown University 

' Morris, C. L., Thesis, Brown University (1941) 

2 Egan, C. J., and Kemp, J. D., J. Am. Chem. Soc., 60, 2097 (1938) 

’ Kennedy, R. M., Sagenkahn, M., and Aston, J. G., /bid., 63, 2267 (1941) 

' Rossini, fF. D., Chem. Eng. News, 25, 230 (1947). 
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THE POTENTIAL USEFULNESS OF NATURAL TRITIUM* 
By W. F. Lipsy 


INSTITUTE FOR NUCLEAR STUDIES AND DEPARTMENT OF CHEMISTRY, UNIVERSITY OF 
CHICAGO 


Communicated February 25, 1958 


The production of tritium by the cosmic radiation in the earth's atmos- 
phere'~® probably can now be taken as being established though many 
quantitative aspects remain unknown. Rainwater definitely contains 
tritium in amounts corresponding to between 5° 10-'* and 40 10 -'> tritium 
atoms per hydrogen atom. This is a small but definitely measurable quan- 
tity, requiring between one quart and five gallons of water as sample. 
The possible potential uses appear to be worth consideration. 

A. Age of Rain and Agricultural Products. —\n the first place it is clear 
that the tritium content of rain water will decrease with time elapsed after 
precipitation, since the cosmic ray source can introduce fresh tritium atoms 
only in the atmosphere where the new tritium atoms react with the at- 
mospheric oxygen to form radioactive water. The time scale is set by the 
decay constant for tritium 12.5 years half life and 18.0 years average life. 
We expect therefore that whereas freshly fallen rain or snow will show a 
definite radioactivity due to tritium, old water will show none, and that the 
age of the water in the sense of the time elasped since precipitation may be 
determinable from the specific radioactivity. Further, agricultural prod 
ucts should be datable by the tritium content of their constituent hydro 
gen, the principle being that a considerable fraction of the hydrogen utilized 
in plant growth originates as rain water and is incorporated in the first few 
months after precipitation. The time of harvest therefore is close to the 
time of precipitation. Consequently, it seems possible that the age of 
agricultural products in the practical sense of the time elapsed since harvest 
probably can be measured by determination of the tritium content of the 
product itself. A series of tests is in progress in our laboratory using water 
of known age and wines of known vintage to test these postulates. 

B. Iydrology and the Identification of Water Sources. A second type of 
application occurs in the identification of water sources. The mean storage 
time of rain or snow water before it appears in a given water supply should 
be given directly by the tritium content of the given supply relative to 
that of average rain in the region involved. The determination of mean 
storage times of between about three and twenty-five years should prove 


possible, these times corresponding to between 15 and 75% loss of tritium 
due to radioactive decay. The determination of tritium in wells should 
reveal whether their ultimate supplies are closely connected to the rainfall 


and therefore are subject to seasonal fluctuations in rainfall. 
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C. Meteorology.—In the field of meteorology several possible applica- 
tions seem to exist. The tritium production probably occurs at a steady 
rate independent of the weather. It seems likely that the production 
rate will vary strongly with altitude in about the same way that the produc- 
tion of neutrons and stars by the cosmic rays does. This means that on 
the average the tritium is produced at an altitude of about 6 or 7 miles, 
and that a maximum in the production rate occurs at about 9 miles. Above 
this maximum the production rate will decrease with fall in atmospheric 


pressure, and below this maximum will decrease exponentially at a rate of 


about 509% for every meter of water increase in pressure. 

1. Test of Vertical Mixing in Air Masses.— Direct sampling of the mois- 
ture in a given air mass at various altitudes should reveal whether vertical 
mixing has occurred. The determination of the extent to which vertical 
mixing occurs may be a matter of considerable interest. It seems not 
unlikely that the average air mass will be found to be mixed. 

2. The Source of the Moisture in Air Masses.—It is clear that deep 
bodies of water will so dilute rain and snow water that the tritium content 
of the waters on the surface of the seas will be very small. Therefore, 
continental water sources such as rivers will be distinguishable from the 
sea, and it should be possible to determine whether an air mass acquires 
any considerable amount of moisture from the land masses as it travels 
over them. An air mass formed over the sea will contain no radioactive 
water except for the radioactivity acquired from the atmosphere. In fact, 
for such an air mass, assuming perfect vertical mixing and neglecting the 
variation of the tritium production rate with latitude,” we can write 


o(7',) = 470 Qt, (1) 


in which p is the moisture content of the air mass (g./cm.*), 7), is the tri- 
tium content of the rain (10~'* T atoms’ H atom), Q is the tritium produc- 
tion rate (T atoms/cm.’ sec.), and ¢ is the time (measured in years) spent 
by the moisture in the air mass in the air since it left the ocean. This 
equation should quite accurately represent the tritium content of the air 
mass so long as it has not acquired moisture from continental sources. 
Therefore, measurement of the tritium content of moisture in an air mass 
along its trajectory over the sea and over the continents may reveal 
to what extent it has acquired moisture from the land. Eq. (1) affords 
a possibility of determining the residence time of sea moisture in air, 
and it seems not unlikely that these determinations may be of some in- 
terest. 

It is true, of course, that the production rate of tritium will probably 
vary with latitude rather strongly,’ the production rate at 50° north geo- 
magnetic latitude being some four times that at the equator. Conse- 
quently, the actual trajectory of the air mass must be used in the applica 





Voe. 39, 19538 CHEMISTRY: PAULING AND CORE} 247 


tion of Eq. (1), and a mean value of Q chosen according to the trajectory 


involved. 

D. The Rate of Escape of He* from the Karth..-When the assay of 
natural tritium in the waters of the earth, now in progress in this labora 
tory, has been completed, a minimum value for the rate of production of 
He* will have been obtained, because He? is the isotope formed by tritium 
during its radioactive decay. This value is only a minimum because it 
would seem likely that He* could be produced in cosmic ray stars* just as 
tritium seems to be. If we dare assume that the cosmic ray intensity 
during the last 1S years has been the same as for the last 100 million years, 
we then can calculate an upper limit for the average residence time of a He* 
atom in the earth’s atmosphere before it escapes from the earth. 


* This research was supported in whole or in part by the United States Air Force 
under Contract So. AF 18(600)-564 monitored by the Office of Scientific Research, Air 
Research and Development Command 

! Libby, W. F., Phys. Rev., 69, 671 (1946). 

% Faltings, V. and Harteck, P., Z. Naturforsch., 5a, 438 (1950). & Harteck, P. and 
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3 Grosse, A. V., Johnston, W. H., Wolfgang, R. L., and Libby, W. F., Scrence, 113, 
1 (1951). 

' Final Report to Air Research and Development Command (P. O. Box 1395, Balti 
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TWO PLEATED-SHEET CONFIGURATIONS OF POLYPEPTIDE 
CHAINS INVOLVING BOTH CIS AND TRANS AMIDE GROUPS 


By Linus PAULING AND ROBERT B. COREY 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY,* CALIFORNIA INSTITUTE OF 
TECHNOLOGY 


Communicated January 30, 1953 


In continuing our search for configurations of polypeptide chains with 
interatomic distances, bond angles, and other features as indicated by the 
known structures of crystals of amino acids, simple peptides, and related 
substances, we have investigated configurations involving both cis amide 
groups and trans amide groups, and have now discovered two pleated sheets 
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FIGURE 1 


The vertically pleated cis-trans sheet of hydrogen-bonded polypeptide chains 
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composed of polypeptide chains with cis and trans amide groups in alterna- 
tion. 

The Vertically Pleated Cis-Trans Sheet.—The configuration shown in 
figure 1 results from the formation of diagonally directed hydrogen bonds 
between adjacent polypeptide chains in which cis amide groups and trans 
amide groups alternate. The cis amide groups have their C’--N axes 
nearly parallel to the fiber axis of the structure (the direction of the polypep- 
tide chains), whereas the trans amide groups are nearly perpendicular to the 
fiber axis. The structure involves 10-membered rings containing two hy- 
drogen bonds, and also 14-membered rings containing two hydrogen bonds. 
The topological features of the structure are represented in figure 2. 


HC 


N 


FIGURE 2 


Diagram representing the topological features of both the vertically pleated cis-trans 
sheet and the horizontally pleated cis-trans sheet 

The symmetry of this sheet is orthorhombic, with vertical and horizontal 
twofold screw axes in the plane of the sheet, and twofold axes perpendicular 
to the sheet. The unit of structure, with dimensions dy = 10.32 A (perpen 
dicular to the fiber axis) and 6) = 11.38 A (parallel to the fiber axis), con 
tains four cis amide groups, which are equivalent to one another, and four 
trans amide groups, which are equivalent to one another. 

One of the atoms (other than N and C’) attached to each of the a carbon 
atoms lies in the plane of the sheet, and the other projects approximately 
perpendicular to the sheet. The first of these positions is unsuited for any 
group except the hydrogen atom, whereas the second is satisfactory for 
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, é 
any amino-acid side chain. The nature of the structure is such that a 
polypeptide composed exclusively of L amino-acid residues could have all 
of these residues in the second position, projecting out from the plane of the 
pleated sheet; this pleated-sheet configuration is accordingly, with respect 
to this consideration, a satisfactory one for proteins. 

The coordinates of the atoms are given in table |. These coordinates 
were obtained in part by calculation, and in part by the measurement of a 
set of models carefully constructed to the scale 10 cm. = | A. The as- 
sumed N--H : - - O distance is 2.77 A. The other interatomic distances 
and bond angles used are those previously reported. ! 

TABLE 1 

AroMic COORDINATES FOR THE VERTICALLY PLEATED C1is-TRANS SHEET 
dy = 10.32 A, by = 11.388 A, co = 1.00 A. Four atoms, at x, y, 2; 1/2 +4, 9.2; 8! 
y, 8; "/2~ x, '/2~ 9», 2 
: . 

C, 0.029 125 13 
0.160 O61 13 
0.264 116 13 
0.150 055 13 
0.025 117 13 
0.045 241 56 
0132 307 96 
—0. 034 263 —0.42 

TABLE 2 


Aromic COORDINATES FOR THE HORIZONTALLY PLEATED CIs-TRANS SHEET 


ay = 11.86 A, by = 10.12 A, co = 1.00 A. Four atoms, at x, y, 2; ¥, '/2 + ¥, 2; 
x, 9,3; "/2>— x, '/2—- 9, 8 


. . 
c. 0.099 —(0.147 Is 
C,’ 0.208 —() 071 92 

0.299 —0.130 76 
N; 0.194 0.058 &9 
Cy 0.085 ).122 07 

0.0385 0.236 0.24 
O. 0.117 0.295 —(0.25 
Ne —0.069 0.255 —0.20 


The existence of pleats in the vertical direction can be seen from the 
vertical projection of the structure given in figure 1. 

The Horizontally Pleated Cis-Trans Sheet.— A second pleated-sheet struc- 
ture that differs from the first one in one essential feature is shown in figure 
3; projections of the structure, also shown in figure 3, show that the pleats 
in this structure are horizontal, rather than vertical. The topological 
description of the structure is the same as for the other cis-trans sheet, and 
is given in figure 2. 
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This horizontally pleated cis-trans sheet involves cis-trans polypeptide 
chains and 10-membered and 14-membered rings, closely similar to the 
corresponding features of the vertically pleated cis-trans sheet. However 
the symmetry elements are different. The symmetry is monoclinic, rather 
than orthorhombic, and the sheet contains centers of symmetry, located in 


Oro 


FIGURE 3 


The horizontally pleated cis-trans sheet of hydrogen-bonded polypeptide chains 


the centers of the 10-membered rings, together with twofold screw axes 
parallel to the fiber axis of the sheet, and a horizontal glide plane of sym- 
metry, which converts one cis-trans chain into an adjacent chain. In this 
structure, as in the other, adjacent chains alternate in direction. 

The coordinates of the atoms are given in table 2. The value of ap (the 
identity distance perpendicular to the fiber axis) is 11.86 A and the value 
of bo (the identity distance in the direction of the fiber axis) is 10.12 A. 
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The feature of the structure that eliminates it as a possible structure for 
proteins is the presence of centers of symmetry. These symmetry elements 
require that half of the amino-acid residues have the L configuration, and 
the other half the p configuration. If it is not required that the side 
chains conform to the symmetry elements of the basic structure, and the 
side chains are introduced in the positions corresponding to the L con- 
figuration, half of the side chains are placed in positions that give rise to 
great steric hindrance, these being the positions from the a@ carbon atoms 
that are in the plane of the sheet. The horizontally pleated cis-trans 
sheet might accordingly be a possible structure for polyglycine, for a pep- 
tide im which the chains contain sequences of two L residues alternating 
with two pb residues, or for polypeptides in which pairs of glycine residues 
alternate with pairs of L (or D) residues; it is not, however, acceptable for 
proteins, except for the unlikely case of a protein with pairs of glycine 
residues arranged in this sequence with pairs of L residues. 

The Possible Occurrence of the Structures in Nature.—The only feature of 
the vertically pleated cis-trans sheet which seems unfavorable to its sta- 
bility, in comparison with other polypeptide structures, is the presence of 
the cis configuration for half of the amide groups. There is some evidence 
that the cis configuration is less stable than the trans configuration ; 
Mizushima and his collaborators’ have emphasized this point. The 
amount of instability has not yet been evaluated. The only peptide that 
is definitely known to contain a cis amide group is diketopiperazine. Evi- 
dence of instability of cis amide groups is provided by the difficulty in syn- 


thesizing the corresponding cyclic anhydrides of other amino acids. No 
x-ray evidence has so far been reported that indicates the presence of the 
vertically pleated cis-trans sheet in proteins; in particular, the identity 
distance 11.38 A along the fiber axis has not been reported for any protein. 


Acknowledgment.—This investigation was aided by grants from The 
Rockefeller Foundation and The National Institutes of Health, Public 
Health Service. 
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TWO RIPPLED-SHEET CONFIGURATIONS OF POLYPEPTIDE: 
CHAINS, AND A NOTE ABOUT THE PLEATED SHEETS 


By Linus PAULING AND ROBERT B. COREY 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY,* CALIFORNIA INSTITUTE OF 
TECHNOLOGY 


Communicated January 30, 1953 


About a vear ago, in the course of the consideration of configurations of 
polypeptide chains with favored orientations around single bonds, we de- 
scribed two pleated-sheet structures.' These structures are suited to poly 
peptide chains constructed entirely of L amino-acid residues or of D amino- 
acid residues. In one pleated sheet alternate polypeptide chains are 
antiparallel, and in the other they are parallel. The amide groups have 
the trans configuration. 

We have observed that closely related structures can be constructed in 
which polypeptide chains of p and L amino-acid residues alternate. The 
configuration of these layer structures is such as to make it appropriate to 
call them rippled sheets. 

The Antiparallel-Chain Rippled Sheet. The antiparallel-chain rippled 
sheet, represented in figure 1, is closely similar to the antiparallel-chain 


TABLE 1 


Atomic COORDINATES FOR THE ANTIPARALLEL-CHAIN RIPPLED SHEET 


ay = 9.44 A, by = 7.00 A, co = 1.00 A (assumed). Four atoms in x, y, 2; '/2 — x, ! 


y, 2; 2.9.3; "/2: + 44/2 — 9,8 
x y z 
287 0.040 0.64 
200 0.216 0.28 
O68 0. 220 0.45 


275 0.359 —().22 


pleated sheet, and the diagrammatic representation given in figure 4 of 
the earlier paper! applies to both structures. One structure is converted 
into the other by the reflection of alternate chains into their enantiomers, 
in the plane of the sheet. The unit of structure was determined by measure 
ment of a model constructed of units precisely built on the scale 10 em. = 
1A. The value of the lateral translation was found to be 9.44 A, and that 
of the translation in the direction of the polypeptide chains 7.00 A. Atomic 
coordinates are given in table |. 

The Parallel-Chain Rippled Sheet. The parallel-chain rippled sheet, 
shown in figure 2, is closely similar to the parallel-chain pleated sheet, and 
has the same diagrammatic representation, shown in figure 5 of the previous 
paper.'. The unit of structure was found by measurement of a model to 
have lateral identity distance ad) = 9.60 A, and identity distance along the 
fiber axis by = 6.50 A. Atomic coordinates are given in table 2. 
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The antiparallel-chain rippled sheet of hydrogen-bonded polypeptide chains 
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The parallel-chain rippled sheet of hydrogen-bonded polypeptide chains. 


FIGURE 2 





Vou. 389, 19538 CHEMISTRY: PAULING AND CORES 


The antiparallel-chain rippled sheet and the parallél-chain rippled sheet 
are satisfactory structures, in that they involve linear hydrogen bonds, the 
interatomic distances and bond angles have the accepted values, and the 
amide groups are planar. The orientations of the bonds around the a 
carbon atoms are such that there is room for a Side chain, projecting out 
nearly perpendicularly to the sheet, if alternate chains are constructed of 
D amino-acid residues and L amino-acid residues. These structures are 
accordingly satisfactory ones for equimolal mixtures of D polypeptides and 
L polypeptides. 

Moreover, mixtures of D polypeptides and L polypeptides in which the 
enantiomeric polypeptide chains are present in unequal numbers can assume 
sheet structures which represent a mixture of a pleated-sheet configuration 
and a rippled-sheet configuration. Irregular sequences of polypeptide 
chains with positive and negative orientations can also lead to reasonably 
satisfactory hydrogen-bonded layer structures. The value of the identity 
distance in the fiber-axis direction would be intermediate between the value 


rABLE 2 
Aromic CoORDINATES FOR THE PARALLEL-CHAIN RIPPLED SHEET 


ay = 9.60 A, by = 6.50 A, 6) = 1.00 A (assumed). Four atoms in x, y, 2; 4 
xX, "/2+ Y, 


—0 O06 0.000 
0.059 0.185 
() ISS 0 205 

—() O27 O 314 


7.00 A corresponding to linear hydrogen bonds for the antiparallel-chain 
sheets and the value 6.50 A corresponding to linear hydrogen bonds for the 
parallel-chain sheets. An intermediate value of this identity distance 
would require that all hydrogen bonds be somewhat bent, and presumably 
also somewhat strained. 

A Note on the Pleated Sheets. In the earlier discussion of the pleated 
sheet configurations of the polypeptide chains! the assumption was made 
that certain orientations around the bonds between the a carbon atom and 
adjacent atoms in the amide group are favored over other orientations. 
The assumed favored orientations around these single bonds led to the 
predicted value 6.68 A for the identity distance in the fiber-axis direction 
(the direction of the polypeptide chains). Although the existence of a 


potential function causing certain orientations around these bonds to be 
favored is likely from @ priori considerations, the magnitude of the effect 
is uncertain, and it may well be that very little strain (less than 0.1 keal. 


mole”! per residue) is involved in rotating the chain from favored to less 


favored configurations. 
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We have now found through the construction of large-scale models that 
other values of the identity distance in the direction of the polypeptide 
chains are indicated for the two pleated sheets, corresponding to slightly 
different orientations about the single bonds to the a carbon atom. The 
new configurations, described by the coordinates in tables 3 and 4, involve 
linear hydrogen bonds that is, the angle N--H---O is equal to 180°. 
The fiber-axis identity distance for the antiparallel-chain pleated sheet is 
7.00 A, and that for the parallel-chain pleated sheet is 6.50 A; these values 
are the same as for the corresponding rippled sheets. 

The two pleated sheets provide satisfactory structures for proteins and 
for polypeptides composed exclusively of L amino-acid residues (or D amino- 
acid residues). We think that it is likely that silk fibroin, for which the 

TABLE 3 
Atomic COORDINATES FOR THE ANTIPARALLEL-CHAIN PLEATED SHEET 
ay = 9.50 A, bo = 7.00 A, co = 1.00 A (assumed). Four atoms in x, y, 2; ¥, '/2 
I /, 


I/y — x, 9, 2; y, 3 


x 
0.0384 O05 
—0.030 173 
0.051 320 0.21 
0.180 326 0.22 
0.024 —(0.005 —2.24 
PABLE 4 
Aromic COORDINATES FOR THE PARALLEL-CHAIN PLEATED SHEE1 
= 4.85 A, by = 6.50 A, co = 1.00 A (assumed). Two atoms in x, y, 2; *, '/2 + y, 3 
x v Zz 
0.012 0.000 —O0.98 
—0. 066 0.186 —0.26 
0.118 0.315 0.28 
0.371 0.295 0.25 


BC —0).093 0.014 —2.45 


observed fiber-axis identity distance is 7.0 A, has the antiparallel-chain 
pleated-sheet structure, and that the 8-keratin proteins, for which the ob- 
served fiber-axis identity distance is about 6.6 A, have the parallel-chain 
pleated-sheet structure, or a structure in which a considerable number of 


adjacent polypeptide chains have parallel orientations. 
The atomic coordinates given in tables 3 and 4 differ only slightly from 
those previously reported. ' 


Acknowledgments. — ‘This investigation was aided by grants from The 
Rockefeller Foundation and The National Institutes of Health, Public 
Health Service. 
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ON THE INCLINATION OF PLASTIC SLIP BANDS IN FLAT BARS 
IN TENSION TESTS 
By T. Y. THoMAS 
GRADUATE INSTITUTE FOR APPLIED MATHEMATICS, INDIANA UNIVERSITY 
Communicated January 21, 1953 


1. Introduction. In the following communication! we have applied 
the Hencky stress-strain relations for small plastic displacements, under the 
assumption of incompressibility,? to the problem of determining the in- 
clination of the plastic slip bands which develop in flat bers subject to ten- 
sion. This elastic-plastic problem is here shown to admit a simple exact 
solution which in a sense is analogous to the simple solution of the uni- 


formly stressed bar in classical elasticity theory. 


As the basic equations of this theory we have*® 
Sap = Nap (stress-strain relations), 


u = 0) (equation of incompressibility), 


a,@ 
Copp = O (equilibrium equations), 


where the proportionality factor g, which is assumed to be positive, is one 
of the initial unknowns of the system. The o,, are the components of the 
stress tensor and u, the components of the plastic displacement. The 
tensors s and » are the stress and plastic strain deviations defined by 


Sap = Cap 3 Poa; Nap = Cap + Yass 


1 ies l — O41 ‘ — Ci; 
where - (tle pt Ug a) P = —1¢g= 


‘ 
~ o . 


The tensor e is the plastic strain tensor and is here identical with the strain 
deviation n since gq = 0 on account of (2). These relations are referred to 
a rectangular coordinate system. The equations (1), (2), and (3) are 
+ 1+ 3 = 9in number and involve the 10 dependent variables o,,4, “4, 
and g. Another equation is provided by the yield condition of which two 
forms are commonly considered, namely 


SapSag = 2k? (von Mises yield), (6) 
51 — $3 = 2k (Tresea yield), (7) 
where k is a positive constant depending on the material. In the condition 
(7) the quantities s; and s; are the greatest and least of the principal values 
of the stress deviation. 


Along the surface of separation of an elastic and plastic region the follow- 
ing general dynamical conditions of equilibrium must be satisfied 
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pilv, G) = pol Ve, G), 


lo,3|v” pPil@y, G) lvl, 


where, as usual, the bracket denotes the discontinuity in the quantity en 
closed. Furthermore v stands for the velocity of the flow and G is the 


normal velocity of the surface of separation whose unit normal vector is 
denoted by v. The quantities p; and 7, are the density and normal velocity 
on one side of the surface of separation while p, and v,, are the correspond 
ing quantities on the other side of this surface. Since we are here con 
cerned with an equilibrium theory of elastic and plastic displacements we 
will have 7, Qand G = 0. Hence the first of the above conditions ts 
satisfied identically and the others reduce simply to 


[oaslv” = 0. (S) 


The conditions (S) can be expressed by saying that the stress vector, 1.e., 
the vector having the components o,,v”, is continuous across the surface 
of separation, Other conditions which must be taken into account in this 
problem involve the vanishing of the stress vector on the outer boundaries 
of the plastic band. 

Both under the von Mises and the Tresea yield conditions our solutions 
lead to an inclination angle of #£45°16' for the slip bands. However in 
the case of the Tresea yield an inclination of +50°46' appears also as a 
possibility. Now the angle of 35°16" has been deduced by methods which 
seem to be more or less established in the literature and which have been 
presented in the book by A. Nadat.° My objection to these methods is the 
following. Apart from the fact that no plastic solution is exhibited, the 
conditions (8), which must be satisfied over the slip surface, as well as the 
outer boundary conditions are neglected. As a consequence it has been 
necessary to introduce an assumption, extraneous to the basic equations of 
the theory, which we tind difficult to understand.’ It may be mentioned, 
however, that such disregard for the basic requirements of theory appears 
to be a characteristic of the work of many of the professional writers in 
this area. 

2. The Ideal Flat Bar. Yo say that a rectangular bar is thin or flat 
wuplies that its thickness is small in comparison with its other dimensions, 
i.e., its width and length. Let us therefore define the ideal flat bar as one 
of finite thickness but whose flat sides are formed by two infinite parallel 
planes. By employing the concept of the ideal flat bar we can limit our 
attention to the boundary conditions on the two flat sides. Results ob- 
tained for such bars may be expected to be valid approximately in the 
interior, 1.e., in a region removed from the narrow sides and the ends, of 
the flat bar of finite dimensions. 

In figure | we have indicated a portion of an ideal flat bar. The lines 
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AB and CD are the intersections, with one of the flat sides of the bar, of 
parallel planes P; and P, perpendicular to the flat sides. Two coordinate 
systems x and y will be used in the following discussion. Both of these 
systems will be supposed, for definifeness, to have their origins at a point 
O mid-way between the flat sides of the bar and also mid-way between the 
Hanes P; and Py. The x, and y; axes (not shown in Fig. |) will be taken to 


Tt 





FIGURE 2 


tT 
FIGURE | 


be identical and perpendicular to the flat sides. The y» axis 1s chosen per 
pendicular to the planes P; and P» and hence the y; axis will be parallel to 
these planes. Since the x), .. and yy, Yo coordinate axes are in the same 
plane it can be assumed that the y,, ye axes are Obtained from the x, x» 
axes by rotation through an angle @ (see Fig. |). Hence the coordinates 


of these two systems will be related by the following transformation 


Y¥, = y, cos 0 yo SIN OA; Xe y, sin 6 + yo COS A; NX, Vs. (DY) 


Now suppose that the ideal flat bar is subjected to a uniform tension 
denotes the components 


7(>O) 1n the direction of the x. axis. Then if 7,4 
of the stress tensor relative to the x system we shall have r.. = 
Denoting the components of the stress tensor by u,, 1n the y 


rand rT. 


0 otherwise. 
system the quantities uw,, can be determined from the values of the 7,4 by 


the tensor transformation relating these two sets of components. ‘Thus 


Ox, OX Ox» OX 


Hines Ov, OV; Oy, OV; 


Computing the derivatives in these relations from (9) we tind that 
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Wun = 751N? 6; wy = 7 Sim Ocos 6; wy, = O, (10) 
foo = 7 COS? 8; wos = 0; wa, = O. 


o 


3. The Plastic Sip Band. Suppose now that the value of 7 is allowed 
to increase until the yield value is reached. It will then be assumed that 
the band or layer between the two parallel planes P; and P, becomes plastic 
in the sense that the displacements uv, in this region satisfy the equations 


(1), (2), and (3) and one or the other of the yield conditions (6) or (7). 
In the case of the von Mises yield we shall have at the yield point +r = 1/3k 
and for the Tresea yield we shall have r = 2k. These relations between 
rand & are an immediate consequence of the von Mises or Tresea yield 


condition and the values 7, of the stress components given in §2. Outside 
this plastic band it is assumed that the material is in the elastic state. 


Remark.--\t is assumed that the plastic displacements ug under consideration are 
relative to the natural or undetormed state. We therefore conceive of the material be- 
tween the above planes ?; and 7» as in the undeformed state before determining its 
displacement uq as a solution of the equations of $1. The possibility of realizing this 
condition can be seen as follows. After the tension in the bar has been brought to its 
yield value we transform the elastically deformed material between the planes P; and P, 
by the transformation inverse to that producing the elastic deformation. By this pro 
cedure the planes ?; and ?, will be transformed into parallel planes P?;’ and P:’ which 
will be perpendicular to the flat sides of the bar and the material between these planes 
will be in the undeformed state. A rigid displacement or rotation ot the undeformed 
layer between the planes ?;’ and P.' will now bring the plane P,’ into contact with the 
elastically deformed material which was previously along the plane ?;. Similarly a 
translation (rigid displacement) of the elastic material formerly in contact with ?» will 
bring this material into contact with the plane P.’. In the following we have omitted 
the primes and have denoted by ?,; and P, the planes separating the elastic material and 
the band which may be assumed to be placed symmetrically relative to the elastic por 
tions of the bar 

After the plastic displacements ug have been determined rigid displacements of the 
elastic portions of the bar may be allowed as well as rigid displacements of the plastic 
layer bounded by the planes P; and P, provided that these displacements do not pro 
duce a separation of material. In fact the elastic displacements and plastic displace 
ments W under consideration can be determined at most to within such rigid displace 
ments by the conditions of the problem. 


We seek a simple solution uw, of the equations (1), (2), and (3) in the 
plastic band satisfying the boundary conditions, 1.e., for which the stress 
vector vanishes (after deformation) on the free surface of the band, and 
such that the conditions (8) are satisfied over the surfaces of separation of 
the elastic and plastic regions. It will be shown in fact that, relative to the 
y coordinate system of $2, such a solution is given by 


Uy = AV2; U2 = bye; uz = —by;, (11) 


where a and } are suitably chosen constants. In the first place it may be 
observed immediately that, regardless of the values of the constants a and 
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b, the equation of incompressibility (2) is satisfied by the displacements 
(11). 

In consequence of (11) the planes P; and P, bounding the plastic band 
will be displaced into planes parallel to P; and P». This is illustrated in 
figure 2 where we have indicated the displacement, in a plane parallel to 
the flat side of the bar, of two points Q and R in the planes P; and P., 
respectively. This part of the displacement, which is determined by the 
first two relations (11), involves a shear followed by a movement in the y» 
direction proportional to the y2 coordinate. From the third equation (11) 
it is seen that planes parallel to the flat sides of the bar remain parallel to 
the flat sides. 

Since the displacement within the band, as above described, envolves a 
slip of material particles along the planes of separation P, and P», it is 
natural to refer to this band as a plastic slip band. In consequence of the 
displacement (11) within the plastic slip band the portions of elastic ma- 
terial on either side of the band will suffer a rigid displacement which will 
not effect their internal stresses. It is important to observe that these dis- 
placements can occur without separation of the material along the planes 
P, and P%. 

4. Calculation of the Plastic Tensors. We shall now calculate the com- 
ponents of the stress tensor o and its deviation s within the plastic slip 
band. This calculation will be carried out relative to the y coordinate sys- 
tem, Observe first that by differentiation of the equations (11) we have 

W,1= 0; Mm,2= a; M,3 = OV, 
U2, 1 0; uU2,2 = 


W3,1 = O; U3,2 = 


From these relations we immediately 
components é,3 from (5). Thus 


én = OF ey 


x2 = bd; ex 


Hence from (1) and the fact that n,, = e,3 for the incompressible case 


under consideration we can write 


| (12) 


Sx = bg; $93 = 0; 533 = sia 


Finally from (12) and (4) we obtain 
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ag 


“2 Oi 


Ov bg p; a 


3. The Boundary Conditions. \t follows from the third equation (11) 
that the quantity — is the plastic elongation in the y, direction. We 
choose this quantity to be equal to the corresponding elastic elongation. ‘This 
amounts to taking the constant 6 equal to vr / where v is Poisson's ratio 
and /: is Young's modulus for the material. Since the value of } so de- 
fined is positive we see that the distance between the planes P; and P» of 
§2 will be increased by a plastic deformation (11) resulting from tension. 

With the above value of the constant ) the displacement in the y; direc 
tion of the sides of the bar, including the outer boundaries of the plastic 
band, will be determined entirely by the elastic deformation. The outer 
boundaries of the plastic band can therefore be made to lie, by the deforma- 
tion (11), in the planes of the flat sides of the elastic portions of the bar. 
Hence the plastic boundary condition is given by o,,6° = U where & is a 
unit vector in the y; direction. In other words we must have o,; = 0 over 
the two outer boundaries of the plastic band. But from (13) we see that 
these equations are satisfied fora = 1,2. Taking a = 3 it therefore fol 
lows that 


p be (14) 


over the outer boundaries of the plastic band on account of the last equa 
tion (13). Now it is immediately seen from (12) and either of the yield 
conditions (6) or (7) that the quantity ¢ must be a constant in the plastic 
band. We choose the quantity p to be the constant given by (14). The actual 
value of p will not be needed in this discussion but it can be found from (14) 
and the value of the constant ¢ which is determined by the relations of §7 
$8. With this choice of p the boundary conditions are thus automatically 
satisfied. .\/oreover we see from (13) that the stress components a, are con 
stant throughout the plastic band and hence the equilibrium equations (3) are 
identically satisfied. 

6. Conditions Along the Slip Planes. Let us now consider the cond1- 
tions (8) which must be satisfied along the two planes separating the plastic 
band from the elastic portions of the bar. These will be the planes P; and 
P, of §2 or more precisely the planes into which P; and P» are displaced by 
the plastic deformation (11). In either case the planes of separation are 
perpendicular to the y. axis so that the vector vin (S) has the components 
(O, 1,0). Hence (S) reduces to 


lo.) = He Oa? 0, (15) 


where the values of the w and o are given by (10) and (15), respectively. 
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Making these substitutions and taking account of (14) the relations (15) 
give 

; a¢g 

rT sin @ cos @ = 


» 


7 cos’ 6 = he 
Dividing corresponding members of (16) and (17) we now bave 


a 
tan @ = ‘ (1S) 

1b 
This equation will determine the inclination @ of the slip band after the 
value of the constant a has been found. To determine a we first set up the 
relation obtained by the elimination of @ from the above equations, This 
elimination is accomplished quite simply by substituting for tan’ @ and 
sec? (= 1 /cos* @) in the identity 


1 + tan’ 0 = sec’ 6 


the values of these quantities given by (17) and (1S). Thus 


a T a’ + 16d" T 
= ’ or (19) 
16b° Phe Sb ¢. 


So far the yield conditions (6) and (7) have not been used explicitly. In 
the following sections the ratio + ¢ is found for each of the two yield con 
ditions. Then the equation (19) will determine the value of the constant 
a and hence from (18) we can find the inclination @ of the slip band. 

7. Von Mises Vteld. When we substitute the values of s,, given by 
(12) into the condition (6) the following expression for ¢ is immediately ob 
tained 


2k ZT 


; = ’ (20) 
Va’ + 4b? V3 Va*> + 40° 


where the plus sign is taken before the radical since it is assumed that ¢ is 
positive (see §1). Substituting the value of ¢ given by (20) into (19) we 
now find 
a? + 166? = 473 bVa? + 46°. (21) 

Squaring both members of (21) we obtain a fourth degree equation in the 
quantities a and 6 which is seen to reduce immediately to the simple rela 
tion a = ++/Sb by which the constant a is determined. Hence equation 
(18) becomes tan #@ = +1°7/2. /t follows that the slip bands have an in 
clination of +35°16' under the von Mises yield condition. 

8. Tresca Yield. The principal values of the stress deviation are given 
by the solutions » of the determinantal equation 
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” 


(bg 


0 ) (—b¢e — s) 


Expanding this determinant and solving for s we find 


b l 
s= -—bg,s = (; + 5 Va + ) ¢. 


Let us denote these three values of s by 5), s2, and s, with the understanding 
that s; 2 s» 2 s;. Assume for the moment that 


b | - 
» = —be; 83 = (| is Va? + ») ¢. 


Then the Tresca yield condition (7) gives 
5 — 8% = Va? + b? ¢ = R=. (23) 
Substituting the value of the ratio r/¢ given by (23) into (19) we have 
a? + 16b2 = 8bVa? + b%. 24) 


Now when we square both members of (24) we find that the resulting 
fourth degree equation in a and 6 leads to the following values of the con- 
stant a, namely 


a=z+v78)b; a= +7245. (25) 
Substituting these values of a in turn into the equations (22) we find 


5, = 2g; 8. = —be; 8s; = —bd¢, 


5, = 36¢; s% = —be; s; = —2be. 


Hence the required inequalities s; 2 s» 2 s; are satisfied, since b and ¢ are 
both positive, for each of the two possible sets of values of the constant a 
given by (25). 

Now it is evident that s; as given by (22) is the greatest of the principal 
values of the tensor s. If we suppose that the least of the principal values 
is given by s; = —bg and proceed as above it will be found that only the 
first set of relations in (25) is obtained. Hence (25) gives all relations pos- 
sible between the constants a and 6. 
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When we substitute the values of a given by (25) into the equation (18) 
we now have 
l V3 
tan@d = + os ; 
V2 V2 
Hence 6 = +35°16' or @ = £50°46'. In other words the plastic slip band 
has an inclination of +35°16' or +50°46' under the yield condition of 
Tresca. 


' Prepared for the Naval Research Laboratory, Washington, D.C 

*Cp. Prager and Hodge, The Theory of Perfectly Plastic Solids, Wiley & Sons, New 
York (1951), p. 36. Usually the equation of the classical elasticity theory relating the 
divergence of the displacement vector and the mean stress is assumed in this theory in 
stead of the equation of incompressibility (2). See, for example, Hill, R., The Mathe 
matical Theory of Plasticity, Oxford (1950), p. 46, and various articles by Russian writers 
in Plastic Deformation, Mapleton House, Brooklyn, N. Y. (1948). This modification of 
the above system of equations (1), (2), and (3) allows for the effects of compressibility 
in the determination of the inclination of the slip bands. We have also treated the com 
pressible case and have found that under compressibility the inclination of the slip binds 
depends in a very simple manner on Poisson's ratio for the material. This result ap 
pears in these PROCEEDINGS, 39, 266-273 (1953) 


* Equations (1) are sometimes written in the form 


l + @D 
Nap = | 9G Saas 
a(r 


and when so written they become identical with the corresponding elasticity relations 
for ¢ = 0. However it has seemed advantageous to write these equations in the form 
(1) for the requirements of our discussion 

4 See a forthcoming paper by Thomas, T. Y., “Singular Surfaces and Flow Lines in the 
Theory of Plasticity,”’ to appear in the April issue of the Journal of Rational Mechanics 
and Analysis (19538). This paper contains a rather extended discussion of the flow equa 
tions of von Mises in connection with the occurrence of flow markings resulting from slip 
as well as those which result from a weaker form of discontinuity. Since the concept 
of velocity is involved in this discussion it was necessary to use the equations of motion, 
containing non-linear inertial terms, instead of the simpler equations of equilibrium 
(3). By making certain approximation assumptions the angle of 35°16 was found to 
be the angle of inclination of slip lines on flat bars under tension although an exact solu 
tion of this flow problem, corresponding to the result of the present paper, was not ob 
tained. 

5 Nadai, A., The Theory of Flow and Fracture of Solids, McGraw-Hill, New York, 
1950, pp. 820-327 

® It is stated that e, = Gin the direction parallel to the edges of the oblique layer when 
it starts to neck down because the adjoining portions of non-yielding material do not 
participate in the deformation. In consequence of this assumption it follows that o, = 
a,/2 and by means of this relation the inclination of the slip bands is found. See Nadai, 
loc. cit., p. 323. 
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THE EFFECT OF COMPRESSIBILITY ON FHE INCLINATION 
OF PLASTIC SLIP BANDS 1.. FLAT BARS 


By T. Y. THomMaAs 
GRADUATE INSTITUTE FOR APPLIED MATHEMATICS, INDIANA UNIVERSITS 
Communicated February 2, 19538 


1. Preliminary Remarks. The condition of incompressibility used in 
our preceding note! in these PROCEEDINGS will now be replaced by the con 
dition of the classical elasticity theory relating the divergence of the dis- 
placement vector and the mean stress. This leads to the usual form of 
the Hencky stress-strain relations for small plastic displacements which, for 
purposes of reference, we shall write as follows 


feu, 


Nap, 7 . : 
1 — 2p 


(1) 


where /¢ is Young's modulus and v is Poisson's ratio for the material. ‘To 
the stress-strain relations (1) we must add the equilibrium equations and 
either the von Mises or the Tresca yield condition. The notation em 
ployed in the above and following equations is the same as that in our pre 
vious note to which reference should be made. 

The above system of equations will be applied to the problem of deter 
muning the inclination of plastic slip bands in an ideal flat bar in a state of 
simple tension.” The discussion in $2 and $3 of the previous note, involv 
ing the selection of the v and y coordinate systems and the planes P, and P, 
bounding the plastic slip band, is applicable and will be understood to ap 
ply in the following treatment. It will be shown that for suitable values 


of the constants a, 6, and ¢ the equations 
uy AV2; Uo = bye; U3 = CVs, (2) 


will give an exact solution of the problem. This statement unplies that 
the boundary conditions for the plastic band will be satisfied as well as the 


required conditions along the planes P,; and P», separating the plastic band 


and the elastic portions of the bar. 

In our previous note we obtained an angle of +45°16' for the inclination 
of the slip bands, both under the von Mises and the Tresea yield conditions, 
with the possibility of an additional inclination of +50°46" under the 
Tresea condition. The theoretically unacceptable attempts in the current 
literature, discussed briefly in $1 of the previous note, to solve this problem, 
have led also to an angle of 35°16/ for the inclination of the slip bands in 
the bar. Now these results are obtained under the hypothesis of incom 
pressibility of the plastic material. The results of the present article, giv 
ing the inclination of the slip bands as a function of Poisson's ratio, depend 
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¢ 


clearly on the effect of compressibility which is allowed by’ the stress-strain 


relations (1). In $4 and $5 we have deduced the inclination formulas 
under the von Mises and Tresea yield conditions, respectively. Formula 
(17), which is valid both for von Mises and Tresea yield, would appear at 
first sight to represent the most likely relationship between Poisson’s ratio 
and the inclination ef the slip bands. For a Poisson ratio of '/; this for 
mula gives an inclination angle of + 30° which is the angle we have personally 
observed on specimens of aluminum alloy tested at the Naval Research 
Laboratory, Washington, D. C. It would be interesting to have exact 
experimental data for comparison with the theoretical relationship between 
Poisson's ratio v and the inclination @ of the slip bands' contained in the 
table in $4. 

2. Formal Calculations. The components of certain of the plastic 
tensors will be needed and these will be calculated relative to the y co 
ordinate system. By combining the quantities u, , resulting from differ 
entiation of (2) we obtain the stress components é,3 as follows 


ei 


C29 = b: 


Using these values of the e,g the strain deviation components nag can 
readily be found and hence, from the first of the equations (1), we obtain 
the components s,, of the stress deviation. Thus 


b+c 
su = = ¥¢;, Jv 
(> 2 - 


From the formula relating the components of stress o,, and stress deviation 
Sag We can now find the o,, in terms of the s,, and the mean stress. But 
this latter quantity can be eliminated by the last equation (1). Hence we 
obtain 


b + ¢ b+ 0k 2b b+ok 
ah = ¢+ 3(1 : On - , 


5 3 = 2y) 5 


3(1 — 2p) 
h + 0k 2 
¢g + : OW ree O13 do, = O 


3(1 — 2p) 


(4) 


6? 
” 


> 


3. Boundary and Separation Conditions. The boundary conditions for 
the plastic band require the vanishing of the quantities Case’ over the free 
boundary of the band (after plastic deformation) where £ is a vector normal 
to this boundary. But by choosing the constant ¢ in (2) equal to — rv F, 
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where 7(>0) is the yield value of the tension, the free boundaries of the 
plastic band will lie in the planes of the elastically deformed flat sides of 
the bar and hence the vector & will be parallel everywhere to the ys axis. 
This leads to the conditions ¢,; = 0 which are satisfied automatically for 
a = |, 2 on account of (4). The remaining condition becomes 


b+ok < 


(1 — 2p) 


(2c — b)g + 0. (5) 


Now the quantity ¢ must be constant throughout the band in consequence 
of the yield condition. Hence the condition (5) can be viewed as a neces- 
sary relation between constants and not merely as a condition restricted 
to the outer boundaries of the plastic band. It may be mentioned also in 
this connection that the stress components o,, must be constant throughout 
the plastic band and hence the equilibrium equations will be satisfied iden- 
tically. 

We now consider the conditions which must be satisfied along the planes 
P, and P» separating the plastic band from the elastic portions of the bar. 
These conditions are expressed by the vanishing of the stress vector over 
the planes P; and P,. In other words we must have 


p \ 
[ouglY = (Hap il Cap)¥ = 0, 


where the v’ are the components of the unit vector normal to the plane P; 
or P, and the wy, are the stress components in the elastic part of the bar 
relative to the y coordinate system. The quantities u.g were calculated 
in our previous note and are given by 


Mu = T sin? 6; Mig >= T sin # cos 6; M3 = 0, 
Me = 7 COS? A; pos = OF was = O, 


where 7 is the yield value of the tension and @ the inclination of the slip 
band. But v® has the values (0, 1, 0) and hence the above condition be- 
comes fa2 = o These relations are satisfied automatically for a = 3. 


an 


For a = 1, 2 they give 


: ag 
7 sin 6 cos @ = ? 


7 cos? 6 = (b — c)¢, 
when use is made of the relation (5). 
It will be found helpful to replace equation (6) by the relation resulting 


from (6) and (7) by elimination of the angle 6. To effect this elimination we 
may first set up the auxiliary equation 
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a ) 
tan @ = , (S) 
2(b — c) 
obtained by dividing corresponding members of the two equations (6) and 
(7). Combining (7) and (8) in an obvious manner we can now obtain 


a” + 4(b — c)? = 4(b — c) -: (9) 


Let us also eliminate / from (5) by means of the relation / = —rv/c used 
to define the constant c at the beginning of this section. This gives the 
following equivalent equation 


v(b + c) = = (1 — 2v)c(2c — Db). (10) 
¢ 


Equations (7), (9), and (10) are three equations for the determination 
of the four quantities a, 6, #, and ¢. The remaining equation is furnished 
by the von Mises or the Tresca yield condition. 

4. von Mises Yield.-Substituting the s,, given by (3) into the von 
Mises yield condition we are led, after some reduction, to the following 
equation 


3a" 
(0° +c? — be + : Jes = 3k? = 17, (11) 


We now eliminate a? from (11) by means of (9). This gives a quadratic 
equation in the ratio +/g which can be written 


(r/¢)? — 3(b — c)r/¢ + 2b? + 2c? — H5be = 0. 
Solving this equation for the ratio +/¢g and labeling one of the solutions 


Case I, the other Case II, we have 


4 T ‘é 
Case I, - = b — 2¢, (12) 
¢ 


C= hm <. (13) 
¢ 


Each of these cases must be treated separately. 
Case I. Substitute the value of the ratio r/¢ given by (12) into (10); 
the resulting equation reduces to 


b = ( - "Ye (14) 


by which the constant } is determined. Hence 
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from (12) and (14). Also from (14), (15), and (9) we obtain 
a’? = 4c*/y, ora = +2c Vv. (16) 
Finally the inclination @ can be determined by substituting the values of 
6 and the ratio rg into the equation (7). This gives 
a l l—yp = 
cos? 6 = , or cos 26 = (17) 
by 1+ p 
Case I]. When we substitute the value of the ratio 7/¢ given by (13) 
into (10) we obtained a quadratic equation in 6. Solving this equation we 


find 


h -(1 — pv) | 
: t Vi+ 14y — 2353p’. (18) 
r by fy 
With b determined by (1S) we obtain the equation determining the con- 
stant a by eliminating the ratio r/¢ from (9) and (13). This gives 
a’ = 4b(b — c). (19) 
Hence from (7), (13), and (1S) we have 


' —(1 — ») V1 + 14y — 239? — 4 
cos’ @ = | ; + P : af ; (20) 
—2(1 9) & 2 1 Pp l4e — 239" — 4p 
The expressions for the determination of the angle @ can be simplified 
if we replace (20) by the equivalent equations for cos 26. We thus obtain 


2p 
cos 26 = j : 
lL+pvt VI + 1l4pv — 23r? 
To consider the relations (21) more closely let us write 


2p 
cos 20’ = (0 <p 
lL+vu+t+ V1 + l4y — 237° 


Dp l 
cos 26, = 0O<vs » 28) 
lL+p— V1 + 14y — 239? 1 


where we have used 6’ and 6, instead of @ to distinguish between these cases. 
It is easily seen that the equation (22) is valid for the entire range of per- 
missible values of Poisson's ratio vy, as indicated, while in the relation (23) 


we must limit our attention to values of » not exceeding '/,. The relation 
(22) corresponds to the selection of the minus sign and (23) to the plus 
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sign before the radical in (18), Also it is easily seen that if @ is determined 
by (22) or (23), and the corresponding value of 6 taken from (18), the value 
of a®, given by (19), will not be negative. However the value of the con- 
stant ¢ given by (13) will be positive only for angles 6 determined by the 
formula (22). Moreover the value of ) will be negative for slip bands 
whose angle of inclination 6; is determined by (23) and hence there will be a 
contraction, in the y. direction, of the band from its undeformed state. 
Such slip bands would not appear to be realized in the case of actual ma- 
terials. 

It is seen from equation (17) that, corresponding to any value of v, there 
will be two slip bands whose inclination angles @ are equal except for alge- 
braic sign. A similar remark applies to the inclination angles @’ which are 


TABLE 1 
» 6’ 
OOO 15°O' 
O25 8 14°21" 
050 : $3 46" 
O75 43 
100 3o° $2 
125 28’ $2 
150 ; 41 
175 22 4] 
200 24 °6' 41 
225 25°23" 
250 34’ 


ore 


aie 


§2°5’ 

300 2 37 < 51°55 
325 41° : 51 °46' 
350 30°37’ 33 51°37’ 
375 31°29’ 38 °2 51°29’ 
400 32°19’ 37 51°20’ 
425 33 °6’ ote SLI’ 
150 33°51’ 36 51°3’ 

$75 ‘ 34’ oe 5OPSS’ 
500 35°16" 35°16 5OP46" 


determined from the equation (22). The principal or first quadrant values 
of the angles @ and 6’ have been listed in the second and third columns of 
table 1. 


Jd. Tresca Yield. The principal values of the stress deviation tensor, 
whose components s,, are given by (3), are readily found to be 


b — 2 | / 26 ~ 0 
s; = F = que" + b?hy; so = ¢; 
) “ . oO 


b— 2¢ | / 
3 = — —-Vqr+ bt}y. (24) 
6 2 
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Let us assume that these values are ordered so that 5; 2 se 2 53; this 
assumption will later be verified. Then the Tresca yield condition gives 


S; — 53 = Va? + hb? ¢ = 2k = 7, 25) 


where r is the yield value of the tension. Eliminating the ratio + ¢ from 
(9) by means of (25) we obtain a quadratic equation in (6 — c), namely 


1(bh — oc)? —~4Va° 4+ (6b — co) +a? = 0. 


When this equation is solved for (b — c) in the usual manner the two rela- 
tionships obtained can be written as follows 


Case I, Va’* + b? b — 2, (26) 
Case IT, Va? + h* 3b — Le. (27) 


We shall examine each of these cases in turn. 

Case I, On account of (25) and (26) the equation (10) reduces to (14) 
thus determining the constant >. Eliminating 6 from (10) by means of 
(14) the equation (15) for the ratio r/¢g is again obtained. Likewise we 
find (16) for the determination of the constant a when we make the sub- 
stitutions (14) and (15) into the equation (9). Also, as before, we are led 
to the relation (17) for the determination of the inclination @ of the slip 
band. For this solution s; > s. = s; so that the required identities in the 
principal values of the stress deviation are satisfied. It is interesting to 
find that the results for this case are the same as those obtained for Case I 
under the von Mises yield condition. 

Case II. From (25), (27), and (10) we find 


2(1 — v)ce? — (1 — v)be — 3vb? = 0. 


Solving this quadratic equation in c* we obtain 


bh b | 9: 
alt veer 
to4 ow 


C= (2S) 


The constant 6 is determined by (28). To determine a we have from (25), 
(27), and (9) that 

a? = 4(b — c)(26 — c), or = +)9V(b —c)(2b —c). (29) 
Hence from (7), (25), and (27) we have 


- —b 
cos? 6 = cos 26 = - (30) 
. > > Bb hoe IC 
Eliminating the constants 6 and ¢ from the second of the relations (30) 
by means of the substitution (28) we find 





VoL. 39, 1953 ENGINEERING: T. Y. THOMAS 


(0 <p 
I + PSv 


i | “ae 


where the range of permissible values of vy has been indicated and we have 
replaced @ by 6” and @, to distinguish between these two cases. It is to be 
observed that the relation (31) corresponds to the selection of the minus 


sign and the relation (32) to the selection of the plus sign before the radical 
in (28). In the case of the relation (52) neither the constant > nor the 
quantity in the right member of (27) will be positive as can easily be seen. 


Hence the possibility of having a slip band whose inclination 6 is given by 
(32) must be discarded. No such difficulty arises in the case of the rela- 
tion (31); in fact for angles 6” determined by (31) the constants 6 and ¢ 
will be positive as will be also the value of a° given by (29). Finally it can 
readily be observed that s; > s. > s, so that the relation (25) gives the cor- 
rect form of the Tresea yield condition for slip bands whose inclination 6” 
is determined by (31). 

The entries in the last column in table | are found from the relation (31) 
and are analogous to the entries in the second and third columns of this 
table. 

' Thomas, T. Y., these PROCEEDINGS, 39, 257—265 (1953 ) 

* Prepared for the Naval Research Laboratory, Washington, D.C 

* See Nadai, A., The Theory of Flow and Fracture of Solids, McGraw-Hill, New York 
(1950), pp. 820-827. 

* Actual slip bands are not to be confused with the bands or markings produced by 
singular surfaces, e.g., characteristic surfaces, which represent a weaker form of discon 
tinuity. See Thomas, T. Y., “Singular Surfaces and Flow Lines in the Theory of Plas 
ticity,’ J. Rational Mech. and Anal., April (1953) 
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NONLINEAR MULTIPOLES 
By L. A. ZADEH 
DEPARTMENT OF ELECTRICAL ENGINEERING, COLUMBIA UNIVERSITY 


Communicated by J. R. Dunning, February 25, 1953 


Linear multipoles, that is, electrical networks with two or more acces- 
sible terminals, have long played a central role in network theory." * In 
this note the notion of a multipole is given a much broader meaning and a 
system of classes of nonlinear multipoles is introduced. An indirect objec- 
tive of this note is to outline a conceptual scheme which, in conjunction 
with the existing theories of oriented graphs* ‘ and with the aid of machine 
computers, could provide an effective tool for the analysis and design of 


nonlinear systems. 
By a multipole, or more specifically, an m + n-pole, we mean a device, 
M, which has m input terminals and output terminals, as shown in 


u, 











FIGURE 1 
Diagrammatic representation 
U; 
o—_e_—_ 





of an m + n-pole 
FIGURE 2 
The sum of M, and M, 


figure |. In the diagrammatic representation, the arrowheads serve to 
distinguish between the input variables, u,,7 = 1,2, ..., m, and the output 
variables, v;, 7 = 1, 2, ,n. The vu, and v, are assumed to be functions 
of time, but for brevity the argument ¢ will be omitted where no ambiguity 
is likely to result. Furthermore, the «,; and v,; are assumed to range over 
specified function spaces, l’, = }u,(t); and 1’, = }v,(f){. In the sequel this 
assumption will often not be stated explicitly. 

An m + n-pole is characterized by a set of 1 input-output relationships 


of the form 


F(t, ta, En ] y (1) 
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where /; represents an operator acting on the ordered m-tuple (1, ts, 


, 
u,). In other words, F; is an operation defined on the product space l, X 
l2 X x U,, with values in |’;. In this sense, the m + n-pole as a 
whole represents an operation on U; &K ls & ... X Uy, to Vi & Ve X 

x V,. A multipole will be called simple if the F; are ordinary functions of 
the u;. The term zero memory is frequently used in the same sense in con- 
nection with nonlinear two-poles. 

The problem of characterizing an m + n-pole in terms of explicit rela- 
tions of the form (1) is of basic importance in the theory of nonlinear 
multipoles. One way of dealing w:th this problem is to consider a system 
of classes of nonlinear multipoles, such as that introduced in this note, and 
approximate a given multipole by one belonging in this system of classes. 
Before entering into a discussion of this approach, it is helpful to introduce 
the following notations and definitions. 

Let M7, and A/, be two m + n-poles characterized by the input-output 
relationships 


M: v, = F(a, tu, 


Ms: v, = Fi (uy, ta, 


respectively. 

M, and Ms are equivalent, 1/; = J/,, if the input-output relationships 
characterizing \/, are identical with the corresponding relationships for 
Mbp, that is, 


F(t, ta, .. +) Um) = F(t, tee, teh j = 1, 2, Nn. (4) 


The sum of M, and Ms, Mo= M, + Mb, is an m + n-pole characterized 
by the input-output relationships 


M: v, = F(a, tee, , Um) + F(a, te, 


The diagrammatic representation for the sum multipole, 7, is shown in 
figure 2. The n + 1|-poles identified by the plus sign represent conven- 
tional adders, an adder being a multipole whose output is equal to the sum 
of its inputs. (Note that the validity of this and other diagrammatic 
representations involving two or more interconnected multipoles depends 
on the assumption that the external characteristics of a multipole are not 
affected by its connection to other multipoles. In engineering terms, this 
is equivalent to assuming that interconnected multipoles do not load one 
another. ) 

The product of ak + n-pole My and an m + k-pole M,; is an m + n-pole 
M, Mos= MoAl,, which consists of 1/7, and A/> connected in series as indi 
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cated in figure 3. From this definition it follows at once that the multipole 
product is associative and left-distributive, that is, 
Ms(MoM,) = (MyM) M, (6) 


and 


(My + Mo)M, = MgM, + MoM,. 


Ann + m-pole, M,~', is called a left-inverse of an m + n-pole 
M,'M=I 


where / denotes a untt multipole, that is, an m + m-pole such that 


v,(t) = u,(t), += I, 2, ,m 


for all u, € L’,. 





FIGURE % 
The product of Mf, and M,. 
FIGURE 4 


Realization of the inverse of A/. 


' such 


Similarly, a right-inverse of A, if it exists, is ann + m-pole, J/, 


that 
MM,” = 1, (10) 


where / is the unit » + n-pole. 
The following two properties of right- and left-inverses are immediate 
consequences of the defining relations (8) and (10). 
(i) If W,~' and My! both exist, then W,>!' = Mp!. 


If \/, ' exists and is unique, then 7, ~' is also the right-inverse of 


(11) 
M. 


When \/,' and Ap ~' both exist, the multipole 


Mo) = M,-' = Mr"! 


is called the inverse of MJ. The determination of the inverse of a given 
multipole is in general a difficult, if not intractable, mathematical problem. 
It is a simple matter, however, to realize approximately the inverse of a 


(11) 
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given physical multipole with the aid of high-gain amplifiers. This can be 
done formally in the manner illustrated in figure 4, which shows the series 
connection of a 2 + 2-pole, \/, and a multipole approximating its inverse, 
M~'. The approximation to the inverse consists of /, two high-gain 
amplifiers A; and Ag, and a pair of subtractors. (A high-gain amplifier is 
essentially a two-pole whose output is equal to G times the input, where G, 
the gain of the amplifier, is a large constant.) It is seen that, as the gain 
approaches infinity, ¢) = 4, — v, and & = uw: — v2 approach zero, the series 
combination of the two multipoles approaches the unit 2 + 2-pole and, con- 
sequently, the multipole labeled 1/~' approaches the inverse of J. 

A multipole is reducible if it can be represented as the sum of two or more 
multipoles having a smaller number of poles. A reducible m + n-pole is 
degenerate if it can be represented as the sum of mn two-poles. The input- 
output relationships for a degenerate m + n-pole are of the form 


m 


v,, = > F,(u,), 
i=l 
where F,, denotes an operator acting on 1,. 
A multipole is /inear if 


F (any! + Buy", ats’ + Buy”, , Ah,’ + Bu,”) = 
aF (uy’, Uy’, ’ tie) + BF (uy", uy", ’ i,” J l, y 


for any real constants @ and § and any time-functions u,’, 4,” 
multipole is nonlinear if it is not linear. 
An obvious but important corollary of (13) is that any linear multipole 
is degenerate. For in consequence of (13), we may write 
Fyi(u, + O + + 0,0 + mw + + (0), ,O+ +0O+ u,,) 
F(m, 0, , 0) + F,(0, tw, ,O) + + F(0, 0, ..., Up) 
m 


> Fy(uj), (14) 
t= 1 


where 
F,(u,) = F,(0, » by, , O). (15) 


Because of the degeneracy of linear multipoles, any linear system may be 
represented as a collection of interconnected linear two-poles arid adders. 
We proceed now to define a certain system of classes of nonlinear m + 1- 
poles for the purpose of systematizing the characterization, realization, and 
optimization® of nonlinear multipoles. The restriction to m + |-poles does 


not involve any loss of generality, since by virtue of (1) any m + n-pole 
may be decomposed into ” m + 1-poles. 


An m + 1-pole is of class 3G, 1, ..., , if its output, o(/), admits a represen- 
tation as an m-fold integral of the form 
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v(t) = / / ao K [u(t — 71), ue(t — 72), ,Uy(t — Tm), 


1}; 72) 2p Tey tlOTy dre... Atm, (16) 


where the u,(t) are the inputs, the 7's are the variables of integration, and 
AK is any real function of its arguments. (A may include delta-functions 
of various orders.) The function A [u(t — 71), ..., Un(l — Tm), Ty 2 
Tm, ¢| will be referred to as the characteristic function of the multipole with 
which it is associated. 

It will be noted that when A does not depend explicitly on f, the m + | 
pole which it characterizes is time-invariant, that is, has the property 


v(t — ty) = Fluy(t — ty), w(t — to), »Un(t — to) | (17) 


for all real constants (), It will also be noted that in the case of any physi- 
cal multipole, that is, one which operates only on the past of its inputs, 
the characteristic function vanishes for all negative values of the r,. 

The class %.; .., is the lowest member of an hierarchy ot 
classes of nonlinear m + 1-poles, of which a general member is denoted by 
Ma give pe Mage, IS defined as the class of all m + 1-poles whose 


input-output relationship is expressible in the form of an (a + 8 + 
+ y)-fold integral 


v(t) = f f ie K [u(t Ti); WE Fug) cy GE Ta) 


Ue(t To), Uo(l — To) gh Taal ys ss 5 
U,(t — Tm)» u,(t as T m2) ’ u,,(t :* Timp) s 
Tity Ti2 » Thay T2ly T225 » T28, + Tmly Tm » Tmys t| 


dry, dtp dtyq dry dt tag ... Atm Atm dtm, (18) 


From this definition it follows at once that MN, »...., © We’ p’,---, x»! 
ifa’ > a, Bp’ > B, ,m’ > yw, and that the class ; of linear m + 1-poles is a 
subclass of every Ny 3. «-- 

To illustrate how an m + 1-pole belonging in this system of classes can 
be characterized in terms of its responses to a family of ‘test’ inputs, we 
consider for simplicity a time-invariant multipole of class %, ;, a 
Let the inputs 4, u%, , U» consist of narrow rectangular pulses of width 
h and amplitudes a, a», , dn, respectively. Furthermore, let these 
pulses occur at the instants f), fo, ,t,. On substituting the expressions 
for the pulses in (16), we obtain the following relation between the charac 
teristic function, A, and the response, v(f), to the set of rectangular pulses: 


K [a,, dh, . » Ba, f Kor t,t tas to, eugene a tm| = v(t). (19) 


h”™ 


(provided v(t) == 0 when one or more of the u;(t) are identically zero). 
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Thus, by varying the amplitudes |a,} and the instants}/,{, one can experi- 
mentally determine the function A and thereby characterize the m + 1- 
pole under consideration via the input-output relationship (16). In other 
words, an m + 1|-pole of class Ni, 1, , is completely characterized by its 


° . ° ° — 
responses to a family of m rectangular pulses whose amplitudes ja,{ and 


times of occurrence }/,} range over all real values. 

Once the characteristic function of an m + 1-pole of class NM, 5... , 
has been determined, the m + 1|-pole may be realized approximately in a 
canonical form such as that illustrated in figure 5. The structure shown 
in figure 5 represents an approximate realization of a physical, time- 


FIGURE 5 


Approximate realization of a three-pole of class WU) .- 


invariant three-pole of class Ji, ; which is characterized by an input-output 
relationship of the form 


v(t) = I I K [u(t — 11), to(t — 2), 71, T2]d 71 dtp. (20) 
0 0 


The structure consists of two tapped delay lines, an adder, and a number 
of simple three-poles (indicated as nodes) whose outputs are expressed by 


Y(t) = fy xla(t — RT), u(t — rAT)], (21) 
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where T denotes the time-delay between two adjacent taps. The output, 
v(t), is given by 


v(t) = Od fe alta(t — kT), u(t — AT)]. 22) 
k=0 =0 

The constant 7’ and the f, , functions in (22) are chosen in such a manner 
that (22) constitutes an acceptable approximation to (20). 

' Cauer, W., “Die Verwirklichung von Wechselstromswiderstanden vorgeschriebener 
Frequenzabhangigkeit,”’ Arch. f. Elektrotechnik, 17, 355 (1926). 

? McMillan, B., “Introduction to Formal Realizability Theory," Be// Syst. Tech. J., 
31, 217 (1952) 

’ Mason, S. J., ‘Signal Flow Graphs,” Proc. 1.R.E., 39, 297 (1951). 

‘Luce, R. D., “Two Decomposition Theorems for a Class of Finite Oriented Graphs,” 
Amer, J. Math., 74, 701 (1952). 

5 The optimization of nonlinear two-pole filters is discussed in a paper by the writer, 
“Optimum Nonlinear Filters,” to be published in J. Appl Phys 


A STUDY OF MUTABILITY IN STREPTOMYCIN-DEPENDENT 
STRAINS OF ESCHERICHIA COLI 
By EpGar L. LABRUM 


DEPARTMENT OF GENETICS, CARNEGIE INSTITUTION OF WASHINGTON, COLD SPRING 
HARBOR, NEW YorRK 


Communicated by M. Demerec, January 380, 19538 


Previous studies of mutability in /schertchia coli have dealt with specific 
mutational changes leading to bacteriophage resistance (Demerec,' Demeree 


and Latarjet,? Witkin®) or to streptomycin resistince (Demerec*°). More 
| i 


recent investigations have been concerned with ‘‘reverse’’ mutations from 
amino acid deficiency to non-deficiency in 45 auxotrophic strains represent- 
ing requirements for | 1 different amino acids ( Demerec,® Demerec and Cahn‘). 
These studies revealed that the various auxotrophic strains differed with 
regard to rates of spontaneous and induced mutation and also with respect to 
patterns of appearance of induced “reversions.’’ For any one auxotrophic 
strain, the pattern of delayed appearance was the same regardless of the 
mutagenic agent used to induce the changes to non-deficiency. Since the 
reversions in at least some of these strains presumably involved different 
genetic loci, these results suggested that loci may differ with re- 
spect to rates of spontaneous and induced mutability and also with respect 
to patterns of appearance of induced mutations. Another finding was that 
although spontaneous reversions occurred readily in the 35 diferent auxo 
trophic strains studied, in 5 of them no induced mutations could be detected 
Such strains were termed ‘mutagen stable.”’ 
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The studies on mutability to be described here were concerned with 
mutational changes from streptomycin dependence to non-dependence in 
different streptomycin-dependent (Sd) strains of /. celt B/r. Crossing 
experiments with the K-12 strain of /. co/1'* have indicated that mutations 
involving streptomycin resistance and dependence occur at a single genetic 
locus. It is possible, therefore, that various streptomycin-dependent strains 
of independent mutational origin represent different alleles of a single locus. 
The purpose of this investigation was to compare the mutation rates and 
the patterns of delayed appearance of induced mutants exhibited by differ- 
ent Sd strains, and also to discover whether or not any of them showed 
mutagen stability. 

Survey of Reversion Frequencies in 70 Sd Strains.— A series of preliminary 
experiments was performed in order to estimate the frequencies of spon- 
taneous and induced reversion in a random sample of 70 strains of independ- 


rABLE 1 
FREQUENCIES OF SPONTANEOUS AND INDUCED MUTATION FROM STREPTOMYCIN 
DEPENDENCE TO NON-DEPENDENCE IN DIFFERENT Sd StRAINS OF EF coli B/R 


STRAIN OF - MUTATION FREQUENCY PER I(? 
E. colt B/r SPONTANEOUS MnClytInpucep 


Sd-21 0 0 
Sd-26 O8 94 
Sd-14 5.7 5, 000 
Sd-24 29 7,380 
Sd-S1 54 2,800 
Sd-SA 250 25,000 
Sd-S5 500 5,000 
Sd-S3 900 & 000 
Sd-S9 4,000 270,000 
Sd-S11 5,000 350,000 
Sd-S7 16,800 206 000 


ent origin. Twelve of these strains were mutants of spontaneous oc 
currence, isolated from independent cultures of /:. coli B or; ten were ultra- 
violet-induced mutants; and the remainder were from the stock culture 
collection of this laboratory. 

Manganous chloride was the mutagen employed throughout the investt- 
gation. The procedure for treating bacteria with the chemical was a stand 
ard one developed at this laboratory,’ and was the same in all expert 
ments. Sd strains to be treated with MnCl were grown in aerated cul- 
tures of nutrient broth containing streptomycin (10 wg. ml), for 24-36 
hours at 37°. The cultures were then centrifuged, washed in 0.3 17 NaCl, 
resuspended in 0.04%, MnChk, and incubated for | hour at 37°. An un- 
treated control sample was resuspended in broth instead of the chemical. 
The viable cell count in both treated and untreated cell suspensions was 
determined by plating suitable dilutions on nutrient agar containing strep 
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‘ 


tomyein (50 wg./ml.) and incubating 48 hours. Suitable dilutions of the 
treated and untreated samples were plated on plain nutrient agar and incu- 
bated for 5 days at 37°. Colonies that grew on the nutrient agar without 


streptomycin represented mutations from dependence to non-dependence. 
The mutants could be either sensitive or resistant to streptomycin. 

The data shown in table | are representative of the results obtained in 
this survey of 70 strains. Most of the frequencies were calculated from 
colony counts of two plates, and thus are only approximations. All strains 
showing spontaneous-reversion frequencies of less than 10 per 10° were 
retested, and their frequencies were calculated from counts of six plates. 

The results showed that none of the 76 strains was mutagen stable; that 
is to say, induced reversions were detected in every strain that produced 
spontaneous reversions. From this it was concluded that if mutagen- 
stable types do exist among Sd strains, they are certainly more rare than 
among the amino acid-deficient strains used by Demerec.® 

Furthermore these data showed that apparently no two of the Sd strains 
were alike with regard to reversion frequency. This may indicate differ- 
ences in mutation rate among the different strains; but there were at least 
two other factors that influenced the observed differences in frequency. 
One was the number of residual cell divisions that the Sd bacteria were 
able to undergo on plain nutrient agar. For example, in several strains 
tested, such as Sd-21 (table 1), no reversions were detected in either the 
treated or the untreated samples when plated on plain nutrient agar. 
Microscopic observation of the plates revealed that the cells had failed to 
divide at all, thus preventing the appearance of any spontaneous or induced 
mutations. This is consistent with the observations of Bertani,'’ which 
indicated that nearly all reversions in the Sd-4 strain that are detected as 
colonies on plain nutrient agar occur during the residual divisions under- 
gone by the cells after plating. Spontaneous reversions that arise pre- 
viously, during growth of the Sd cells in streptomycin broth, either are 
unable to survive or have a selective disadvantage. Background mutants, 
therefore, can be expected to represent an extremely low proportion of the 
reversions observed when Sd cells are plated on plain agar. In contrast 
to these strains that showed no residual growth, many strains such as Sd- 
S11 and Sd-7 (table 1), which underwent six or seven residual divisions on 
medium lacking the antibiotic, exhibited high frequencies of mutation. 

Besides intrinsic differences in mutation rate and differences in extent of 
residual growth, there was variation in the proportion of cells surviving 
MnCl, treatment in different experiments. Survival varied from 10 per 
cent to SO per cent, even though care was taken to perform all experiments 
in the same manner. Past experience with the use of MnCl in this labora- 
tory has indicated some correlation between survival and frequency of 
induced mutations. 
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Mutation Rates in Different Sd Strains.— Experiments were performed 
to determine whether or not 15 Sd strains did differ significantly with respect 
to mutation rates. Each strain was treated with MnCl, as described pre 
viously, and four plates were spread for each dilution of the control and 
treated samples. The number of residual divisions that each strain was 
able to accomplish on nutrient agar without streptomycin was determined 
by microscopic observations. Knowing the number of residual divisions 
and the number of bacteria plated, it was possible to estimate a spontaneous 
mutation rate per cell per generation. Rates of induced mutation were 
calculated on the basis of the number of bacteria plated, before any cell 
division occurred, and the values were corrected for spontaneous, mutations. 
It was probable that all the strains used in these experiments passed through 
a sufficiently large number of residual divisions so that all induced mutations 
were expressed. Each experiment was repeated on three different days. 

The results, presented in table 2, show a high degree of variability among 
the mutation rates determined for each strain at different times. It is 
evident, however, that significant differences do exist among the various 
strains. Sd-U3 and Sd-S4, for example, were characterized by relatively 
low spontaneous mutation rates, whereas strain Sd-S12 exhibited an ex- 
tremely high one. Most of the other strains had rates intermediate be- 
tween these two extremes. 

Pattern of Delayed Appearance of Induced Mutants. Demerec! reported 
that cells of /. colt treated with either x-rays or ultra-violet radiation must 
pass through 10 to 12 divisions before all induced phage-resistant mutants 
appear. This ‘‘delayed effect” or ‘“‘phenomic delay"’ was also observed by 
Davis'' in the isolation of induced nutritionally deficient mutants of /. 
coli by the penicillin method. The phenomenon will be referred to in this 
paper as the “delayed appearance of induced mutants.’’ Current views 


of the possible mechanisms responsible for delayed appearance have been 
discussed by Witkin'* and Demerec.* 
By plotting frequency of induced mutants against number of cell divi- 


sions elapsed since treatment, one obtains a curve representing the rate at 
which induced mutants appear, or the “‘pattern of delayed appearance.”’ 
In cases where there is considerable variation among the frequencies obtained 
in different experiments, the percentage of the total yield of induced muta- 
tions manifested at each cell division can be plotted instead of the frequency. 

From the 70 Sd strains used in the summary described previously, a dozen 
strains were chosen as being possibly suitable for these experiments. The 
strains selected were ones that pass through a minimal number of residual 
divisions on plain nutrient agar. 

It was found that if the nutrient agar was “enriched” with low concen- 
trations of streptomycin the number of residual divisions of some Sd strains 
could be increased. Thus, by plating the MnCl.-treated and control sam 
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TABLE 2 


RATES OF SPONTANEOUS AND MnCl-INDUCED MUTATION FROM STREPTOMYCIN 
DEPENDENCE TO NON-DEPENDENCE IN DIFFERENT Sd Srratns OF E. coli 


STRAIN OF NO. OF RESIDUAL MUTATION RATES 
KE. cola B/t DIVISIONS SPONTANEOUS (XX 10%) INDUCED (& 10%) 


Sd-U3 6 2.4 10 


‘ 
» 


- 


Sd-S10 


Sd-S11 
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ples on a series of agar plates containing increasing concentrations of strep- 
tomycin, one could determine the reversion frequency after each successive 


division. The actual range of streptomycin concentrations that was 
effective in this respect differed for the various Sd strains, but in no case 
did it exceed | wg per ml. 


Fig. 1. Rate of appearance of 
induced reversions in the strepto- 
mycin-dependent strain B/r 
Sd-34 of E. coli. The data plot 
ted represent the results of three 


ro) 
° 


ra 
° 


independent experiments 


ae 
° 


Fig. 2. Rate of appearance of 
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induced reversions in the strepto 
mycin-dependent strain B/r 

Sd-138 of FE. coli, The data plot 
ted represent the results of three 
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Fig. 3. Rate of appearance of 
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mycin-dependent strain B/r 
Sd-S1 of E. coli. The data plot 


my hecentecsttel — ' ted represent the results of three 
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An unexpectedly large increase in number of colonies occurred when the 
Sd cells were plated on medium containing streptomycin in concentrations 
of 0.4 to | wg. per ml., and it was noted that these colonies varied greatly 
in size. When some of the smaller colonies were picked and streaked on 
both plain and streptomycin agar, it was seen that they were not true rever- 
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sions but apparently “‘semidependent”’ mutants, which could grow to some 
extent on very low concentrations of streptomycin. Since the mutation 
frequencies to be determined in these experiments were to represent only 
complete reversions, it was necessary to correct all colony counts. This 
was done by picking a random sample of at least 50 colonies from each set 
of plates that contained a diiferent concentration of streptomycin, and 
streaking them out on both plain and streptomycin agar. From the results 
obtained it was possible to estimate what percentage of colonies represented 
true reversions for each concentration of streptomycin used in both control 
and treated series. Most of the Sd strains originally selected for these 
experiments gave rise to such a large proportion of the ‘‘semidependent’’ 
type of mutants that they could not be used. 

The three strains best suited for the experimental procedures were Sd-Sl, 
Sd-34, and Sd-138. In figures 1, 2, and 3 percentage of induced mutants 
expressed is plotted against number of divisions for each of these three strains. 
The curves thus obtained represent the patterns of delayed appearance. 
By comparing the curves it can be seen that in each case 100 per cent ex- 
pression was reached at about the fifth or sixth division. Fifty per cent 
expression occurred at about 3 to 4'/, divisions. It was therefore concluded 
that the three Sd strains exhibited similar patterns of delayed appearance 
of induced mutants. 

Discussion. In a previous analysis of 163 Sd mutants, Demerec! found 
that practically no two strains were alike when compared on the basis of 
six different criteria, including spontaneous mutation rate and proportions 
of streptomycin-sensitive and streptomycin-resistant reversions produced. 
The data presented in table 2 extend these observations to include induced 
mutation rates. Considerable variation was encountered among separate 
experiments done with the same strain on different days; however, the 
rates for some strains were consistently low, whereas others showed consist- 
ently high or intermediate rates. Although no absolute correlation was 
found between spontaneous and induced rates, it was evident that there 
was, in general, a tendency for strains having low spontaneous rates to 
show relatively low induced rates, and those with high spontaneous rates 
to exhibit high induced rates. 

Determination of the pattern of delayed appearance of induced mutants 
for three Sd strains (Figs. 1, 2, and 3) indicated that in each strain 100 
per cent expression of induced mutations was reached at about the fifth 
or sixth division. In other respects, however, these three strains were not 
similar. The reversions produced by Sd-Sl were all streptomycin sensi- 
tive; those of Sd-13S were all resistant; and Sd-34 produced approximately 


63 per cent resistant and 37 per cent sensitive reversions. They also dif- 


fered with regard to the number of residual divisions they were able to 
undergo on plain agar. 
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The pattern of delayed appearance exhibited by the Sd strains appedrs 
to be an unusual one since 100 per cent of the total yield of induced muta- 
tions was obtained earlier than in most other mutation systems previously 
studied. With mutations involving resistance to coliphage T1, 100 per 
cent expression was not reached until the tenth or twelfth division.' Most 
of the auxotrophic mutant strains that have been analyzed®? also required 
more than six divisions to achieve 100 per cent expression, although in two 
amino acid-deficient strains all mutants had appeared by the end of the 
first division. 

If one views the different Sd strains of independent origin as representing 
different alleles at a single locus, then it would seem likely that each different 
allele has its own characteristic rates of spontaneous and induced mutation. 
The fact that in three Sd strains all induced mutants had appeared by the 
fifth or sixth division, however, suggests that the pattern of delayed ap- 
pearance is possibly determined by the genetic locus involved rather than 


by the particular allele. 
Summary.— Frequencies of spontaneous and induced reversion from strep- 


tomycin dependence to non-dependence were estimated for 70 different 
Sd mutants of /. coli B/r. The results indicated that apparently no two 
strains were alike in this respect, and that no mutagen-stable types existed 
among the 70 strains. 

Estimation of mutation rates from streptomycin dependence to non- 
dependence in 15 strains indicated wide differences among some strains. 

The pattern of delayed appearance of induced mutants was determined 
for three different Sd strains. In each case, 100 per cent expression of 
induced mutants occurred by about the fifth or sixth division. 

These results indicated that Sd strains of independent origin may be 
characterized by different rates of mutation to non-dependence, but that 
the pattern of delayed appearance of induced mutants 1s possibly similar 
for the different strains. 


Acknowledgment.— The author wishes to express his gratitude to Dr. M. 
Demerec for his helpful advice and discussions during the course of this 
investigation. 
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INDUCED POLLEN LETHALS FROM SEEDS OF DATURA 
STRAMONIUM TREATED WITH THERMAL NEUTRONS* 


By J. L. SPENCER, W. R. SINGLETON, AND A. F. BLAKESLEE 


UNIVERSITY OF MASSACHUSETTS, BROOKHAVEN NATIONAL LABORATORY, AND SMITH 
COLLEGE GENETICS EXPERIMENT STATION 


Communicated January 19, 1953 


The cytogenetic effect of thermal neutrons was shown by Conger and 
Giles! to be about eleven times that of x-rays per unit dose. The difference 
in effectiveness is believed to be due largely to the internal origin of the pro- 
tons and alpha rays which result from capture reactions. 

A preliminary study of the effects of thermal neutrons on plants of 
Datura stramonium grown from treated seeds was undertaken for the 
purpose of discovering the induced cytogenetic aberrations, especially 
pollen lethals. Cartledge and Blakeslee*® in a study of genetic changes in 
aged seeds have shown that pollen abortion in Datura affords a convenient 
and delicate index to the mutation rate. They distinguished two types of 
pollen abortion. In one, the chromosomal type, the pollen grains abort at 
an early stage resulting in small, shriveled, and empty grains. The other is 
marked by subnormal to very small size and usually by differences which 
were interpreted as degenerative changes in the character of their contents. 
Of ten plants listed as belonging to the gene type, Dr. Satina found all 
were free from obvious chromosomal abnormalities thus confirming their 
genic nature. Of 11 classified as chromosomal types all except three showed 
chromosomal abnormalities. These discrepancies may be due to securing 
chromosomal figures from different sectors from those from which the pol- 
len was recorded but more likely was due to sone of the lethal genes causing 
abortion of pollen at an early stage. In the majority of cases, however, it 
appears possible to distinguish pollen abortion types due to genes from those 
due to chromosomal abnormalities. Since further evidence in support of 
this distinction is given in the succeeding paper* we shall continue to use the 
terms chromosomal and gene types for the two rather distinct types with 
the admission that the distinction may not always be ‘ustified. In the 


present study the chromosomal type of aborted pollen occurs in rather 
definite proportions including about 25 per cent, 50 per cent, 75 per cent, 
or rarely 100 per cent of the grains. The gene type causes the abortion of 
either 25 per cent or 50 per cent of the total grains. Mutations of both 
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types arise naturally and are induced readily by exposure to ionizing radia- 
tions, such as x-rays. The induced mutations generally occur in limited 
sectors of the plant, although they often involve the entire plant. Single 
plants and, indeed, even single flowers, may show several types of pollen 
abortion due to as many mutations. 

From the foregoing it is apparent that the two types of pollen abortion 
may be used as a rough index of the rate of both chromosomal and gene 
mutations. Such an index avoids the necessity of growing the large num- 
ber of plants that would be required in the second generation to disclose 
recessive ‘visible’ mutations. Previous work has shown that pollen abor- 
tion types are several times as frequent as recessive morphological types 
disclosed in the second generation. 

For these experiments we employed the standard line | of Datura stra- 
monium which has been selfed for over 40 generations and brought through 
a haploid generation. The seeds from three capsules that had developed 
from flowers having less than one per cent bad pollen were each divided into 


TABLE 1 


POLLEN ABORTION IN DATURA FROM THERMAL NEUTRON TREATED SEEDS 


LENGTH NUMBER NUMBER PLANTS EXAMINED 
EXPOSURE OF OF PLANTS WITH WITH 
HR SEEDS SEEDLINGS MATURING NORMAL POLLEN ABNORMAL POLLEN 


2 876 742(85%) 647 (74%) 128 (24%) 395 (76%) 
$ 890 683 (76%) 321 (36%) 30 (14%) 191 (86%) 


6 911 205 (22%) 32 (3.5%) 2(11%) 17 (89%) 
Control 150 139 (93°) 114 (760%) 182 (100°, ) O(0&%) 
Additional controls (Line 1) 82 


four lots. One lot of seed from each capsule constituted the control. 
The other three groups were exposed for two, four, and six hours, respec- 
tively, in a thermal column to a flux of approximately 4.6 & 10° Ne 
em.?/sec. (The cadmium ratio was 5000: 1 and the gamma contamination 
was about 50 r/hr.) 

The pollen was examined from a single flower from each fork of the plants 
tested. Datura plants usually have only two forks although three or 
rarely four may be present. If the flower showed abnormal pollen, pollen 
of adjacent flowers was examined in order to delimit the sector exhibiting 
the pollen abnormality. 

Approximately 93 per cent of the control seeds (untreated) produced 
seedlings, while only 85 per cent, 76 per cent, and 22 per cent of the seeds 
exposed for two, four, and six hours, respectively, resulted in seedlings of 
sulicient development to show true leaves (table 1). The S76 seeds sub- 
jected to a dosage of two hours gave rise to 647, or 74 per cent, mature 
plants. This percentage was slightly less than that yielded by the con- 
trols (76 per cent). The four-hour exposure of S90 seeds ultimately re- 
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sulted in $21, or 36 per cent, mature plants, but only 32, or 3.5 per cent, 
of the 911 seeds treated for six hours gave rise to mature Datura plants. 

The mutation rate, based on pollen abortion, was found to be 76 per 
cent, S6 per cent, and 89 per cent for plants from seeds treated for two, four, 
and six hours, respectively. No control plants were found that yielded over 
four per cent bad pollen; however, previous accounts indicate that a muta- 
tion rate of about 0.5 per cent may be expected. 

For the determination of the types of pollen abortion and sectoring in- 
duced by the dosage of thermal neutrons reported here, the pollen from 3517 
flowers was examined (table 2). An analysis of the 2034 flowers from plants 
that arose from seeds exposed for two hours to thermal neutrons showed 
685 which’could be considered normal (less than 10 per cent bad pollen). 
The four-hour treatment resulted in 142 normal flowers out of the 897 
examined, while the six-hour dosage yielded only seven flowers with less 
than 10 per cent aborted pollen out of the 84 analyzed. In other words, 


TABLE 2 
TYPES OF POLLEN ABORTION FOUND IN DATURA FROM THERMAL NEUTRON TREATED 
SEEDS 
POLLEN CONDITION 
CHROMOSOMAL 
AND 
LENGTH OF NUMBER OF CHROMOSOMATI ABNORMAL GENE TYPES 


EXPOSURE FLOWERS TYPE GENE TYPE (SHRIVELED TOTAL 
Hk ANALYZED NORMAL (SHKIVELED) (SMALL) AND SMALL) ABNOPMAL 


2 2034 685 (340% ) 1064(52°,) 158(89)) 1382(6%%) 1349 (66%) 

1 807 142 (16>) 603 (67%)  57(7%) 95 (11%) 755 (84%) 

6 84 7 (8%) 57 (68°) 3(4%) 17 (20%) 77 (92%) 
Control 502 502 (100%) 0(0%) 


only 34 per cent, 16 per cent, and 8 per cent of the flowers from plants 
that arose from seeds treated for two, four and six hours, respectively, 
were comparatively free from pollen lethals. All of the flowers from con- 
trol plants were found to be normal with respect to aborted pollen. 

The presence of pollen lethals resulted in flowers that contained aborted 
grains of either the chromosomal type, or the gene type, or a combination 
of both types. The latter condition, both shriveled grains and small grains 
in the same flower, indicated the presence of at least two mutations, one 
of each type. The two-hour dosage resulted in chromosomal types in 
52 per cent of the flowers, gene types in 8 per cent, and both types in 6 per 
cent, or a total of 66 per cent (1849) abnormal flowers. The four-hour 
exposure yielded chromosomal types in 67 per cent of the flowers, gene 
types in 7 per cent, and II per cent with both conditions. A total of S4 
per cent (755) of the flowers examined from plants from seed exposed four 
hours was abnormal. The six-hour dosage resulted in 92 per cent (77) of 
the flowers being abnormal: 68 per cent demonstrated chromosomal types 
of abortion, + per cent showed gene types, and 20 per cent both types. Of 
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course the number of flowers in each category would depend upon the num 
ber examined from each sector. However, the percentages stated give 
approximation of the severity of treatment and of the resultant induction 
of sectors. 

The evident lethal effects on the seeds and seedlings of Datura caused by 
the severity of treatment at the longer exposures of four and six hours made 
it difficult to reach definite conclusions regarding the genetic effects of 
different dosages. However, if flowers with both shriveled and small grains 
are considered to have two types of mutations, then there have been in 
creases in the percentages of both types with further increments of time. 


If the double types, chromosomal and gene, are taken into consideration 
there is no significant change in the ratio of chromosomal to gene types with 


increased dosage. 
Large pollen grains were found in some of the flowers of exposed plants, 
but not in the flowers on the control plants. Grains of diploid (2N) size 


TABLE 3 
PoLLEN ABORTION SECTORS IN DATURA FROM THERMAL NEUTRON TREATED SEEDS 
NUMPER OF PLANTS HAVING DISTINCT 
NUMBER OF PLANTS HAVING SECTORS BETWEEN 
LENGTH OF POLLEN LETHALS NORMAL NORMAL CHROMOSMAI 


EXPOSURE WITH WITHOUT CHROMOSOMAL GENE GENE 
HR SECTORS SECTORS" ryPer TYPE TYPE OTHERS? 


2 288 107 (27%) 96 19 32 141 
} 132 59 (381%) i 1 3 10 
6 9g 8 (470%) : l 1 1 
All flowers of plant show a single pollen abortion type. 

” Sectors recorded, but individual flowers show both pollen abortion types 


were recorded in the following percentages: two-hour treatment, 7 per 
cent of the plants; four-hour treatment, 10 per cent of the plants; and six 
hour treatment, 21 per cent of the plants. Pollen of larger than diploid 
size was rarely noted. 

The number of plants that developed from treated seeds and that dis 
plaved but a single type of pollen abortion (i.e., either chromosomal or 
gene) in all their flowers increased from 27 per cent in the two-hour group 
to 47 per cent in the six-hour group (table 43). This would seem to indicate 
that the severity of the longer exposures induced mutations affecting pollen 
abortion which involved the entire plant. The pollen has its origin in the 
second layer of the embryo. Perhaps adjacent cells of this layer were 
killed by heavy dosage and only one cell which was lightly hit survived to 
give rise to pollen grains. With very few exceptions plants which demon- 
strated a single pollen abortion type throughout possessed a chromosomal 
type of mutation. A significant positive correlation appears to exist be 
tween pollen abortion and morphological aberrations; however, the severity 
of treatment at both the four- and six-hour levels tended to obscure this 
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pattern by making it difficult to delimit sectors clearly. Although the 
data at hand precluded the possibility of making rigorous comparisons be- 
tween plants arising from seed subjected to thermal neutrons and plants 
from x-irradiated seed, no developmental or morphological disturbances 
or lethality patterns in ‘Jatura unique to either kind of treatment were 
detected. Caldecott, Frolik, and Morris‘ found that individual barley 
seed subjected to identical doses of x-rays were variously affected, as deter- 
muned by seedling heights, but injury was much more uniform with identical 
doses of thermal neutrons. These investigators also reported that ‘‘dif- 
ferential” delayed killing was characteristic in barley exposed to 10,000 to 
20,000 r units of x-radiation, but that similar lethality patterns did not 
occur in thermal neutron-treated material. 

Further study of thermal neutron induced pollen lethals in Datura is 
in progress. 

* Contribution from the Department of Botany, Smith College, New Series No. 48 
This work has been aided by a grant from the Atomic Energy Commissiou 

' Conger, A. D., and Giles, N. H., Genetics, 35, 8397-419 (1950) 

2 Cartledge, J. L., and Blakeslee, A. F., Proc. Nari. Acab. Scr., 20, 103-110 (1934) 

* Yost, H. T., Jr., Singleton, W. R., and Blakeslee, A. F., Proc. Nati. AcAb. Sct., 39, 
292-297 (1953) 
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THE EFFECT OF THERMAL NEUTRON RADIATION 
ON THE CHROMOSOMES OF DATURA* 
By Henry T. Yost, JR., W. R. SINGLETON, AND A. F. BLAKESLEE 
BIoLoGIcAL LABORATORY, AMHERST COLLEGE, AMHERST, MASSACHUSETTS; BROOK 


HAVEN NATIONAL LABORATORY, Upton, NEW YorRK; AND GENETICS EXPERIMEN] 
STATION, SMITH COLLEGE, NORTHAMPTON, MASSACHUSETTS 


Communicated January 19, 1958 


The work of Conger and Giles' has demonstrated the great effective- 
ness of thermal neutron radiation in producing chromosomal aberrations. 
This finding has led to the investigation of the effects of thermal neutrons on 
plant growth’ and on polyploid material.’ In the light of this work, it 
seems advisable to investigate the genetic effects of thermal neutron radia- 
tion and, if possible, to correlate them with chromosomal aberration. The 
general problem is the relationship between chromosomal aberration and 
the expression of genetic changes. This problem may be investigated by 
discerning whether increased chromosome breakage yields proportionate or 


disproportionate genetic change. The studies presented herewith are pre 


liminary to the solution of this problem. 
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One of the great difficulties in studies correlating gene and chromosomal 


effects is the time which it takes for recessive mutants to appear. Indeed. 
Tradescantia, on which the original work was done,' has no genetic markers; 
consequently there are no data available on the genetic effects of thermal 
neutrons. In Datura there is available an ideal test organism. A meas 
ure of the over-all mutation rate may be gained from the pollen abortion of 
the mature plants developing from irradiated seed. The basis for this 
estimate has been described elsewhere.‘ °’ In addition, there are two 
distinct types of aborted pollen: small, empty, shrivelled grains and very 
small grains with more or less degenerate contents. (For brevity these two 
types will be referred to as empty-shrivelled and small.) Cartledge and 
Blakeslee® have indicated that there is a correlation between these two types 
of abortion and the gene and chromosomal types of mutation. They have 
suggested that the small grains are associated with point gene mutation, 
and the empty-shrivelled type with chromosomal mutation. This distine- 
tion was confirmed cytologically by Dr. Satina from study of a limited num- 
ber of cases of each type. It was deemed important to test this hypothesis 
with the neutron radiation, since, if it were true, it would provide a rapid 
survey method for a much larger number of plants than could be investi 
gated cytologically. The work reported here bears specifically on this 
problem. 

Experimental Procedures and Results. Seeds of Datura stramonium, 
Line 1, were exposed to radiation in the thermal column of the nuclear re 
actor at the Brookhaven National Laboratory for a period of two hours. 
The flux was approximately 4.6 & 10°N,,/em.*;/sec. The cadmium ratio 
was about 5000: 1 and the gamma contamination was about 50 r/hr. The 
control seeds and the radiated seeds came from the same parents. All 
seeds were germinated in the greenhouse and subsequently transferred to 
the field. The mature plants were studied for pollen abnormalities and 
sector of good and bad pollen delimited for each plant. Buds were then 
obtained from both the good and bad sectors, and aceto-carmine smear 
preparations of the anthers were made. The metaphase and anaphase of 
the first meiotic division were scored for abnormalities. 

Segmental interchanges are the most easily detected chromosomal 
aberrations in Datura. In this case, there is the formation of a circle or 
chain of 4, 6, 8, ete., depending upon the number of chromosomes involved. 
Translocations involving the satellites are also scored readily. In general 
they give rise to a configuration similar to a bow-tie, with the translocated 
parts joined at the center. Large deficiencies can be scored as well as de- 
ficiencies or duplications of whole chromosomes. No attempt was made to 
determine the frequency of inversions, and it is probable that small defi 
ciencies would be missed unless the fragment remained in the cell. It is 
evident from these studies and from others® 7 that segmental interchanges 
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give a reliable measure of the rate of breakage. Wherever there was any 
doubt about the nature of a configuration, several determinations of the 
same sector were made. Ina large majority of the cases the configurations 
were confirmed by Dr. Satina, to whom we are greatly indebted. 

The results of the comparison of pollen abortion types are summarized in 
table 1. It is immediately evident that there is a correlation between the 


rABLE 1} 


CYTOLOGICAL CHECK ON POLLEN ABORTION TYPES FROM NEUTRON-IRRADIATED SEEDS 
(2 Hours4.6  10* Nyu/Cm.2/Sec.) 
POLLEN ABC ETION, POLLEN ABORTION 
CHROMOSOMAL TYPES GENE TYPE 
POLLEN POLLEN 

PLANT NO ‘ CHROMOSOMES PLANT NO % CHROMOSOME 
SIOL75-18 i 4) + 11 bvts.:* 12 bvts 510175-2 50 12 bvts. 
510175-338 § 12 bvts. 5IOL75-18 5O 12 bvts. 
5LOL75-40 2! ‘ 510175-63 50 4 
510175-52 df 12 bvts 510175-93 dO 12 bvts. 
DIOL 5-55 5 4) + 11 bvts. 4+ frag * 

12 bvts DIOL75-111 25 12 bvts. 
5IOL75-56 4 510176-32 50 12 bvts. 
510175-68 f 4 510176-97 25 12 bvts 
SIOLT5-92 5 Chain of 6 51O176-118 50 12 bvts. 
5IOL75-114 2 i 510175-20 25 12 bvts 
SIOL76-1 Sat. trans 510176-69 25 4 
5IO1L76-14 2! ‘ 510176-64 a 12 bvts 
5LOL76-25 1 510176-70 50 12 bvts 
SLOL76-380 : ‘ 

510176-82 5 ‘ 
510176-38 a 4 
510176-52 Sat. trans. 
510176-63 2! ‘ 
510176-67 25 12 bvts. 
DIOL76-83 Sat. trans. 
SLOLT6-86 12 bvts. 
510176-97 
5SIOL76-5 
510176-26 
510176-189 Jat. trans 
JSIOL76-118 

Retested. 


types of pollen abortion and the types of mutation. Chi-square tests indi- 
cate that there is a highly significant difference between the two types of 
abortion with regard to the proportion of chromosomal aberration present 
in each case. The P value determined for this comparison is 0.0005. 
There can be little doubt that the two types of pollen abortion show dif- 
ferent chromosomal patterns. Table 2 shows the results of studies made on 
good sectors and plants; that is, on plants showing less than 10°; pollen 
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abortion of any type. It is evident that there is little indication of chro- 
mosomal aberration in this group. It is interesting to note that one of the 
normal sectors showing a segmental interchange is from the same plant 
(510175-63) which showed a segmental interchange in the small pollen 
group of table 1. This is an indication that there is reason to doubt that 
the small pollen is associated in any way with the segmental interchange. 
Taking both tables together, it 1s evident that there are some cases in which 
segmental interchanges do not give rise to pollen abnormalities. However, 
since in all eases in which empty-shrivelled sectors showed segmental inter- 
changes the good sectors on that same plant showed no segmental inter- 
changes, there is good evidence that the two effects are related. There is 
good reason to believe, on the other hand, that in cases of small pollen show- 
ing segmental interchanges there is no relation between the two other than 
one of chance association. 

Plants number 510175-18 and 510175-55 are worthy of separate con 
sideration. In each of these an extra pair of chromosomes was apparent on 
the first determination. Subsequent determinations, however, failed to 


TABLE 2 


CyTroLoGICAL CHECK ON PLANTS SHOWING NORMAL POLLEN 


POLLEN ABORTION, °, NUMBER OF PLANTS CHROMOSOMES 
0-1 30 12 bvts. 
0-1 Ss 4 
2-5 17 12 bvts. 
5-10 6 12 bvts. 





' SI0176-69; SLOLT6-49, 





show this extra pair. In both cases, the extent of the damage to the chro- 
mosomes was apparently limited to one bud, or at least to a very small 
sector. 

Discussion. The data just cited indicate that it is permissible to group 
the aborted pollen into two distinct types. These two types are different 
morphologically and genetically. Grains which are empty and shrivelled 
are due to an early abortion and seem to be positively correlated with chro 
mosomal aberrations. There are several cases in which no chromosomal 
aberration was observed; nevertheless, the only types of aberration observed 
were of the translocation type, and it is evident that inversions and small 
deficiencies could go undetected. Therefore, it seems quite safe to term the 
empty-shrivelled pollen a chromosomal type. On the other hand, the small 
pollen group showed very few chromosomal abnormalities. It is true that 
small deficiencies would be easily missed in this case, and there is no basis 
for a distinction being made on the old question of point gene mutation ver 
sus the infinitesimal deficiency. Still, the lack of segmental interchanges 
makes a broad distinction permissible, and the small pollen may be classed 
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as gene type. Such a distinction makes it possible to set,up a rapid survey 
method for the whole field of plants, classing them for number of mutants 
and mutant type simultaneously from the morphology of the pollen. This 
should be of great aid in the investigation of future radiation problems in 
Datura. The question of deficiencies in general deserves further attention 
in the near future. There is little evidence into which group small deficien- 
cies would fall. 

Considering the general problem of the relation of the gene and chromo- 
some mutant types developing from neutron radiation in the light of these 
data, it is evident that the gene type are few in number compared to the 
chromosomal type. In fact almost 70°, of the plants showed chromosomal- 
type pollen. At low doses it would be expected that the gene type would 
far outnumber the chromosomal type of mutation, especially since the aber- 
rations scored were of the ‘“‘two-hit’ type. At higher doses, however, the 
reverse relationship might hold. That is, more of the treated cells would 
tend to have chromosomal aberrations, and the gene type would be ob- 
scured. ‘There is some indication of this given in table 1. It will be noted 
that the great proportion of the types showing segmental interchanges 
have 50% bad pollen. In previous studies, it has been demonstrated that 
segmental interchanges in Datura may give rise to 25% bad pollen, but 
generally have good pollen.’ There is a striking difference between the 
data collected with thermal neutron radiation and those from x-rays, 
heat, aged seed, etc. This difference may be explained on the basis that 
the dose of neutrons was great enough to give a large number of cells with 
both chromosomal aberrations and gene mutations, and these cells give 
rise to flowers with 50°, bad pollen as a result of the two types of aberration. 
On the other hand, the possibility still remains that the segmental inter- 
changes induced by thermal neutron radiation are different in extent from 
those induced by other agents. It is quite likely that the aberrations scored 
consist of complexes of breaks rather than simple segmental interchanges 
resulting from two breaks. These data clearly support the findings of 
Conger and Giles! that the thermal neutron radiation is extremely effective 
in producing chromosome breaks. The earlier work of Frolik and Morris* 
indicated that there was little difference in the effectiveness of x-rays and 
neutrons in the induction of sectors in corn. Our data seem to indicate 
that the neutrons are far more effective. The solution of the general 


problem, however, must await further studies at lower doses of neutron 


radiation. 

Final consideration must be given to the two cases which showed extra 
chromosomes. These plants were numbers 510175-1S and 510175-55 
(table 1). The original preparations derived from these plants showed cells 
in which extremely complex aberrations had been induced. In plant 55 
particularly there was evidence of very extreme damage and only the mini- 
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mum found is reported here. The cells were aberrant to the extent that 
it was hard to find metaphase figures showing clear configurations, although 
the slide was good. 

An extra chromosome, at least in the second layer,‘ is generally marked 
by some easily detectable morphological abnormality. One would expect 
that the presence of an extra bivalent should have produeed very clear 
morphological differences. Nevertheless, when the plant was examined 
there was no morphological abnormality seen. For this reason more slides 
were made of this particular sector. Upon examination there was no 
evidence found of chromosomal abnormalities. There are two possibilities 
to explain this occurrence. One is that a single cell was affected and sur- 
vived to give rise to the only bud examined in each plant at a much later 
date. If this be the case, it is reasonable to suppose that the cell never 
became involved with more than a limited portion of the plant, since there 
was no morphological abnormality detected. This is scarcely a plausible 
explanation, since innumerable cell divisions must have taken place after 
the induction of the aberrations. On the other hand, it is possible that 
these two buds represent cases of delayed effects of the neutron radiation. 
It is possible that activated molecules were induced by the neutron capture 
and that later decay might result in the release of energy sufficient to cause 
chromosomal abnormalities. The major difficulty with such an hypothesis 
is that one must assume a relatively stable condition, since it is of very 
long duration, and since the ‘active state’’ must pass through several cell 
generations. On the whole, however, the latter of the two possibilities 
seems the more plausible. Certainly this suggestion is worth further 
study in the near future. 

Summary.— Data have been presented which indicate that pollen abor 
tion types in Datura may be classed in two distinct groups with respect to 
the type of mutation involved. Those showing empty-shrivelled grains 
result from chromosomal aberration, and those showing small grains result 
from gene mutation. With thermal neutron radiation there is a high pre- 


ponderance of the chromosomal type of aborted pollen. The data suggest 


the possibility of extremely delayed effects induced in the chromosomes 
of Datura after thermal neutron radiation. 


* Contributions from the Department of Botany, Smith College, New Series, No. 49 
Phis work was supported by contract with the U.S. Atomic Energy Commission 
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ON SOME VARIATIONAL PROBLEMS OCCURRING IN) THE 
THEORY OF DYNAMIC PROGRAMMING 


By RICHARD BELLMAN, IRVING GLICKSBERG, AND OLIVER GROSS 
THe RAND Corporation, SANTA MONICA, CALIFORNIA 
Communicated by J. von Neumann, February 5, 1953 


1. Introduction. The purpose of this note is to present some results 
of our investigation of a class of interesting and important variational prob- 
lems arising in connection with problems involving the control of a sys- 
tem over a time interval. We may, on one hand, desire to maintain a 
dynamic system of mechanical, economic, or other origin, in or near a speci- 
fied state at minimum cost, where the cost is compounded of two parts, 
the first part measured in terms of the deviation of the system from the de- 
sired state, and the second part measured by the cost of the resources used 
for control. Alternatively, we may, given a fixed quantity of resources, 
wish to maximize or minimize the output of a system. 

The mathematical difficulties encountered in treating problems of this 
type depend upon the mathematical model used to represent the system, 
the functionals employed to measure the cost of deviation and the cost of 
control, and the constraints imposed upon the permissible type of control. 

We shall present here only some results we have obtained for systems 
whose state at any time / is described by a column vector x(t) satisfying a 
linear vector-matrix equation 


dx 


= Ax + f(t), x(0) = ¢, 
dt 


where A is constant. 

An equation of this type enters in connection with problems of continuous 
control. If we take account of time-lags, in place of (1.1) we obtain dif 
ferential-difference equations, the simplest type of functional equation 
allowing for hereditary effects. If in place of continuous control, inter- 


mittent control is used, (1.1) will be replaced by a difference equation. 


With reference to (1.1), let v(¢) be the desired state of system at time ¢, 
and let F(x y) be the functional measuring the total cost of deviation 
over an interval |[0, 7'], G(f) be the functional measuring the cost of control 
as exercised through the foreing term /(/), and //(c) be a function measur- 
ing the cost of control as exercised through starting the system in state ¢. 

We may then minimize the total cost J/(f, c) = F(x — y) + Gif) + (eo), 
or minimize f(x — y) subject to G(f) <a), (ce) < dy. There is an equiva- 
lence between these two types of problems, the second of which leads di- 
rectly, via Lagrange multiplier techniques, to Sturm-Liouville problems of 
non-classical type. If time-lags are permitted, we are led to consider 
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Sturm-Liouville problems for diiferential-diference equations. These 
topics will be discussed in subsequent publications. 

We shall present here some results we have obtained for the cases where 
F and G are linear or quadratic functionals and the constraints upon f are 
of simple type such as 0 < f; < m,. For simplicity we shall assume that 
cis fixed. 

2. Linear and Quadratic Functionals..-The methods used to obtain 
the following results are interesting mixtures of abstract and analytic tech- 
niques reinforcing each other at critical points. What constitutes the 
difficulty of the problem, of course, is that free variations are only occa- 


sionally permissible. 

THEOREM |. Let the scalar quantity wu satisfy the first order equation, 
du/dt = —u+f(t),u(0) = 1. Then the minimum of f°’ (1 — wu)? dt over 
all f subject to0 <f <a,>1,fu' fdt < a, < T, is furnished by 


0 0<t<as; 
f(t) = a, a3 St< ads 


az ( ds aa t < ep 


’ 


where ds, (4, d; are known constants determined by complicated transcendental 
equations. 
The function, f, which under the same conditions yields the minimum of 


So" (du/dt)? is 


f(t) = ag(t — ay)? + ar(t — as), 0O<t<as 
- 0 tes. b= J. (2.2) 


where again dg, d;, and as are constants determined by T. 

Tueorem 2. /f u satisfies the above equation, the minimum of J,’ f? dt 
subject to the constraint —ay < 1 — u < dy,0 < ay <1, forO <t < T, ts 
furnished by, 


for O<T < log b, (b = 1 as), f=9, 
for = logb<T < loge, (c = b/(1 — V1 — 6%)), 
f=2(b-e ee! 

for log ¢ < 7, f = - V1 b*)e', O<t < loge, 

b ; lore SH < T 


3. Linear Functionals..- The problems arising in connection with maxi- 
mizing or minimizing  % (x(t), dG(t)) subject to constraints of the form 
: + (f, ¢) dt = a, 0 < f; < m,, where (y, 2) denotes the inner product, 
may readily be converted into Neyman-Pearson type and do not possess 
any new features of interest. 

A more interesting problem is one of the following sort: 
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THEOREM 3. Consider the equation 


dx; ~ ; : 
=: > ¥ ayy, t= 1,2,...,N, x0) = c; (3.1) 
dt j=0 


where the a,, are non-negative. The choice of the functions y, which maxt- 


N 
mizes J = Sf,’ [s: (x; — v0 | dt, subject to the constraints 0 < y, < x, 
Li=1 


is given by vy, = x,,0 <t < 7), where the 7 are determined by the ¢,,a,,and T, 
and y, = 0,7; <t <T. 

This problem arises in economics in cotinection with maximization of 
profit, as was pointed out to us by H. Markowitz. Many other interesting 
problems arise if we replace the functional f¢’ (x, dG) by the functional 
So Min (x,(t)) dt. We shall discuss these ‘bottleneck’ problems subse- 


i 
quently. 

4. Quadratic Functionals.lf we choose all the functionals occurring 
to be quadratic, measuring the cost, referring to (1.1), by J(f) = 
Se (x — yx — yy dt +a fi" (ff) dt, the problem of minimizing J(f) 
over all f « L°(0, 7) is one which may be solved completely by standard 
variational techniques. Using the known representation theorems express- 
ing the solutions of linear inhomogeneous differential and differential- 
difference equations in terms of the solutions of the homogeneous equa- 
tions, ' * the problem reduces to that of minimizing, in the scalar case, 


Jif) = Se" (alt) + Se’ k(t, f(t) dh)? dt +a, fy’ fe dt, (4.1) 


over all f ¢ L°(0, 7). 

The corresponding problem for difference equations may be included in 
the above by a judicious use of Stieltjes integrals. 

It is now easy to establish 

THeorem 4. [f a(t) ¢ L200, 1), and Jo! R(t, ty) dt; is uniformly bounded 
for O <t < T, there is a unique function f « L?(0, T) which furnishes the 
minimum to J(f). It satisfies the integral equation 


So R(t, t) [a(t)) + mM R(t, t>)f (to) dts| dt, + ayf =: (), (4.2) 
If as in the case of differential equations, 
a(t) = Jy’ R(t, ty)b(t) dt 


for b ¢ L?(0, 1), the value of J(f) for the minimizing f ts a, So! f(tyb(t) dt. 

A similar result holds for the discrete analog. For those problems aris- 
ing from differential equations, the integral equation may be converted 
into a differential equation of twice the original order, with a two-point 
boundary condition. A similar but more complicated result is obtained in 
the case of the differential-difference equations. Since the integral equa- 
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tions which occur are of Fredholm type, it is important to establish an a 
priori existence and uniqueness theorem. 

The related probiem of minimizing Jy’ (a(t) + Jo! R(t, tif(t) dty)? dt 
subject to fj ' f2 dt < ay, may be made to depend upon the above problem, 
under certain assumptions concerning &(f, f;). In passing, let us note the 
following interesting result, which admits of considerable generalization, 

Turorem 5. Let g(r) = Min fy’ (¢ ~ x, ¢ — x) dt, subjectto fy" (f, f) dt 
=y<cr* = Tic, c), andle 


So (en f(t), © — x) dh: = Mr) f(d) (4.3) 


be the equation satisfied by the minimizing function f(t), where \(r) ts a La- 
grange multiplier. Then dg(r)/dr = X(r). 

5. Further Problems.--We have in what has preceded considered only 
some of the simplest representatives of the problems encountered in plan- 
ning over time. Frequently, in applications, there are forcing terms g, 
arising from exogenous causes beyond our control, so that the equation 
takes the form 


dx 


y ~Axtf+e (0) =<. (5.1) 


There are now two points of view that we may take. We may consider 
g to be a random function, with known expected value and auto-correlation 
function, and seek to minimize Exp J/(f, g) = J/(f). In some applications 
g 
this is a reasonable procedure. 
In connection with other applications, it is better to introduce the con- 
cepts of game theory. In this way, arise the problems of determining 


Min Max So’ (x — y,x — y) dt, Max Min So! (x Poe y)dt, (5.1) 


4 & 


: >] 7 : ; 
subject to Jy’ (Uf, f) dt < a, So’ (g, g) dt < as, and many interesting and 
difficult questions concerning games over function space. 


' Bellman, R., “On the Boundedness of Solutions of Non-linear Differential and Dif 
ference Equations, Trans. Am. Math. Soc., 62, 3857-386 (1947 ) 

2 Bellman, R., “On the Boundedness of Solutions of Non-linear Differential- Difference 
Equations,” Ann. Math., 50, 347-355 (1949). See also, “On the Theory of Dynamic 
Programming,” Proc. Nati. AcAp. Scr, 38, 716 719 (1952) 
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FOURIER TRANSFORMS OF TIME SERIES 
By S. BOCHNER 
PRINCETON UNIVERSITY 


Communicated January 19, 1953 
We take all bounded intervals A = A,yg: (@ < x < 8) on the open line 
(— ©, 


o) and we consider complex-valued finitely-additive random in- 
terval functions /( A) for which always 


E\|F(4)|?} < 


= 0. 


(1) 
We call F(A) a Khintchine function if for any two intervals A;: 
As: 


(a, By yy 
(as, 8») and any real t we have 
Ee} F(A,) F(A)! Ie | F(A‘) F(As')} 
where J,, is the translated interval (a@,, + ¢, 8, + 0,m = 1, 2 
it orthogonal if 


2. We call 
I} F(A,) F(As)! Q for AyeAs = 0, 
and orthonormal it furthermore 


iD } F(A) 7! = ci\A| = c(8 — a) 


forall J = A,,, and we then usually put c 


= 1. It follows easily that F(A) 
is orthonormal if and only if it is an orthogonal Khintchine function. 
THEOREM |. (i) /leuristically, that is, subject to criteria, with “any 


Khintchine function F(A) there is associable an orthogonal function G(A), its 
Fourier transform, such that 


; > a en ntBy ws er tte ; 
F(Ags) = ' dG(y) 
‘ 2mriy 


symbolises Stieltjes integration with respect to the interval 
function G(A). 


(2) 
where ‘“dG(y)" 


(11) Conversely, with ‘any’ orthogonal function G(A) there ts associable 
a Khintchine function F(A) such that 


. 


ee 2 rtd. = : 4 : 
Gia.,) = J dF(x), 


~Drix 
and the two transformations (2), (3) are inverses of one another. 


(111) Furthermore, if g(x), W(y) are numerical functions sufficiently small 
at + 0 such that 


g(x) = fe" wy) dy 


(4) 
then (2) implies 
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S. olx — &) dF(t) = f*. 27" Wn) dG(n) 
In particular we have 
SS. ¢(—8)dF(é) = S. Vn) dG(n) (5) 


and, formally, the validity of this relation for sufficiently many pairs {¢, | 
is equivalent with the formula (2) for all intervals (a, 8). 
(iv) Finally, for the random point function 


f, (x) = S ‘a G(x — £) dF(E) 


we formally have 


EMfo (x + 0 fx)} = Si. 2" Wn)? dV (n) 


where 

I (A) = E}\G(A))*} (6) 
ts a non-negative finitely additive interval function which is of bounded vart- 
ation in every finite interval (only). 

Perhaps the best known “‘rigorous”’ case of this ‘‘heuristic’’ theorem arises 
if we take a random point function f(x) which has finite second moment 
and is L.-continuous and theri put 

F(Aas) = Se? f(x) dx. 
This set function is a Khintchine function if and only if the correlation 
Eh f(x + t) f(x)} = R} 
is independent of x (and is continuous in /), and the known theorem is in 


part as follows. 
THEOREM 2. ° A random function f(x) is such a Khintchine point func 


tion if and only if we have 
S-2. 2M dG(y) 
where G(A) is an orthogonal function whose second moment is bounded, 
P(A) = EfjG(A)|?} < R(O). 
We then also have 
R(t) = fe" dT (n) 
and G(A) is uniquely determined if it is L»-continuous from the left, say. 
Another rigorous case, a seemingly different one, is as follows. 
THEOREM 3.4 An orthonormal function F(A) has an orthonormal Fourier 


transform G(A), and relation (5) can then be validated for pairs of functions 
lo, W} both in Lo(— ©, ©), 
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There are other classes of functions } F(A)! which are invariant under 
Fourier transformations, as for instance the class of functions which are 
Gaussian in the sense that for any disjoint intervals 4), ..., A, the joint 
distribution function of the 8 random variables 


F(A), ..., F(Ax) (7) 


is a Gaussian one. 

Such an F(A) need not be a Khintchine function, but if the characteristic 
function of the joint distribution of (7) has the special form 

exp [—u(| Ajay? +... + (Ag) ax?) | (S) 
for some fixed u > 0 then F(A) is even orthonormal, and it 1s, in fact, a so- 
called Wiener-process. Now, the family of such | F(A)} is again invariant 
under Fourier transformation. As a matter of fact if we take a fixed aver- 
aging weight dp(u) > Gin 0 <u < © with J,” dp (u) = 1 and if we re- 
place (8) by 
So" exp {—u(j|4i} a? +... 4+ (Ay! a?) ] dp(u) 

then such a class | F(A)! is invariant likewise.! 

We are now going to bring Theorems 2 and 3 together by means of a 
statement which is the main remark of the present paper. 

THeoreM 4. [f | fe(x)}, 0 < € < @ is a family of Khintchine point func- 
tions which is such that for any € > 0, €' > 0 we have 


gtd ¢ 5 Tae ti | w/e(x—§&)? 4 
; “) = Fi f (&) dé 


then there exists an orthogonal function G(A) such that 


fizpu Je" eo Gly) 


and 


Et ft (x + 0) fe(x)} = Ji er e~ "dT (y); 


meaning for instance that if we put 


G* (4.3) = J,’ e~*™” dG(y), 
then we have 
f (x) = S fo" iC (y) 


in the sense of Theorem 2. 


Now, if we take a Khintchine function f(x) to which Theorem 2 applies 


Poca ei 
f* (x) = | =~ eW 8/ts —O* H( 8) de (12) 
» Ve 


and form 
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then we obtain a family to which Theorem 4 applies, and the spectral func- 
tions G(A), P(A) in Theorem 4 are the same as those in Theorem 2. 

On the other hand if we take an orthonormal F(A) to which Theorem 3 
applies then the family 


om - 4 re ee 
f (x) = / nig “Ine are 
rn VE 


falls under Theorem 4, and the spectral functions G(A), ['(A) of Theorem 
4 are now those of Theorem 8, and thus Theorems 2 and 3 have been 
brought together rather closely. 

Notwithstanding, this, if we take an arbitrary family |/‘(x)} to which 
Theorem 4 applies, then it can be represented as an integral (12) by means 
of a Khintchine point function f(x) if and only if '(A) has bounded variation 
in (— ©, @), 


S. dv (y) < R(O). 


Now, for P(A) = ¢ Al, ¢ > O, this variation is always unbounded, and 
thus, in their proper versions, Theorems 2 and 3 appear to be very different 
indeed, nevertheless. 
The distance between Theorems 2 and 3 can also be stated in certain 
quantitative terms. If }y, ¥} are connected by (4) then for the function 
fo (x) = JL o(x — £&) f* (&) dé, 
we obtain from (10), (11) the relations 


fixe) = Sf. &™-" Wy) dG(y) 


EVE (xe tafe (x)} = fL% 77" |x(y)|? dl (y), 
and if G(A) and y¥(y) are such that 
S—. Wly))? dV(y) < @ (13) 
then the limit 


lim,.o f¢ (x) (14) 


exists in an appropriate sense. Now, if (4) has bounded total variation, 
then (13) holds for y(y) = 1 and this being the Fourier transform of the 
Dirac function, the limit (14) is then a point function f(x) as in Theorem 2. 
However for '(A) = ¢ A) we have, say 


3 “| 
/ + | , aly) < @, (15) 
: y? 


l 
and thus (13) holds for Y(y) = 1/7y sin 2ray, this being the transform of 
the function ¢(x) which is = 1 in |x| < a, and 0 outside, and this then leads 
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to a boundary set function #(A) asin Theorem 3. In general, if for a given 
family f‘(«) there exists an exponent k for which 
3, Aly) < 
Ms 

then the limit (14) exists as a kth order set function as was introduced by 
this author in the theory of generalized Fourier integrals*® (and later carried 
further in the work of Laurent Schwartz). But this kind of generalization 
would be worth undertaking only if requirements of stationarity and 
orthogonality would not be introduced from the very beginning as we have 
done for the present, and problems of smoothing and inversion of smoothing 
would be approached as well, as we will do on another occasion. 

Returning to Theorem 4 we note that the semigroup-requirement (9) 
could be replaced by other ones, as for instance by 


e's f' (&) dé 


+ (x £)° 


errixy eri! dG(y) 


instead of (10). The two-variable function f(x) 1s then harmonic in the 
complex variable x + 7 €in the half-plane « > 0, with f(x), if existing, being 
then its boundary value, and this set-up certainly has its attractions for 
further studies. But the function f‘(«) of Theorem 4 is a random-valued 
solution of the diffusion equation 


O°f*(x) Of'(x) 
ie — 


Ox” Oe 


and stochastically this kind of approximation might be more appealing 
perhaps. 

Finally, to complete the picture, we will state the integral-time-counter 
part to Theorem | and we note that to a large extent it will be “mgorous” 
from the very beginning. 

THEOREM 5.''* A (two-way) sequence |f,\ of random variables is a Khint- 
chine sequence if and only if there exists an orthogonal function G(A) on the 


closed circle 


such that 
fo = Se” dGy) 


Dually, if GCA) is a Khintchine function on (16) then (17) defines an orthog 
onal sequence and we have 
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S\ e dR 
where 
Rit ke NG( Ag t t) G( Ap). 


Conversely if an orthogonal sequence \f,{ is such that for every interval in 
(16) the series 
, st Bete e 
G(A., S 
Lo “ 
Jrin 


converges suitably, then this G(A) is a Khintchine set function for which (17 
holds. 
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ON THI BASIS THEOREM FOR FINITE ABELIAN GROUPS: 1V* 
SY JESSE DouGLAS 
CocumBia Universiry, New York 
Communicated February 1, 1958 


|. The existence of a basis for any finitely generated vector space 1s 
immediate; one has only to select a minimal subset of generators. The 


Irief proof depends essentially on the fact that the domain of sealars ts a 


held; this permits solution of any linear relation among vectors for any one 
of these vectors whose coetlicient 40. 
Why is more required to establish a basis for any finite abelian group Gr? 
Fundamentaliy, it is because the corresponding module .\/ (G, A) is, in 
general, only that and not a vector space. Here the “external operation 
is 06 — a@ (additive notation), and the operator domain A consists of the 
integers a considered modulo m, where m denotes the maximum period! of 
the elements 4 of G. 

Thus A ts, in general, only a ring, nota field. However, in the particular 





308 MATHEMATICS: J. DOUGLAS Proc. N. A. S 


case where m is a prime Pp, A is indeed a field, the abelian group becomes a 
finite vector space, and the proof of the basis theorem runs very simply as 
follows. 

By hypothesis, p@ = O for every element 60f G. Let S = (A, 6, .. ., 4) 
be a minimal set of generators of G. Then (we assert) this set must be 
independent, and so constitutes, according to definition, a basis of G. “‘In- 
dependent”’ means that in any linear relation 


a); + AA» + > +S + a6, = 0 


-ach term a; = 0; in the present case this is equivalent to a; = 0 (mod p) 
fori = 1,2,...,7r. Otherwise, x exists such that xa; = 1 (mod p) for some 


a, (use of the field property, existence of a reciprocal); then, by multipli- 
‘ation with x, 


X00, +... Ht X0;-30;-1 +O, + x0;4:0;4; +... + x00, = 0; 


thus @, is a linear combination of the other 6’s, and the proper subset .S — 6, 
of S is still a system of generators—but this is contrary to the supposed 
minimality of S. 

The proof offered in the present note consists in a slight adaptation of 
this argument for the prime period case so as to afford an inductive proce- 
dure applying to the prime-power period case, which is essentially the 
general case. 

2. Let p denote any prime and X\ any positive integer. We shall say 
that a finite abelian group G is “‘of p*-type”’ or is a “p*-group”’ if every 
element 6 of G obeys the equation 


p’o = 0. (1) 


A primary abelian group is one which is of p*-type for some p and  . 
Because of the standard fact that every finite abelian group is a direct prod- 
uct of primary groups, we may restrict the proof of the basis theorem to 


primary groups. 

3. Turorem. very primary abelian group has a basis. 

The proof is by induction on X. For \ = 1, the theorem has been proved 
in the introductory section, 

Assuming, then, that the theorem is true for all groups of p* '-type 
(\ > 1), we shall infer its truth for the p-type. 

Let G, denote any p-group. Consider the subgroup G,_, of elements 
pé where @ ranges over Gy. Then by definition (1) we have p*~'(p0) = 0; 
hence G,_, is a p*~'-group. By hypothesis, therefore, G,_, has a basis 


B = (p6;, p02, .. ., POm). 


To the set of elements 6; whose p-multiples form B adjoin all elements 
of G, which obey the equation 
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pw == {), 
In this way we form the set? 
» = (A, As, ee sy Oui Why Way . + oy Wy). 

We assert that S 7s a system of generators forG,. Tosee this, let @ denote 
any element of G,; then pd belongs to G,_, and is therefore expressible in 
the form 

po = a, Pb; + Ap. + cee + AmPOm.- 
Hence 
Pid — G10, — Ade — ... — Abn) = YU; 
thus the element represented by the parentheses obeys (2) and is accord- 
ingly one of the w's: 
@ — A)\0; — ABe — 2... — Anbm = Wy 
so that 
o= 0; + 8» +- 4 + Am m + Wt, 
which establishes our assertion. 

1. Now, the set S, being finite, must contain a minimal subset So of 

generators of G,; by rearranging notation we may suppose 
So = (O;, A, ow ey &. WW, Way « « oy We). 

We assert that this set is independent, and therefore constitutes a basis of 
G). 

Indeed, let 


ay), + or + a, + bya + © ae + bw, = (0 (33) 


be any linear relation among the elements of So. 

Then each coefficient a,, b; must be divisible by p. For suppose, e.g., 4a; 
non-divisible by p. Then an integer x exists such that xa, = | (mod p’). 
Multiplying (3) by x and taking account of (1), we get 


xa,0, +... xa;_,0;-, + 0 + x0;4,0;4,... + xa,0, 
4 xbywy + a as 4 xh A, = (). 


Hence @,, expressible as a linear combination of other elements of So, may 
be omitted without loss of the generative character of the set- but this is 


contrary to the minimal assumption regarding So. 
Each coefficient in (3) being therefore divisible by p, this relation takes 
the form 


CP), + ig oe ora + cp, + d\ pw + her + d pws -_ (). 
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Here each w-term = 0 because of (2), and then each 6-term = 0 because of 
the independence of the elements p6,, . . ., p6,, which are part of the basis 
B of G,_,. The proof of our main theorem is now complete. 


* Work on this paper supported in part by contract with the Office of Ordnance Re- 
search, U. S. Army. Previously published papers by the author on the same topic: 
these PROCEEDINGS, 37, 359-362, 525-528, 611-614 (1951). 

' Equal to the l.c.m. of the periods, since, in an abelian group, the period of any ele 
ment is a divisor of the maximum period 

* The basis B may be assumed not to contain the identity, 1.e., p6, # 0 for all 7; hence 
no @is equal to an w, and the elements of S are all different. 


GAMES WITH EQUILIBRIUM POINTS* 
By RICHARD OTTER AND JOHN DUNNE 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF NOTRE DAME 
Communicated by J. von Neumann, February 5, 19538 


This paper deals with the combinatorial aspects of games of strategy, in 
particular with the relation between the information partition and the set 
of strategies. In his book J. von Neumann! has shown that a zero-sum two- 
person game is strictly determined (s.d.) provided there is perfect informa- 
tion. Utilizing the notion of equilibrium point introduced by Nash,’ 
Kuhn* has generalized this theorem to the general finite game. On the 
other hand, it is easy to see that there are less restrictive conditions on the 
information partition which allow a game to be s.d. regardless what the 
pay-off functions are. In this paper we state a theorem which character- 
izes these games structurally. Only the proof of sufficiency is given; that of 
necessity, being a good deal simpler and less informative, is omitted. 

A rooted tree is a finite, non-empty, partially ordered set 7° in which (7) 
there is an element 7 satisfying r < v for every v € 7 and (11) for every ve 
7 — {r{ there is only one descending chain v = v, > 2,1 >... >% =r 
such that v, and v,, are adjacent in 7 fort = 1,2,...,”. Wecall the ele- 
ments of 7 vertices, the element 7 mentioned in (1) the root of T and read 
the relation « < v as ‘uw is lower than v."" For each v « 7, let A(v) denote 
the set of elements in 7 adjacent to and higher than v, and let 6(v) denote 


the number of elements in A(v). The family @® = (8(7)t > 0), is a parti- 


tion of the set 7. We put & = 3(0) and call its elements end-points and 
put B = 7 — Eand call its elements branch- points. 

Suppose now that for each v e B we select a fixed enumeration of the 
set A(v), 1e., A(v) = fu, wo, ..., Ug}. We express the fact that u; is the 
ith vertex attached to v by writing uw, = (v, 7). If we now express v in the 
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same way in terms of the vertex immediately below it and proceed in this 
way down to the root 7 to which we assign the symbol (0), we get that 
each vertex « may be expressed in the form 1 (O, 41, 42, ---,%,) and we 
sav we have a coordinate system in 7°. 

In order to define a game we follow McKinsey! in requiring that there be 
given (2) a rooted tree 7 with a coordinate system (17) the player partition 
» = (P,)po, of the set B, (271) the information partition % of the set B, 
(av) a probability measure wy) on Ly, the set of strategies for chance, and (v) 
afamily y = (g))g-, of real valued functions defined on /. We call the 
quadruple (7, Y, 9, wo) a game structure and the pair (I, y) we calla 
game. 

We consider the indices k = 0, 1, 2,...in | 1 correspondence with a set 
of players, the player corresponding to index 0 is called chance, the rest 
personal players. 

The information partition ‘> ts required to be a refinement of the parti 
tion Y and of the trace of (& on B. We also require of © that if w < v then 
wand v are notin the same set from‘). The equivalence relation that exists 
between two vertices “u,v when they belong to the same information set D 
« Dis written as uw == 7, and sometimes as u =v (k), provided D ¢ P,. 

In order to explain what is meant by 2, mentioned above we let Y con 
sist of all those mappings o of B into 7 such that 

(a) a(v) € A(z) 

(b) at ofa) (u, 1), o(v) v, J) and u v then 7 j. We con 


sider © as subset of I] A(z) and observe that it is determined when 7), @, “ 
reB 
are given. Its elements we term sfrafegies. The set = pry, (2) ts called 
the set of strategies for player k and we evidently have = Hi =,.° We 
ko 
may write for oe X, (ov, O1, O2,...) (op)p-» Where a, prea). It 
should be pointed out that the previously given definitions of a game have 


suffered from an inadequate description of the probabilities involved. — All 


that has been supposed given 1s polo) ao(v) v,2)) for each v « Py and 
1=1,2 Biv). This data alone does not suffice to determine uy which 


9 —yp et 2 


one needs to know in order to compute the expected pay-off. 
For each o ¢ Y the last term of the sequence 7, a(r), o(a(r)),... which is 
meaningful is an end-point m(a). The correspondence o — x(a) deter 


mines the mapping wof Yon /:. We then put 
Pet Gt, Oh. <6 A £,( lav, 1, G2, ...)) dy ( oy 
and call these the expected pay-off functions. We puttorve7 
S(v) iaeLivs x(a)j 


S,(v) pre, (S(v) 
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The following lemmas are very useful in the theory. Their proofs are 
simple and we omit them. 
LeMMA |. (a) We have o « S(v) tf and only if u < vimplies a(u) < v. 
(b) We have o, € S,(v) if and only if u <vand ue P, implies o,(u) < 2. 
LEMMA 2. S(v) = II S,(v). 
k 


"0 

LEMMA 3. (a) S,(u) A S,(v) = pif and only tf for some u,v) € Py with 
My v, we have (mu, 1) < u and (wy, J) < v where i # j. (b) S,(u) ¢ 
S,(v) of and only if for each v1, € P, with 1, < v there is au, ¢ P, with uy = % 
and such that both (uy, 1) < wand (2%, 1) < v for some 1. 

For o, € X, let us define the “track of 0,” by tr(o,) = joe 7 oy € Si(v)} 
then the previously given definitions of equivalent strategies amount 
to this: o, equivalent to o,’ if and only if tr(o,) = tr(o,"). We may also 
call two information patterns , 0’ (for the same 7’, ®) equivalent whenever 
tr(Z;) tr(Y,’) for k = 0, 1, 2,... where &,, 2,’ are determined by (7, 
y, DB) and (7, W, 2’) respectively. 

We calla set C ¢ P, U Py, k # 0, connected for player k if any two of its 
vertices may be connected; that is, if u, v « C then there is a sequence “= 
ly, Wy, ..., Uy, = v such that wu, e« C and either 4; adjacent “4; or u, 
u,-1 (k) and S,(u,) A Si(ui31) ¥ @fort = 1,2,..., 0. 

Aset A ¢c P, UY Py, (k # 0), is an info component for player k means 
that (2) A is connected, (77) its minimal vertices belong to P;, and (111) K 1s 
maxinial, t.e.,1f A’ > A and satisfies conditions (2) and (72) then A’ = A. 

LemMMa 4. For any game structure 1 = (7, &, SD, wo) there is an informa- 
tion partition ‘)' equivalent to all of whose sets are connected for the player 
to whom they belong. 

Proof: We put u =’ v (k) if and only if « =v (k) and u andv may be 
connected, Evidently the partition ’ defined thereby is a refinement of 
, hence tr(X,) ¢ tr(2,’). Suppose that tr(o,’) is not an element. of 
tr(X,) then there exist uv, v with u == v (k) and with (u, 2), (v, 7) € troy’) 


but 7 # j, hence uv ==‘ vis not true. If S,(uw) 9 S;,(v) = @ we may replace 


u,v by the w, 7, of Lemma (8a). By working down toward the root in this 
way we see we may suppose S,(“) N S,(v) # @ but this means u ==" v which 
is a contradiction. Thus ‘D’ is equivalent to and it is clear that its sets 
are connected. 

TueoremM. Let VT = (7, &, D, wo) be a game structure where By = 6 and 
uy(So(v)) > O for everyu eT. In order that the game (V, y) possess an equilib 
rium point for every family y of payoff functions, it ts necessary and sufi 
cient that whenever for some k > O we have u == v (k) and S,(u) A S,(v) # o 
we also have S;(u) = S,(v) for any 7 ¥ 0, k. 

Before proceeding to the proof of the sufficiency we shall establish a 
lemma: 

LemMa 5. Under the condition of the theorem the info components have 
the following properties: 
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(a) If J, K are distinct info components belonging to personal players 
j,k, respectively, and tf there is a vertex v « J higher than some vertex v1, «eK A 
P, then there is an integer 1 such that for every vertex ue J we have u > (uy, 1% 

£ i : 


for some u, == (Rk). 
(b) If u <v and both lie in some info component J of player j then so do 


all vertices between u and v. 


(c) If M ts the set of minimal vertices of an info component then M 1 tr(a,) 
k>1 


ts either empty or itis Mand M  tr(oy) ts either empty or it ts a single ele 
ment of AM. 

Proof: (a) We first suppose that 7 # k. Now the two vertices u, ve J 
may be connected in J. u = Wy, Wi... Wr v. If w, wij) and 
S(w) A S(w,.) # @ then by the condition of the theorem S,(w;) 
S.(wei). If w, is adjacent to w,;, we have S,(w,) = S,(w;1.) because 
Ww), W,;-, are not elements of P,. Now since 7 € P;, and m4 < v say (v%, 2) 
< v then by Lemma (34) there is a v, (Rk) such that (am, 7) <u. 
The case } = k follows from this and part (4) by similar arguments. 

(>) Put W = jweTiu << w <j then Wo ¢ P, MN Po. For suppose 
ve Won P,, k # O, j, and let A be the info component in which 2, lies. 
Then J, A satisfy the hypotheses of part (a) so we find there is a ™ such 
thatv >v, >u>u,and2, u, (Rk). This contradicts the definition of an 
information pattern, so I’ ¢ P,; VY Py. Since J is a maximal connected 
set we have IV’ ¢ J. 

(c) This follows directly from part (a) and Lemma 1. 

Proof of Theorem (Sufficiency): Since the question of whether or not 
(o1, 0, ...) is an equilibrium point may be settled by looking at tr(a;), 
tr(a-),..., We may suppose by Lemma 4 that the info sets for each player 
are connected for him. We proceed by induction on the number of info 
components in the game. If there is just one info component belonging, 
say, to player 1, then fork > 1, P, = @ so X, consists of a single element 
a,. So if a, maximizes fi(o;, 0, ... ) then (a), 0, ...) is an equilibrrum 
joint. The case of no info components is even simpler. 

Suppose now there are V > | info components and that the theorem is 
true for games with less than .V info components. Using Lemma (5a) 
we easily construct an info component A’ such that the set /:’ = jue 7 
K’ u>vforsomev eA‘! is subset of /¢. Let 1/ denote the set of vertices 
minimal in A’, then y is not an element of /. For indeed any other ver 
tex is higher than ry and soit yr ¢ J/ ¢ A’ then 7 K’ V FE’ which contra 
dicts V > 1. If we now construct 7” = A’ VU k’ UV fr} then 7” isa tree 
whose root is r (Lemma (5))) with /:’ for its set of end-points. We con- 
struct a game (I'’, y’) = (7", 0, 9’, wo’, y’) by relativizing the original 
game to 7” and by assigning r to chance. Thus ®&’ is the trace of Pon AK’ 
except that Po’ contains 7, ‘0’ is the trace of % on A’ together with }r{. 
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We put 05 Prpmk: (Xo) then XY)’ = Oy K Al is the set of strategies for 
chance in I’. We put for ay’ € Xo’, 
Me pr pags (Bo) A Sy(u)) 


by (au) = po’ (Oo, W) 
po( SoC. A) ) 


This is a probability measure on Ly’ because by Lemma (5¢) (Sy(w))yey 18 
a partition of Sj(.MW/). And finally we put yy’ = (g,")e>, where g,’ is the 
restriction of g, tok’. Evidently, the game (I'’, y’) has just one info com 
ponent, namely A’. Let (a,', o’,... ) be an equilibrium point. 

Now put £” = (fF — FF’) v MV, A” = B — K’and finally 7” = A” vu 
fe”, then 7” isa tree with yas roc’ and /:” as its end-points. We construct 
a game (I, 7”) by letting &”, “2” be the trace of &, ‘D, respectively, on 
A”. Then Yy" = prp,(X») is the set of strategies for chance in ['” and we 


put wo"(o0") = wol(pr’) p,o0")). Finally, we put y” = (g,"),~, where 


s g,(e”) fe"ek — E' 
C9, tay’, Ges oS) ife”« M 


The game (I',” 7") satisfies the hypothesis of the theorem and has .V-1 info 
components so there is an equilibrium point (0,", 02", ... ). 
We identify the pair (¢,’, ¢,”) with a strategy a, « X, by putting a, = 


é,” fork > 2and 


fa,'(v for ve P,’ = Kk’ n P, 
la,"(v forve P” = (B — K’) nP, 


a\(v) 


That o, « >, fork > 21s obvious. For k = |, we ebserve that if a,(%) = 
(u, 1) and o\(v) = (v, 7) where u =v (1) then because of the connected 
ness of the information sets either u,v are both in A’ or both in A” soz = /. 

It is now an elementary computation to show that (o;, o2,...) 1s an equi 
librium point for (1, y). 

CoROLLARY. T/f in addition Py = @ then in order that (V, y) possess an 
equilibrium point for every family y it is necessary and sufficient that u # 
vand u==v (k) imply S,(u) A S,(v) = @. Hence the information partition 
iy equivalent to one with 1-element sets, 1.e., perfect information. 


* This work was supported in part by the Office of Naval Research, in part by the 
National Bureau of Standards 

' von Neumann, J., and Morgenstern, O., Theory of Games and Economic Behavior, 
Princeton University Press, 1944, pp. 112-128 

* Nash, J. F., these PROCEEDINGS, 36, 48-49 (1950 

‘Kuhn, H. W., /bid., 36, 570-576 (1950) 

' McKinsey, J., /ntroduction to the Theory of Games, McGraw-Hill, 1952, pp. 119-120 

’ We are following the customary conventions for product sets, e.g., if Py = @ then 
Y, is taken as some |-element set. The notation is for the most part that of Bourbaki. 
N., Theorie des Ensembles, Hermann, Paris, 19390 
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ON THE NUMBER OF HUSIMI TREES, I 

By FRANK HARARY AND GEORGE E. UHLENBECK 
DEPARTMENTS OF MATHEMATICS AND Puysics, UNIVERSITY OF MICHIGAN 

Communicated by J. von Neumann, December 31, 1952 


1. Introduction. It was pointed out by one of us! that the theory of 
the so-called virial development of the equation of state of a non-ideal 
gas, leads to combinatorial problems for certain types of linear graphs, 
which can be considered as direct generalizations of the well-known 
“trees” of Cayley. We have called the simplest such types of graphs 
Husimt trees since they were first introduced in a paper by Husimi.* A 
Husimi tree is a connected graph in which no line les on more than one 
cycle. It is characterized by the numbers m, my, ny of lines, triangles, 
quadrilaterals, ete., out of which it is built up. A pure Husimi tree is one 
which consists of only one type of figure. Otherwise we speak of mixed 
Husimi trees. A pure Husimi tree consisting only of lines 1s a Cayley tree. 
If it consists only of triangles we will call it a cactus. 

There are two kinds of combinatorial problems connected with each 
type of graph. The first and simplest kind of problem is to find the 
number of topologically different graphs if a// points are distinguished from 
each other. For Cayley trees of p points the answer* is simply p” 
For Husimi trees of p points and of type me, ms, ,n, the answer? is: 


p! 


Dp Pp: 
HW ((@ — 1)!]"n,! 
which clearly specializes to p” * for Cayley trees. The second and more 
difficult kind of problem is to find the number of topologically different 
graphs if the points are not distinguishable.‘ For Cayley trees the problem 
has been solved® in all details. The method has always been to start 
with finding the number of “rooted” trees, in which one point is preferred, 
and then to enumerate the different types of “free” trees in terms of the 
different rooted trees. We have generalized this procedure to Husimi 
trees, and the main purpose of this note is to present an outline of the 
method and some particular results. 
2. The Number of Rooted Husimi Trees. To simplify the writing we 


will exemplify the method for the case of cacti. Let A(xv) = SCA,x" be 
0 


the counting series (generating function) for the number A, of rooted cacti 
consisting of m triangles (A) = 1). Starting from the root we can classify 
the different cacti according to the number ¢ of principal leaves, which are 
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the triangles having the root point in common, Let A(x) be the counting 
series of the cacti with r main leaves then obviously: 


A(x) = 1 + De A™ (x). 
r 1 


To obtain A(x) one has to find in terms of A(x) the counting series (x) 
of an unordered pair of rooted cacti since from the two remaining points of 
the single triangle, which starts out from the root, any combination of two 
rooted cacti can emerge. Although &(x) can be found very simply, it is 
useful to already apply here the notation and results of the Hauptsatz of 
Pélya’s paper (Ref. 5, p. 163). Let Y be a permutation group; then we 
denote by Z(\) the so-called cycle index,® and by Z(\, /(x)) the function 
which results if one replaces f, by /(x*). Let S, denote the symmetric 
group of degree n. By an immediate application of the Hauptsatz, we 


obtain: 


Side x - : 
A(x) xt(x) = xZ(Go, A(x)) = = [A(x?) + A(x) ]. (1) 
The extra factor x in (1) takes account of the single main leaf triangle. 
A second application of Pélya’s Hauptsatz then gives A(x). One has 
to choose at each of r places a pair of rooted cacti, and therefore: 
A(x) WZ(S,, E(x)). (2) 
The factor x” takes account of the r main leaves. Clearly (1) is a special 
case of (2). By a simple calculation, one finally obtains: 


A(x) = 1+ Yi wZ(S,, &(x)) 


r 1 
exp | So facet) + ate} 
X] 93" ! 
as the functional equation for A(x). 
The generalization to rooted pure Husimi trees consisting of m-gons is 
accomplished by redefining the function £(«) in equation (3) appropriately. 
If 7\"(x) is their counting series, then the functional equation is 


— 
S 


8 


; 7 
Em u(x") | where 


| & "| 
(x) = = (7 (x))"-! + (T™(x*))>— * (T™(x)) 


< 
Sm 


The extension to mixed Husimi trees is very simple. It is necessary to 
introduce counting series with as many variables as there are kinds of 
polygons in the tree. For instance, let 7(m», 3) be the number of rooted 
trees consisting of m2 lines and n; triangles; one then has for the counting 


series: 














Vol. 39, 19538 MATHEMATICS: HARARY AND UHLENBECK 


T(x, y) a (ns, n3)x"") 
0 


ne, Wa 


and the functional equation: 


: — |x Ys, ee 
T(x, y) = exp >> | Pe Fa } T(x, y*) 4+ rv. | (5) 
r=1Ul ~! 
3. The Dissimilarity Characteristic. To count the number of free 
Husimi trees we will generalize the method which Otter’ developed for 
the case of Cayley trees, and which makes use of the notion of the dis- 
similarity characteristic. Consider a Husimi tree // with points belonging 
to different sets, say of different color.’ The notion of similarity for 
points, lines or cycles, will be defined in terms of those automorphisms of 
the tree //, which also preserve the color of the points. Call p*, /*, c* 
the number of dissimilar points, lines and cycles in //. Let a be the 
number of symmetry axes of //, which do not pass through a cycle. Clearly 
a can only be zero or one. By a symmetry of a cycle we will always mean 
a symmetry of the whole tree with respect to the cycle. Let c;, and ¢; 
be the number of dissimilar cycles, which have at least one point to point, 
respectively line to line symmetry axis. Then one can show: 


1 = p* — (I* — a) + (c* — ch + cy). (6) 


For Cayley trees this reduces to Otter’s theorem | p* — (l* — a)’ 
which can also be expressed in the form 
* . 

1 = p* —hj, (Ha) 
where /,; is the number of dissimilar lines having no symmetry. Analo- 
gously for pure Husimi trees consisting of m-gons one can express the 
dissimilarity characteristic in terms of the numbers of dissimilar m-gons 
which have a certain type of symmetry. For instance for cacti there are 
the three symmetry types: 


e ° * * * ~ 
where the letters refer to rooted cacti. Let ¢,, ¢:, ¢; be the number of 
dissimilar triangles of the three respective types of symmetry occurring 
in a given cactus //. Then clearly 


I* = 3¢ 


* » 


tf = cf oe 


+ 2c, + ¢. 
? + Cs 


* 
J 
* 
1 
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” * . - . . 
and sincea = Oandc,, = ¢, = V (as is the case for all pure Husimi trees, 


consisting of odd polygons) one obtains from (6): 


om od e 
Ll = p* — 2, — Cy. (6b) 


For a pure Husimi tree consisting of quadrilaterals one has five types of 
symmetry, illustrated by: 








A 8 Ay- 
























































. * * . - . . . . ‘a 
Lt cy; , ¢, denote again the number of dissimilar quadrilaterals of these 


symmetry types occurring in a given tree, then one has 
* 

= 4c, + 2c, 
* 

= ot ¢ 


( 


ul 
Thus one gets from (6), since a is of course again zero: 
. * . 
| = p* — Bcf — 268 ~ ct + (6e) 
In this way one can go on and we have obtained for every pure Husimi 
tree consisting of m-gons a formula of the form: 
(m) 
= | * > 
p™ - = aC, (6d) 
' i 
where s(m) is the number of symmetries of an m-gon and the a; are integers 
expressible as simple number-theoretic functions of m. 
4. Functional Equations for the Number of Free Ilusimi Trees. Let 


im) =, m 
(™(x) = 1% ag 
" 0 


be the counting series for the number /," of free Husimi trees consisting 
of nm m-gons. The problem is to find a functional equation for £C"'(x) in 
terms of the counting series 7°”'(v) for the rooted trees. Now it is clear 
that by summing the dissimilarity characteristic, written in the form 
(Gd), over all pure Husimi trees of 1 m-gons, one gets: 

aim) 


= St ” » 2 a; V m (7) 


tr 
1 


where the },", is the total number of dissimilar m-gons of the 7th sym 


Lia 
metry type in the collection of Husimi trees of m m-gons. From (7) one 
then gets: 
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(x) = T™(x) — YO a, Yi"(x) (8) 
' 1 
where the ¥,"(x) are the counting series for the ¥,",. 

By using the methods of Pélya® one can find in any specific case the 
counting series },"(x) in terms of 7‘"(x) by inspection. For instance, 
since in a Cayley tree one gets a line with no symmetry by hanging on its 
two end-points two different rooted Cayley trees, one gets for Cayley trees 
(f°? (x) = tx); T(x) = T(x)) from (6a): 

t(x) = T(x) — Y¥,(x) 
with:!° 
= x /Z(%, T(x)) — Z(Sx, T(x))] 


= = (T(x) — T(x*)], 
which leads to the result of Otter: 
U(x) = T(x) — = [T*(x) -- T(x?)]. (9a) 
Analogously for cacti, the counting series for triangles with no symmetry 
is given by: 
Y}?(x) = x [Z(M%s, T(x)) — Z(Ss, T (x)) | 
eee ‘ 
7-4 [A%(x) — 3A(x) A(x?) + 2A(x) | 
) 
where we have written A(x) for 7‘*(x). Since triangles of the second 
symmetry type must have two isomorphic rooted trees on two end-points 


which are different from the tree on the third end-point, one sees easily 
that: 


¥3° (x) = x[A(x) A(x?) — A(x) |. 


Writing 6(x) for ¢? (x), one gets therefore from (6b): 
: , 
6(x) = A(x) — 3 [A%(x) — A(x’) |. (9b) 


Turning now to pure Husimi trees consisting of quadrilaterals, we have 
to consider the symmetry types shown in Fig. 2. Clearly, any two differ 
ent pairs AA and BB of rooted trees will lead to a quadrilateral of the third 
and fourth symmetry type. Therefore: 


V(x) = V(x) = x[Z(X, Ole?) — 2(Sz, O(x?))] 


= * [Q%(x*) — Q(x*)] 
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where again’ we have replaced for simplicity 7°"(x) by Q(x). For the 
second symmetry type, one has to count the number of arrangements of 
two different trees A and C together with a pair B, B, where B is “‘inde- 
pendent’’® of A and C. Therefore: 


¥,°(x) = xQ(x?)[Z(%, Q(x)) — Z(S, Q(x))] 
ee . Bes 
= 51Q™)O(x*) — Q*(x*) ]. 


Finally for the first symmetry type one has to begin with the number of 
different quadruples of rooted trees A, B, C, D. This is counted by 
Z(My, Olx)) — Z(Ss, Q(x)). This last expression must then be multiplied 
by 3 since the four different rooted trees can be arranged in three non- 
equivalent ways at the four vertices of a quadrilateral. However this is 
not enough. If two of the rooted trees are identical one can still arrange 
them in a completely unsymmetric fashion on the corners of a quadrilateral. 
We therefore have to add the number of arrangements of two different 
trees (say B, C) together with a pair (say A, A) of identical trees. This 
is counted just as for the second symmetry type by 


Q(x?) [Z(%, Q(x)) — Z(S2, O(x))]. 


But now we have too much, because A can not be equal to B or C, since 
otherwise the quadrilateral would have a symmetry. So we must subtract 
the number of arrangements A, A, A, C where A is different from C. This 
is clearly counted by Q(x°)Q(x) — Q(x‘), and so we get altogether: 


Y,;“(x) = : [O4(x) — 207(x)QO(x?) — O7(x?) + 20(x*) J. 
From equation (6c) one then gets the final formula (1° (x) = q(x)): 
q(x) = Q(x) — : [3Q%(x) + 207(x)Q(x?) — 80%(x7) — 20(x4)] (Ye) 
In this way one can go on, but the details become very involved. Using 
the same method we have also derived the equations for pure Husimi 


trees consisting of pentagons and hexagons, and the results (after correct- 
ing an error pointed out by R. Z. Norman) are: 


9 
p(x) = P(x) — = x[P5(x) — P(x*)] (9d) 
») 


h(x) = H(x) — _ [SIT8(x) + 37 (x IP ( x?) — 4IP (x?) — 27P(x*) — 2H (x) 


(9e) 
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From these results one can further easily obtain a formula for the 
counting of mixed Husimi trees. In this case, one easily sees that the 
dissimilarity characteristic (6) can be written in the form: 


1=p*-DYdci. (10) 


where ¢) , is the number of dissimilar cycles of k points and the ith sym- 
metry type (for k = 2, this refers to lines not contained in a cycle). There- 
fore one can combine the formula (9) in an additive fashion except that 
now one has to use counting series with more variables. For instance, 
for free Husimi trees consisting of lines and triangles one finds for the 
counting series /(x, y) the functional equation: 


(x, y) = T(x, vy) - : [T?(x, y) — T(x’, y*?)] — 5 [T3(x, vy) — T(x, y*)] 


where 7 (x, y) is determined by (5). 

In a subsequent paper we hope to present a more systematic treatment 
of the general case of pure trees consisting of m-gons, and an asymptotic 
evaluation of their numbers. 


1 Uhlenbeck, G. E., ‘“‘Some Basic Problems of Statistical Mechanies,’’ Gibbs 
Lecture, Am. Math. Soc. (1950); (in process of publication). Compare also Riddell, 
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6. The cycle index of a permutation group is defined (see Polya, reference 5, p. 159) 


as follows: 


l , 
Z(%) = j S hj, de, tol te? ‘ fi 
nS 


where / is the order and s the degree of XY; fi, fo, , fs are s variables; Aj, . . , fy is the 
number of permutations of ¥ containing j; cycles of length 1, j, of length 2, ete. The 
round summation sign is over all s-tuples (j,) such that 1); + 2). + + sip = s. We 
will be primarily interested in the symmetric group, for which 
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s! 
is oe aad 

7. For free Husimi trees all points belong to the same set. The dissimilarity 
characteristic relation (equation (6)) is however also valid in the more general case. 

%. Harary, F., and Norman, R. Z., ‘The Dissimilarity Characteristic of Husimi 
Trees,”’ Annals Math.; (in process of publication ). 
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ZMn, f(x)) — Z(San, f(x)), 
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We also make frequent use of Pélya’s result, p. 164, that for two independent figure 
collections, the corresponding counting series is the product of the individual ones. 


10. In this and all subsequent formulas for the ¥i(x), the extra factor x is needed to 
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STUDIES IN THE CONFORMAL MAPPING OF RIEMANN 
SURFACES, 1 
By Maurice Hets* 
BROWN UNIVERSITY AND UNIVERSITY OF PARIS 
Communicated by Marston Morse, February 20, 19538 


1. We are concerned with an analysis of the exact form of the Lindeléf 
inequality for directly conformal maps of one Riemann surface with posi- 
tive boundary into another. Let F and G denote Riemann surfaces with 
positive boundary, let Gp and G,, denote their respective Green’s functions, 
let f denote a directly conformal map of F into G (not necessarily univalent 
onto), and let m(p) denote the multiplicity of f at pe F. The exact form of 
the Lindelof inequality may be taken as the assertion that for such F, G, 
f for each q ¢G, the residual term 4, in 


Gel f(p), gq] = +? n(r)Ge(p, r) + u,(p) (1.1) 


fir) =q 


is non-negative harmonic on F. We observe that u, is the greatest har- 
monic minorant of Gel f(p),q]. 

We shall use the following terminology.” A non-negative harmonic 
function on a Riemann surface will be termed qguasi-bounded provided that 
it may be represented as the limit of a monotone non-decreasing sequence 
of non-negative bounded harmonic functions all having the given surface 
as their domain; a non-negative harmonic function on a Riemann surface 
will be termed singular provided that the only bounded non-negative 
harmonic function on the surface dominated by it is identically zero. Par- 
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reau’ has shown that a non-negative har:nonic function on a Riemann sur- 
face admits a unique representation as the sum of a quasi-bounded har- 
monic function and a singular harmonic function having the given surface 
as their domain. Our first result is: 

THEOREM A: (1) (p, gq) > 4,(p) ts upper semi-continuous on F X G. 
Let v, denote the quasi-bounded component of u,. Then (p, q) — v,(p) is 
lower semi-continuous on F XG. (2) Either v, = 0 for all q ¢ Gor else for 
nog. (3) The set of gq for which u, — v, > Ois an F, of capacity zero. 

Part (3) generalizes in several respects a corresponding result proved by 
Frostman! for Blaschke products. The proof of (1) depends upon the up- 
per semicontinuity of g —~ u,(p) for each p and the lower semicontinuity of 
q > v,(p) as well as uniformity properties of non-negative harmonic func- 
tions that are assured by the Harnack inequality. (2) depends upon 
known properties of the Green’s function. The proof of (3) lies deeper as 
might be expected. The key idea of the proof lies in showing that for a 
Green's potential P on G, where 


P(q) = S'Ge(q, s)du(s), (1.2) 


uw being a non-negative mass distribution on G, with compact kernel, the 
greatest harmonic minorant of P o f is 


S udlu(q) (1.3) 


and the quasi-bounded component of (1.3) is 
JS vdu(q) (1.4) 


2. Maps of Type-Bl. We term f of type-Bl provided that v, vanishes 
(for allg). Incase F = G = (2 < 1), this is tantamount to saying that 


for a <1, 


j— a 
| — af 


is a Blaschke product save for a set of a of capacity zero. 

We say that f is of type-Bl at q(eG) provided that there exists a relatively 
compact Jordan region Q, g «2 ¢ G such that (1) f~'(2) is not empty, (2) 
the restriction of f to a component of f~'(2) is of type-Bl relative to Q for 
each component of f -'(2). We say that f is /ocally of type- Bl provided 
that f is of type-Bl at each point of G. It is to be remarked that in these 
latter definitions the hypothesis that F and G have positive ideal boundaries 
may be dropped. We have: 

THEOREM B: (1) /f G has (as previously) positive ideal boundary and f 1s 
of type-Bl, then for each region  ¢ G, the restriction of f to a component of 
f '(Q) ts of type Bl relative to 2 for each component of f (2). (2) If Ghas 
positive ideal boundary and f ts locally of type-Bl, then f ts of type-Bl. 
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An essential component of the proof is the following lemma. 

LemMA C. Given u subharmonic on F and f locally of type-Bl, let U be 
defined on G as the upper function of sup u(p). If U< + @, then Us 

f(p) =a 
subharmonic on G. 

A consequence of theorem B is that the classical theorem of Iversen may 
be extended to maps which are locally of type-Bl. Lemma C admits a 
number of applications to other problems. By way of illustration we cite: 

THEOREM TD: Let G have positive ideal boundary and let u be singular 
positive harmonic on G. If f:F — G is of type-Bl, then u » f is singular on 
F. 

Theorem D leads to the following consequence for convergent infinite 
Blaschke products: the set of Fatou boundary values of modulus one of a 
convergent infinite Blaschke product is the unit circumference and each 
such Fatou boundary value of modulus one is attained at infinitely many 
points of the unit circumference. 


* The author is a Fulbright Research Fellow and President's Fellow of Brown Uni 
versity for the academic year 1952-1953. 

! Frostman, O., Potentiel d’équilibre et capacité des ensembles avec quelques applica 
tions a la théorie des fonctions. Lund Thesis (1935). 

2? Parreau, M., Sur les moyennes des fonctions harmoniques et analytiques et la classi- 
fication des surfaces de Riemann. Paris Thesis (1952). 


THE PARAMETRISATION AND ELEMENT OF VOLUME OF THE 
UNITARY SYMPLECTIC GROUP 
By F. D. MuRNAGHAN 
INSTITUTO TECNOLOGICO DE AERONAUTICA, SAO José pos Campos, BRASIL 
Communicated January 21, 1953 


We have given previously in these PROCEEDINGS! a system of parameters 
for, and the complete element of volume (including the in-class factor as 
well as the well-known class factor) of, the n-dimensional rotation and 
unitary groups and we complete, in the present note, this problem for the 
classical groups by giving a system of parameters for, and the complete 
element of volume of, the 2k-dimensional unitary symplectic group. 

The 2k-dimensional symplectic group, over the field of complex num- 
bers, may be presented as the collection of 2k X 2k matrices, XY, with com- 


: ae 0 —K, 
plex elements, which are such that X’/,X = /:, where 14 = E 0 ‘). 


_ 


Writing Y asa 2 X 2 block matrix ° ah whose elements are k X k ma- 
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trices, Y is symplectic if, and only if, B’A and D’C are symmetric and 
D'A — CB=Kk,. Thusa 2 X 2 matrix is symplectic if, and only if, it 
is unimodular. It is known that all symplectic matrices, for any value of 
/, are unimodular but the proofs we have seen of this fundamental fact are 


AA i ceieatinaeae tomate © Vere 
unnecessarily complicated. n tac 0 E, BD = BD so tha 


(det D’) (det Y) = det D. If, then, D is non-singular det Y = 1. If D 


E, O\. 
is symplectic and unimodular, no 


e e k ( 
is singular we observe that 
ek, E, 


; kK, 0 
matter what is the complex number ¢, and we form the product EE 
€ ftp [hy 


a& 4 . which is symplectic (being the product of 
= . Mie S Sy ccuc . y > oO y ) 
BD cA+BeC+D _— aed 

two symplectic matrices). Since eC + D is singular for at most k values 
A Cc 


i ees ~ is unimodular and this 
€/ € 


of ¢ there exists an ¢e for which ( 


BD 


The unitary symplectic group of dimension 2k is the collection of 2k X 


— c\. 
implies that ( ) is unimodular. 


2k matrices which are at once symplectic and unitary. Its typical matrix 


BA 
It is a k(2k + 1)-parameter group, & of the parameters being class param 


er : A «3B oe :, 
is of the form ( i) where B’A is symmetric and A*A + B*B = fy. 


eters and 2k* in-class parameters. Thus when k = 1, in which case the 
group is the unimodular unitary 2-dimensional group, we have | class and 
2 in-class parameters; when k = 2, we have 2 class and 8 in-class param- 
eters; when k = 3, we have 3 class and IS in-class parameters and so on. 
It is natural, then, that the problem of determining the in-class factor of the 
element of volume of the group is a more serious one than the already solved 
problem of determining the class factor. Any unitary symplectic nm.atrix 


; A-B : , 
X = ; } may be transformed, by means of a unitary symplectic 


B A 


matrix V, into diagonal form: V* X Vo = = 5). If the diagonal ele- 
ments of D are ee, ..., e’” the k angles am, ..., a, (each of which may be 
taken to lie in the interval 0 < a < mr) may be taken as the class parameters 
of the group and the class factor of the element of volume of the group is 
well known to be 

h 
aC sin” a, IL }1 COS( ay a,)} | 1 — cos(a, + a,){day... day 


/ 
I 


or the product of this by any positive constant. 
In order to obtain a convenient svstem of in-class parameters we 
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sider the particular unitary symplectic 2k X 2k matrices (pq), each of 
which involves two parameters 6 and 7, p and q > p being any pair of dis- 
tinct numbers from the set 1, ..., & or, else, p being any number of this 
set and q being p + k. The notation will be clear from the following ex- 
amples for the case k = 3: 

0 0 ? 0 —s 0 

cet —$§ 0 0 0 

S ce , l 0) 0 

0 0 , s Q ce tt 0 

0 0 0 0 0 

0 0 cer 0 0 1 


where ¢ = cos 8 and s = sin 6. We place a subscript under the paren- 
theses (pg) when we have to distinguish different pairs (8,7); thus we 
obtain (23) by replacing 8 and 7 by 6, 7; in the formula above, and so on. 
1 
Then any 2 X 2 unitary symplectic matrix X can be written as (12) 6 where 
6 is a diagonal unitary symplectic 2 & 2 matrix; any 4 X 4 unitary sym- 
plectic matrix may be written as (12) (24) (12) (13)6, where 6 is a diagonal 
. Se 1 
unitary symplectic 4 X 4 matrix; any 6 X 6 unitary symplectic matrix may 
be written (12) (23) (86) (12) (23) (12) (25) (12) (14)6, where 6 is a diagonal 
i) 8 7 4 3 2 l 
unitary symplectic 6 X 6 matrix and so on. The factors preceding 6 may 
be conveniently separated into sets of 1, 3, 5,..., (2k — 1) each (counting 
from the right); thus, when 2k = 4, we write {(12) (24) (12)} {(13) $6; 
4 3 2 1 
when 2k = 6 we write {(12) (23) (86) (12) (23)} {(12) (25) (12)} {(14) 46 
9 8 7 6 5 4 3 2 1 
and so on. Each of the angles 8 lies in the interval 0 < 6 < 2/2 while 
each of the angles + may lie in any one of the 4 quadrants. Turning to 


‘ DO 
the formula Y = | F D 


7 6 5 


) V* we suppose V factored in the manner described 


: ' DO\. a 
and observe that since 6 commutes with ’ >) it may be omitted. Thus 


"is the product of 1 + 3+ ...+ (2k — 1) = k* factors, each involving 2 
parameters 8, 7, and we take these 2k* parameters as the in-class param- 
eters of the group. 

The in-class factor of the element of volume of the group possesses the 
cumulative property (with respect to the dimension of the group) that the 
element of volume of the (2k + 2)-dimensional group contains as a factor 
the element of volume of the 2k-dimensional group. The coefficient of 
(dg) (dr), i. e., of the product of the differentials of the k? 8's and k? 7's, 
is a function of the 6’s only and always contains sin 26, sin 26)... sin 
28,2, which we shall denote simply by II(sin 28), as a factor. For k = 1 
the in-class factor of the element of volume is 


dV» = (sin 28)d6 dr 
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(it being always understood that the element of volu:rne may be multiplied 
by any convenient positive constant). For k = 2 it is 


dV, = (cos! Bs) (sin® By) IT (sin 26) (dB) (dr) 


there being 4 p's and 4 7's. For k = 3 it is 
dV. = (cos® 85) (cos Bs) (sim? By sin’ 6; sin® B,) (sin® Bs) IT (sin 28) (dB) (dr) 
there being 9 8's and 9 7's. Fork = 4 it is 


dV, = (cos* Bw) (cos 85) (cos? By) (sin? By5 sin’ By4 sin® By, sin? By, sin’ By) X 
(sin? B, sin' 8; sin® 35) (sin® By) Il (sin 28) (dB) (dr) 


there being 16 3's and 16 7's; and soon. IndI>, the part involving the 
cosines starts out with cos” 8), ),.,), each factor of this part being 
cos” By; pean J = Rk, Rk — 1,...,2. The part involving the even powers of 
the sines starts out with sin? 8,.—,, and consists of various factors corre- 
sponding to the described grouping of the factors of |’. For example, when 
k = 5, in which case we have 50 in-class parameters, the first (and largest) 
of these groups of factors is (12) (23) (84) (45) (5.10) (12) (23) (34) (45) and 
25 24 23 22 21 20 iv 18 17 
the corresponding part of the factor which involves the even powers of the 
sines 1s 
sin” Boy sin? Bo; sin® Bo» sin® Bo; sin® By sin’ Bis sin® By. 


' Murnaghan, F. D., these PROCEEDINGS, 38, 69-73, 127-129 (1952) 


HEAT CONDUCTION ON ARBITRARY RIEMANNIAN 
MANIFOLDS 
By ID. C. SPENCER 
DEPARTMENT OF MATHEMATICS, PRINCETON UNIVERSITY 
Communicated by S. Lefschetz, January 12, 1953 
Let JJ = A" be an arbitrary orientable Riemannian manifold of class 
C’ and of dimension n. We consider the exterior differential forms on A, 
and we use the notation of de Rham? in which the Laplace-Beltrami opera 


tor is the negative of the usual one. The scalar product of two p-forms 


¢, wis 


(9,0) = Su eAt¥ = (y, 9) 
while the Dirichlet scalar product is 


Dig, y) = (dg, dy) + (d¢, by) 
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where d is the operator of exterior differentiation and 6 (applied to p-forms) 
is equal to (—1)"" +" * '*d*. We define the norm (g) by the usual form- 
ula |g ? (y, ¢), and we write Dig) = Dig, ¢). Wesay that the p-form 
¢ isin the domain of the operator d if ¢ is of class C! and ¢  < @, 
dp || < ®. We define the domain of 6 in a similar way, and we say that 
¢ is the domain of A if ¢ 1s of class C? with |} ¢ || < ©, ||} Ag!) < ©. The 
closures of these domains are defined in the usual Hilbert space sense. 
We introduce the following spaces where the superscript p refers to the 
degree of the forms: 
A” 
Gg eC, ¢ with compact carrier}, 


{ o!, B’? = }y~ gin the closures of d, 6}. 
' 
' 
pb { yg in the closures of d, 6, A}, 
‘ 
' 
' 
' 


7) 
¥ 


G? ¢g¢D’, Dig, ¥) = (¢, Ap) for every y « D’}, 
NP ve D’, Dig, ¥) = (Ag, ) for every y € D’}. 


If B! is the subspace of B’ formed by taking the closure, in the norm 
Dig) + | ¢ *, of the space C’, it turns out that G? ¢ BP. 
Let s be a positive number, and write 


Dog, ¥) = De, v) + sle, W), Ag = A+ s. 


The Green's operator G, is the one-one linear mapping of A’ onto G’ 
whose inverse is A,, while the Neumann's operator N, is the one-one linear 
mapping of A” onto NV’ whose inverse is Ay. 

THEOREM. For each positive s the Green's and Neumann's operators 
exist. They are symmetric operators and they satisfy D&G, ¢) < ‘¢ll?/s, 
DAN, ¢) <4 ¥¢ mT Ss; 

Let 

F’(\) = le ¢ «€ A”, (N, — A)e = O} 
be the space of eigenforms belonging to the eigenvalue A. Since N, is 
bounded, therefore continuous, F?(A) is a closed space and, moreover, we 
have 


FRA) = tel ee N’, Ase = ¢/A}. 


In fact, if Pe N’, AW = 0, then 0 = (A,y, Y) = D,(W) so W vanishes identi 
cally ifs > 0. If ge F’(d), then) ¢ ?/A = (Avg, ¢) = Dy(¢) > 5\\¢/)?, and 
it follows that A > 1/s with equality if, and only if, Dig) = 0. Hence 


Fe = F?(s) = ty|ge A¥, dp = 5g = O}. 


Since the Neumann's operator N, is bounded and symmetric, therefore 
self-adjoint, we have by the spectral theorem 


1/s 
N, = z d dE,(A) 
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where /,(A) is the canonical spectral field associated w‘th N,. 
Finally, let ¢ > 0 and define 


W, = JS’ expl(as — 1)t/A] dE,(A) 


Mor each fixed ¢ > 0, W, is a symmetric operator which is bounded by 1. 
Moreover, since W,:-W,, = W,,,,, these operators form a semigroup. 
We see that W, is the identity operator, W., = lim W,(f—  ) is the opera- 
tor F = F’ which projects orthogonally onto the space F’, and that (A + 
0 O01)W, = 0. Given ae A’, Wa is a non-increasing function of ¢ and 
W,a is a solution of the heat equation with initial distribution a. 

Analogous considerations apply to the Green's operator, but the space 
F” is replaced by the subspace 


FP = tel ge BY, de = b¢ = O}. 


The heat operator in this case, call it V,, differs in general from W, in that 
it tends, as / > ©, to the operator F, = F? which projects orthogonally 
onto the space F?. 

We say that J/ is a totally bounded manifold (shortly, a 6-manifold) if 
M7 is a compact subdomain of some larger Riemanman manifold R. If A/ 
is a b-manifold, there exists a Green's operator Gy which is a completely 
continuous linear mapping of A’ onto G’? — F! whose kernel is F? and 
whose inverse mapping is A. Moreover, in this case, F! is the subspace of 
harmonic p-forms on R which vanish identically outside J/. 

The opposite of a 6-manifold is a w-manifold (totally unbounded mani 
fold) which is characterized by the property that, for each p, 0 < p <n, 
every p-form in the closure of the operator A belongs to N’. If A/ is a 
u-manifold, the spaces G? and N’ coincide, so G, N, for s > Oand W, 
V,. Let 


H? = {¢ A’, Ag = 0}. 


On a u-manifold the space //?”, s > 0, contains only the form 0, and H? = 
Il, = F® = F?. Furthermore, the operators d, 6 and A commute with 
N, and W,. It follows, in particular, that the period 


P(a, C”) Sc Wia 


of any given closed form @ ¢« A” on a finite (differentiable) p-cycle is inde 
pendent of ¢, and therefore 


P( a, C*) Se Fa = P(Fa, C’). 


Finally, we observe that, on a u-manifold, the closure under the norm of 
the space of forms dg, ¢ « A” ', is the same as the closure of the space of 
I 


forms dg, ¢g « C’~', and similarly for the operators 6, A. A u-manifold 


therefore has negligible boundary in the sense of Gaffney,' that is (da, 3) 
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(a, 63) for every ain the domain of d and every 8 in the domain of 6. In 
fact, a u-manifold is one which has negligible boundary in Gaffney’s sense 
and satisfies the further condition that any form ¢ in the domain of A satis- 
fies dg < , 6g < », éddg < , débg) < © and hence Ag = dég + 
é6dyg. Compact (closed) manifolds are a subclass of u-manifolds; the heat 
equation on compact manifolds was considered by Milgram and Rosen- 
bloom.# 

There are manifolds which have negligible boundary but are not u- 
manifolds. For example, let R be a compact Riemannian surface and let 
¢ be the potential function on R which has a (logarithmic) source at one 
point, a sink at another point, and is regular harmonic everywhere else. — If 
we remove from R the source and the sink points, we obtain a manifold A 
which has negligible boundary, but 7 is not a u-manifold since ¢ lies in 
the domain of A(Ag = 0) but not in the domain of d. 

' Gaffney, M. P., “The Harmonic Operator for Exterior Differential Forms,’’ Proc 
Nati. Acapb. Scr., 37, 48-50 (1951) 

*de Rham, G., and Kodaira, K., “Harmonic Integrals,"’ mimeographed notes, 
Institute for Advanced Study, Princeton, N. J. (1950). 

* Milgram, A. N., and Rosenbloom, P. C., ‘“Heat Conduction on Riemannian Mani- 
folds I, II,’ Proc. Nati. Acap. Scr., 37, 180-184 (1951), and 37, 435 488 (1951) 


ABSOLUTE VALUED ALGEBRAS* 
By FrepD B. WRIGHT 
UNIVERSITY OF CHICAGO 
Communicated by A. A. Albert, February 26, 1953 


Let A be a not necessarily associative algebra over the real field R, which 
is a normed linear space under a norm satisfying, in addition to the usual 
requirements, the equation |\xy) = x): y),forevery x, yof A. Such an 
algebra is called absolute-valued. A. A. Albert has shown! that a finite 
dimensional absolute-valued algebra with a unit element is isomorphic to 
one of the following: (1) the real field R; (2) the complex field C; (3) the 
quaternion algebra Q; (4) the Cayley-Dickson algebra D. He has further 
shown® that an absoiute-valued algebraic algebra is finite dimensional; 
hence the preceding result applies. Taken with the Gelfand-Mazur theorem 
this latter result implies that the same theorem holds for power-associative 
absolute-valued algebras with unit element. This suggests the possibility 
of extending the result to arbitrary absolute-valued algebras. The theorem 
of $2 shows that the theorem indeed holds for absolute-valued division alge- 
bras. 

$7. Preliminary Considerations. For the reader's convenience, the 
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relevant definitions and results of the above cited papers will be set forth 
in this section, but no proofs are given here. 

An algebra A over a field R is a vector space over R closed with respect 
to an operation called multiplication, denoted by xy, which satisfies iden- 
tically the following for all x, y, s of A and a of R: (x + y)z = xz + ys; 
a(x + y) = sv + sy; (ax)yv = a(xy) = x(ay). The algebra is said to have 
a unit element e if there is an element e of A such that ex = xe = wx for 
every x of A. The unit element will always be denoted by the same symbol 
as the number |. We may then identify R with the subalgebra of A gener 
ated by 1. 

The linear transformations R, and L, on the vector space of A are defined 
for each a of A by x ~ xa = xR, and x ~ ax = xL,. The algebra A 1s 
called a division algebra 1f both R, and 1, have inverses as linear transfor- 
mations for every a # Oin A. 

Throughout this paper, RX is taken to be the real field. Then an algebra 
over R is called absolute-valued if there is a mapping x ~ |x of A into R 
satisfying: |x| > 0; |x| = Oif and only if x = 0; |xy) = |x/-lyl; 

x+y < ix) + ly); lav) = |a|-!x| for every x, y of A anda of R. 

Two algebras A and B over R are isotopes of each other if they have the 
same underlying vector space and if multiplication vey in A is defined in 
terms of the multiplication xy in B by the relation voy = (xP)(yQ), where 
P, Q are fixed invertible linear transformations on the common underlying 
vector space. 

Every element x of A generates a subalgebra Al{x| of A. A is called 
algebraic if A[x] is finite dimensional for every x of A. 

The principal results of Albert are stated here as lemmas 1-3. 

LEMMA |: An absolute valued division algebra is an isotope of an absolute 
valued algebra with a unit. 

LEMMA 2: An absointe valued algebraic algebra with a unit element 1s 
finite dimensional. 

LeMMA 3: A finite dimensional absolute valued algebra with a unit element 
ts isomorphic to one of the following: (1) the real field R; (2) the complex 
field C; (3) the quaternion algebra Q; (4) the Cayley-Dickson algebra D. 

$2. Absolute Valued Division Algebras. 

LemMMA 4: An absolute valued division algebra with a unit is an inner 
product space. 

Prooj: tis sufficient! to show that for every x, yof A with x = y 
1 the inequality x + y°> + x — y\* > 4 holds. In the identity (1 4 
x)? — (1 x)? = 4x, let |x} = 1; since |1/ = 1, it follows that 


lItwi?+il—wte4 (1) 


r *¢ - . ] 
Now if a, 6 are any two elements of A with a = 6) = I, letx = aR, ; 


’ 


then x = 1. Furthermore, (1 + «)R, = / + a and (1 WR, = b—a. 
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It is clear that 1 : b+ ai and! 1 x| = - . Substitution 


in (1) gives 
b+al?+ib—al?>4. (2) 


THEOREM: = /tvery absolute valued division algebra with a unit ts algebraic. 
Proof: ‘There is now an inner product (x, y) defined on A with (x, x) = 
x *. By standard procedures this inner product space may be completed 
to be a real Hilbert space //, in which A is contained as a dense linear sub- 
space, and in which the inner product on // agrees with the inner product on 
A, for elements of A. Let / be the algebra generated by 1 in A; M isa 
one-dimensional subspace of // and is hence closed. We write // = M @ 
AM’, where AM’ is the orthogonal complement of Min //. Every element x 
of A has a unique expression as x = a + by, where a, ) are in R and y is an 
element of A contained in AJ’, yy = 1. Since (1, y) = Oand (y, y) = 1, 
it results that (1 + y, 1 + y) = 2and (1 — y,1 — y) = 2. Butifu,v are 
any two elements of A, (u, u)(v,v) = (u)*iv ? = |uv)* = (uv, uv). Then 
it results that (1 — y’, | — y®) = 4, or (1, y?) = —1. But y? can be writ- 
ten y? = ¢ + dz, where c,d are in Rand (zg = 1, zisin both A and M’. 
Thenc = —1,d = 0; ory? = —I. 

For arbitrary x of A, with x = a + by as above, define x* = a — by; 
then x + x* = 2a,x*x = a? — b*. Butx®? — (vw + x*)x + x*x = 0 is an 
identity, so that x? — Bax + (a* — b*) 0, which shows that A is algebraic. 

Then the use of lemmas 1-3 above gives a complete description of all 
absolute valued division algebras, even without unit element. 

$3 Remarks. In the proof of the theorem, use is made of the assump- 
tion that A is a division algebra only in passing from (1) to (2). It is 
natural to suppose that (2) can be proved independently of (1). The 
ptolemaic inequality’ seems made to order for this. 

The connection of this result with the Gelfand-Mazur* theorem seems 
worth noting. This result says that an associative normed (the weaker 
condition xy) < ox + y is assumed) algebra over the reals, without 


topological divisors of zero, is either the reals, complexes, or quaternions. 
The difficulty in extending this result to algebras which are not assumed to 
be associative is the lack of an adequate definition in such algebras of the 


spectrum. It is also to be noted that there exist normed division algebras 
over the reals which are not absolute valued. All known examples are 
finite dimensional; it may be conjectured that every normed division al 
gebra over the reals is finite dimensional. 

* This paper was prepared under contract with the Office of Naval Research 

' Albert A. A. Annals of Math. 48, 495-501 (1947) 

2 Albert, A. A., Bull. Amer. Math. Soc., 55, 763-768 (1949) 

’ Kaplansky, I., Duke Math. J., 16, 899-418 (1949) 

' Schoenberg, I. J., Proc. Amer. Math. Soc., 3, 961-964 (1952) 
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AN ANALYSIS OF Til SYNTHESIS AND DISTRIBUTION OF 
THE CONTRACTILE PROTEIN, MYOSIN, IN THE 
DEVELOPMENT OF THE HEART 


By James D. ERBER1 
INDIANA UNIVERSITY * 
Communicated by B. H. Willier, February 10, 1958 


The development of the heart presents a series of problems in morpho 
genesis and differentiation which afford an unexcelled opportunity for the 
analysis of the interrelationships of specific structural and functional prop 
erties during ontogeny. As the heart is formed, so it begins to beat. In 
the chick embryo of 7 somites the paired primordia of the heart have already 
begun to fuse; local contractions have been recorded in the heart of the 
9-somite embryo." * The onset of contractility has been correlated with 
the time at which significant histogenic changes have been recorded, 
e.g., the first appearance of cross-striations at the 10-somite stage*® and the 
first detectable glycogen in the embryo of “about 2S hours.’ The early 
contractions have been described as of a fundamental myogenic nature 
since they occur well in advance of the time at which the first nerve fibers 
approach the heart. 

The early morphogenesis and onset of contractility of the heart are 
foreshadowed by a precocious restriction of the heart-forming cells to rel 
atively well-defined areas of the blastoderm. Willier and Rawles’ and 
Rawles® have described two large bilaterally situated areas of the head 
process blastoderm endowed with heart-forming poteney. No critical 
information has been presented concerning the location of the heart 
forming areas earlier than the head-fold stage. Rudnick’ has concluded 
from a review of the pertinent data that the heart mesoderm must be verv 
early in a highly organized state, and is early invaginated to migrate lat 
erally to its definitive position. A similar conclusion was reached by Ba 
con® in his analysis of the development of the heart in A mblystoma. 

In contrast to the fund of information from ‘the field of experimental 
morphogenesis, information concerning the chemistry and cellular physt- 
ology of the development of the heart is almost completely lacking. The 
void may be due, at least in part, to the difliculties encountered in studying 
the contractile proteins in minute quantities, as illustrated by the expe 
riences of Herrmann and his coworkers,®* '’ who have not been able to make 
determinations of actomyosin in skeletal muscle of the chick embryo prior 
to the ninth day of incubation, although fibrillations and cross-striations 
hive been observed in the 7-day embryo, a time at which muscular re 


sponse to mechanical stimulation and histological differentiation of motor 


centers in the central nervous system also have been described. In the 





304 ZOOLOGY: J. D. EBERT Proc. N. ALS 


present investigation the difficulty in detecting minute quantities of cardiac 
contractile proteins is largely overcome through the application of the pre- 


cise techniques of immunochemistry. 

In earlier reports'' '* it was shown, by absorption and precipitin 
tests, that antisera produced in rabbits against extracts of adult chicken 
heart will react with extracts of the early chick embryo prior to the mor 
phogenesis of the heart. The early blastoderm must therefore contain 
antigens identical, or closely related, to those of adult heart. This con- 
clusion is supported by the finding that when early blastoderms are cul- 
tured in media containing antiheart sera, the development of the heart is 
differentially suppressed. 

It becomes of interest to determine next the relationship of the specific 
antigenic components of cardiac muscle present in the early blastoderm to 
the contractile proteins, actin and myosin. We ask the following questions: 
(1) Do actin and myosin develop simultaneously? Is the development of 
each protein an independent event or does the development of one precede 
that of the second and have a causal role in its synthesis ? 

(2) Are there differences in the contractile proteins of cardiac and 
skeletal muscle which play a significant role in the precocious restriction of 
cardiac contractile proteins to the heart-forming areas of the blastoderm, 
and in the early onset of contractility? Only after the facts relative to the 
time of synthesis and distribution of the cardiac contractile proteins are 
available will it be possible to attack the basic problem of the factors con 
trolling their synthesis. The specific objective of this report is to present 
findings relative to the time of synthesis and progressive localization of 
cardiac myosin in the early chick embryo. 


METHODS AND RESULTS 


1. Preparation of Myosin. Cardiac and skeletal myosins were extracted 
from heart and hind limb muscle, respectively, of adult New Hampshire 
Red chickens of both sexes, following the method described by Mommaerts 
and Parrish,'’ in which myosin and actomyosin are simultaneously pre 
cipitated from the muscle extract by dialysis, most of the actomyosin being 
irreversibly insoluble. The remaining actomyosin is precipitated from the 
redissolved myosin by adjusting the ionic strength and pH. <Actin-free 
myosin is then crystallized by dilution with glass-distilled water. Twelve 
to 36 fresh adult chicken hearts were used in each extraction (100 to 350 
g. of fresh cardiac muscle). Two minor deviations from the technique 
described by Momunaerts and Parrish should be noted: (1) following 
dialysis for 1S hours, the ionic concentration of the muscle extracts 
is frequently not reduced sufficiently to permit maximal precipitation; 
therefore, 200 ml. of glass-distilled water per 100 g. of muscle is added 
directly to the extract following dialysis. (2) Mommaerts and Parrish 
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describe an average yield of 1.2 g. of myosin from 100 g. of muscle. Our 
vields have been consistently lower: from chicken cardiac muscle, 750 mg. 
of myosin per 100 g., and from chicken skeletal muscle, 410 mg. of myosin 
per 100 g. of muscle extracted. 

It must be emphasized at the outset of the discussion that although the 
extraction technique employed yields actin-free myosin which approaches 
purity, the final product is not monodisperse. Mommmaerts and Parrish 
state that the purity of their preparations is of the order of So-90 per cent. 
Preliminary studies of cardiac myosin indicate that the degree of purity of 
the preparation is only slightly lower than that described by Mommaerts 
and Parrish. 

2. Production of Antisera. Despite the brief positive report by Kesz 
tyus, Nikodemusz,and Szilagyi'' that myosin 1s antigenic, our early attempts 
to produce potent antisera in rabbits against cardiac myosin from adult 
chickens were unrewarding.'’ Hence it seems advisable to record here, 
briefly, the more important details of the immunization procedure which has 
proved to be reliable in our laboratory during the past year. First, it ts 
essential that the myosin used for injections be either freshly prepared or 
stored for less than one week at —20°C. In the latter case, it should be 
thawed only once. The use of freshly extracted myosin for each injection 
is prohibitively expensive of both time and materials, hence we have made 
extractions weekly. The myosin should be a bluish-translucent gel; no 
shreds are desired. The myosin is stored and in’ected in the phosphate 
KCI solution described by Mommaerts and Parrish, at pH 6.8, tonic 
strength 0.5. 

Individual rabbits vary greatly in their antibody response to injected 
myosin. To date, of 22 rabbits injected with cardiac myosin, only 12 
have produced antisera of sufficient potency to be of value. In general, the 
immunization schedules are so arranged that each adult male rabbit re 
ceives a series of intravenous injections on alternate days over a two-week 


period, totaling 18-20 mg. of protein, is rested for 4-6 days and bled. 


After 20-25 days further rest, a second series of injections is started, in 
which a first injection of 4-6 mg. of protein is given intraperitoneally, fol 
lowed by a total of 48-60 mg., intravenously, on alternate days over a two- 
week period Once again the rabbits are bled following 4-6 days rest. 
Our rabbits have consistently built up circulating antibodies against myosin 
rapidly, and lost them rapidly. Hence the rabbits are bled 4-6 days rather 
than the more general 7-10 days after the last injection of a series. Sera 
are sterilized by passage through a Seitz filter and stored at —20°C. No 
preservative is added. 

3. Specificity of Antisera. The titer and specificity of each serum was 
determined by the precipitin technique, employing homologous or heter 
ologous antigen. A preliminary titration was made with a constant 
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volume of full-strength serum and progressively decreasing quantities of 
homologous antigen. The end-point of the antigen dilution series was then 
titrated with serial dilutions of serum. In this manner, inhibition due to 
excess antigen is avoided, the serum end-point being greatest with the last 
active dilution of antigen.” ~~ A total of 54 antisera from 22 rabbits has 
been assay d. Appropriate controls have been carried out with sera of unin- 
jected rabbits. The effective serum dilution titers of the antisera tested 


ranged from 1:64 to 1:S192. Only those sera with titers of at least 1:512 


with homologous antigen were retained (31 antisera). Potent sera were 
then tested for reactivity with heterologous antigen (skeletal myosin). 
The results clearly show that the antisera against cardiac myosin are highly 
tissue-specific, 1.e., that cardiac and skeletal myosins differ markedly in 
antigenic constitution. The highest titer obtained in tests of 27 antisera 
against cardiac myosin with nine separate fresh preparations of the heter- 
ologous antigen, skeletal myosin, is 1:256. There is a consistent difference 
of 4-5 tubes in the reactivity of each serum with homologous and heterol- 
ogous antigen. In further tests, antisera were first completely absorbed 
with adult skeletal myosin, and then tested with homologous antigen. 
Much of the specific reactivity of the sera was retained. The reduction 
in reactivity due to elimination of antibodies to reactive groups in the 
heterologous antigen ranged from | to 4 tubes. Therefore, due to the 
marked difference in antigenic constitution of cardiac and skeletal myosins, 
it is possible to prepare, by absorption procedures, antisera highly specific 
for adult cardiac myosin. In the sections to follow, all references to anti- 
sera against adult cardiac myosin refer to these highly specific, absorbed 
antisera. 

#. Alyosin in the Developing Ileart. \t will be of interest to examine 
next the evidence pertaining to the antigenic specificity of cardiac myosin 
during development, as determined by an analysis of the reactivity of 
antisera against adult cardiac myosin with cardiac extracts from embryos of 
36 hours to 1S days incubation. All embryos used in these experiments 
were of the New Hampshire Red breed, from eggs incubated at 87.5°C. 
After removing the embryos from the eggs as quickly as possible, the hearts 
were removed (in an iced Petri dish), freed from excess blood with filter 
paper, weighed rapidly on a micro-torsion balance and minced finely with 
scissors on chilled waxed paper. The mince was homogenized in a Ten 
Brock glass tissue grinder in buffered KCI solution. At IS days, 60 hearts, 
and at 13 days, 144 hearts were extracted by the technique described for 
adult cardiac muscle, with a workable yield of purified myosin. At the 
earlier stages, however, it was necessary to employ a crude extract of car- 
diac muscle (containing both myosin and actomyosin), since a measurable 
vield of purified myosin could not be obtained from the small quantities of 
material available. Hearts of embryos of from 36 hours to 12 days in 
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cubation were homogenized in the buffered KCI solution. The homogenate 
was clarified by filtration through cheesecloth followed by centrifugation 
at 3600 r. p. m. for 80 minutes at 1°C. Both the purified cardiac myosin 
from older embryos and the clarified cardiac extract from younger embryos 
were tested for reactivity with antisera against adult cardiac myosin em- 
ploying several modifications of the precipitin technique, including (1) the 
ring test, (2) the Hanks technique,’ ' and (3) the micro (capillary) 
method of Ten Cate and Van Doorenmaalen,'* kindly demonstrated to me 
by the authors in their laboratory. All preparations from cardiac muscle, 
both the purified myosin and the clarified extract from younger embryos, 
reacted positively with all antisera tested. Cardiac extracts and purified 
cardiac myosin from embryos at the following stages have been tested 
36, 48, and 72 hours, 5, 7, 9, 12, 13, 15, and IS days. At least three sep- 
arate fresh extracts have been tested at each stage, each extract being re- 
acted with at least three potent anticardiac myosin sera. Further, 
by the absorption technique, employing repeated absorptions with em 


TABLE | 
CARDIAC Myosin IN THE EARLY EMBRYO 
HAMBURGER TOTAL NUMBER OP NUMBER OF TOTAL NUMBER NUMBER 
AND NUMBER OF FRESH ANTISERA OF OF 
HAMILTON EMBRYOS EXTRACTS USED WITH TEST POSITIVE HIGHEST 
STAGE EXTRACTED rESTED FACH EXTRACT SERIES REACTIONS TITER 
l 264 5 F 15 0 
248 5 : IS l 1 
236 5 ¢ 15 15 26 
381 7 : d 17 17 256 
140 3 : d Ss S 912 


bryonic cardiac antigens of any of the stages tested, it has been clearly 
demonstrated that all activity of the antisera for adult cardiac myosin 
can be removed. It is concluded, therefore, that the embryonic heart, 
from 36 hours through 1S days incubation, contains antigens or reactive 
groups identical, or closely related to those found in adult cardiac miyosin. 
5. Cardiac Myosin in the Embryo Prior to Morphogenesis of the Heart. 


The question next arises as to the time at which the reactive groups of 


cardiac myosin first can be detected in the early embryo. In this series 
of experiments, embryos in stages | to 5 of the Hamburger-Hamilton series 
(pre-streak through head-process stages) were extracted following the 
method deseribed for the extraction of cardiac tissue of embryos from 36 
hours to 12 days incubation. Three dozen or more embryos at each stage 
were rapidly removed from the yolk, trimmed and washed to remove 
adhering yolk, and homogenized in ice-cold KCl-phosphate solution in 
the Ten Brock grinder. The homogenate was allowed to stand for | hour 
at 2°C., with constant stirring, and then clarified by filtration through 
cheesecloth, followed by centrifugation. The resultant clear extract was 
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used as the test antigen in precipitin tests with antisera against adult 
cardiac myosin, chiefly by the Hanks’ and capillary methods. The re- 
sults, which are summarized in table 1, are based on 73 separate serum- 
dilution series of reactions, for which a total of 1269 embryos has been ex- 
tracted. It is clear that reactive groups identical with or closely related 
to cardiac myosin first can be detected in the embryo by the methods em- 
ployed at stage 3, the mid-streak blastoderm. All tests of pre-streak, and 
17 out of 1S tests of initial streak embryos were negative. At the latter 
stage one exception was noted—a trace reaction in the first 2 tubes of | 
series. The experiment was repeated with 3 additional fresh extracts of 
embryos at the same stage; all were completely negative with the same 
antiserum; therefore the trace reaction is not regarded as significant, 
likely being due to an error in staging one or more embryos in the extract in 
question. 

6. Distribution of Cardiac Myosin in the Early Embryo. ~Rawles* has 


PABLE 2 


ANTERO-POSTERIOR DISTRIBUTION OF CARDIAC Myosin 


NUMBER 
HAMBURGER TOTAL OF NUMBER OF TOTAL SECTION NUMBER 
AND NUMBER OF FRESH ANTISERA NUMBER OF OF HIGH 
HAMILTON EMBRYOS EXTRACTS USK WITH OF TEST EMBRYO POSITIVE EST 
STAGE EXTRACTED TESTED KACH EXTRACT SERIES TESTED REACTIONS TITER 


4 601 11 2-4 29 Ant 29 64 
Mid 29 256) 
(node } 
Post 29 256 
Ant 0 ae 
Mid $i 256 
(node) 


Post 


mapped rather precisely the shape, average dimensions and position of the 
heart-forming areas in the chick blastoderm at the head process stage. 
Her cogent analysis has provided the experimental target for the present 


investigation. Attention was directed first to a study of the distribution 
of cardiac myosin in the antero-posterior axis of the embryo at stages 4 
and 5. The blastoderm at stage 4 or 5 was divided into 3 segments, 
designated anterior, mid (containing the node), and posterior, by 2 trans- 
verse cuts made with glass or freshly sharpened steel needles across the 
blastoderm at the levels defined by, Rawles as the anterior and posterior 
limits of the heart-forming areas, 0.6 mm. in front of, and 0.45 mm. behind, 
the node, respectively. Anterior, mid, and posterior segments, respectively, 
were pooled in cold buffered KCI solution and extracted following the 
procedure described for whole embryos. The clarified extracts were tested 
for the presence of reactive groups of cardiac myosin by the capillary 
method. The results are summarized in table 2. It is clear that in the 
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embryo at stage 4, cardiac myosin is found in all 3 sections of the embryo 
in the antero-posterior axis. The data further suggest, but do not prove, 
that the amount of cardiac myosin in the mid and posterior segments is 
greater than in the anteriormost region of the stage 4embryo. In order to 
establish the observed difference as a significant one it will be necessary to 
so modify the present micro-techniques that strictly quantitative data can 
be obtained. Efforts in this direction are being continued. In the embryo 
at stage 5, cardiac myosin is confined to the mid segment in the antero- 
posterior axis, being absent from the anterior and posterior segments. 
Cardiac myosin is limited in its distribution in the antero-posterior axis, 
therefore, to a region slightly over one millimeter in length. The possible 
mechanism of the progressive localization of cardiac myosin will be dis- 


cussed in a later section, but it should be emphasized at this point that the 


rABLE 3 
MeEpbI0-LATERAL DISTRIBUTION OF CARDIAC Myosin 
HAMBURGER TOTAL NUMBER NUMBER OF TOTAL SECTION NUMBER 
AND NUMBER OF OP FRESH ANTISERA NUMBER OF OF HIGH 
HAMILTON EMBRYOS EXTRACTS USED WITH or TEST EMBRYO POSITIVE KST 
STAGE EXTRACTED TESTED RACH EXTRACT SERIES TRSTED REACTIONS TITER 
4 271 6 2 or 3 13 Left lat : 64 
Median 3 4 
(streak } 
Right lat 3 4 
Left lat. 
Median 
(streak ) 
Right lat 
Left lat 
Median 
(streak ) 
Right lat 11 128 


demonstrated localization is of a specific contractile protein and not merely 
of prospective cardiac celis destined for the synthesis of that protein. 

Next let us consider the distribution of cardiac myosin in the medio- 
lateral axis. Again the data of Rawles provide the morphological basis for 
our experiments. Groups of embryos at stages 4 and 5 were transected as 
previously described, 0.6 mm. before, and 0.45 mm. behind, the node. The 
mid-piece, however, was further subdivided into 4 sections, by 2 cuts per- 
pendicular to its cut surfaces, 0.2 mm. on either side of (and roughly paral'el 
to) the node. The resultant sections, designated left lateral, median, and 
right lateral, respectively, were pooled and extracted, and tested for the 
presence of the reactive groups of cardiac myosin as previously described. 
The data are are summarized in table 3. 

In the embryo at stage 4, cardiac myosin is approximately equally 
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distributed in the 4 sections tested in the medio-lateral axis, and has not 
been localized to the heart-forming areas. In the embryo at stage 5, al- 
though cardiac myosin is largely localized in the lateral areas, traces of the 
reactive groups are found in the median section (0.4 mm. wide) containing 
the streak. Thus the data are not in complete agreement with the finding 
by Rawles that the median area does not produce heart when grafted to 
the chorio-allantoic membrane. Although the reaction ‘s weak, it 1s con- 
sistently found and is believed to be significant. It appeared likely that 
cardiac cells are being invaginated even as the streak regressed. In order 
to determine whether the localization of cardiac myosin to the prospective 
heart-forming areas is complete at a later stage, the experiments were ex- 
tended to include stages 6-7 (stage 6 is the transitional head-fold stage, 
difficult to “‘schedule’; hence stages 6 and 7 are pooled). The results, 
shown in table 3, demonstrate that in the embryo at stages 6-7 the median 
(streak) region does not contain cardiac myosin detectable by the immuno- 
chemical methods employed. 


DISCUSSION 


The labile protein which has been characterized from the physicochemical 
point of view (at least within reasonable limits), and designated “‘myosin™’ 
by most workers in the field of muscle chemistry has been the myosin ex 
tracted from rabbit skeletal muscle. The further questions as to whether 
myosin 1s species- and organ-specific have received little attention. A 
recent study” of skeletal myosins prepared from tissues of the rabbit, rat, 
and mouse has led to the conclusion that the myosins from these species 
are practically indistinguishable on the basis of a comparison of electro 
phoretic properties, sedimentation rates, yields, viscosities, and specific 
refractive increments. On the other hand, the very fact that myosins have 
been shown to be antigenic argues for species differences of sufficient mag- 
nitude to account for their observed ability to stimulate antibody produc- 
tion. Kesztyus, Nikodemusz, and Szilagyi!’ report that rabbit skeletal 
myosin is antigenic in dogs, but that it is not isoantigenic. Our own findings 
show that chicken cardiac myosin is antigenic in the rabbit. 

With respect to the problem of the organ-specificity of myosin, Snellman 
and Gelotte®! have reported that cardiac and skeletal myosins of the rabbit 
show several differences, e.g., in solubility and in the sedimentation 
constant, and in the firmness with which actin and myosin are bound 
(less in cardiac than in skeletal actomyosin). The results of the present 


study clearly show that cardiac and skeletal myosins are markedly different 


in antigenic constitution, a difference which has been of importance in es- 
tablishing the specificity of the antisera employed in detecting cardiac 
myosin in the early embryo. The question remains unanswered, however, 
as to whether antigenic differences reflect important physical differences in 





Vow. 39, 19538 ZOOLOGY: J. D. EBERT 341 


the proteins, and is currently the objective of a comparative physical study 
in this laboratory. 

Next let us consider the question as to whether embryonic cardiac myosin 
is qualitatively the same as that in the adult. On the basis of differences 
in behavior of actomyosins from newborn and adult animals (guinea pigs, 
rats, mice) upon the addition of adenosinetriphosphate, Ivanov and Kasin 
ova*’ concluded that actomyosin in the newborn animal differs qualitatively 


10 


from that in the adult. Csapo and Herrmann" do not support this con 
clusion, arguing that, at least in their comparison of chick embryonic and 
adult actomyosins, only quantitative changes occur. Our findings support 
the latter view, inasmuch as the reactive groups of cardiac myosin capable 
of combining with specific antisera to adult cardiae myosin do not change 
qualitatively from the embryo of 36 hours to the embryo of IS days or, 
in fact, to the adult fowl. 

It will be of interest to consider next the questions of the time of syn 
thesis and localization of cardiac myosin in the early embryo, with partic 
ular reference to the concept of organ-forming areas. It is necessary to 
recognize at the outset of the discussion that while the data clearly demon 
strate that a major synthesis of cardiac myosin is initiated only after the 
prospective mesodermal cells of the epiblast have invaginated, the evidence 
does not permit a discussion of the problem of mechanism of synthesis. 
It is clearly possible that ‘‘traces”’ of cardiac myosin may be present in the 
earlier embryo-quantities which cannot be detected even by the highly 
sensitive methods employed. The present findings do present a sound 
chemical basis for further experiments designed to test the related hypothe 
ses of the regulation of the growth and differentiation of an organ by its 


own products’ °* and of pattern formation involving the specific suppres 
sion of embryonic differentiation by metabolites of differentiated tissues. *4 
We ask, for example, whether the reported failure of cardiac tissues to 
differentiate in the presence of extracts or products of fully differentiated 
heart muscle** ** results from a suppression of the synthetic activities of 
the tissue (e.g., in the synthesis of contractile proteims), or from a failure in 
morphogenesis of synthesized products. 

Tests of developmental potencies of small pieces of the chick blastoderm 
through transplantation to the chorio-allantoic membrane clearly show 
that each of the isolated pieces possesses a relatively restricted develop- 
mental capacity. Particular structures develop only from rather definite 
areas of the blastoderm; in short, the organ-forming areas are localized 
to a remarkable extent. Recently, the nature and significance of the con 
cept of the organ-forming areas has been the object of sharp eriticism.”” 
Yet most of the investigators who have employed this approach have been 
acutely aware of its major limitation, namely that each of the isolated 


pieces is triploblastic, and have, in fact, emphasized the importance of 
e 
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tissue interactions. Moreover, it has been recognized that the degree of 
specificity of localization of the several organ-forming areas differs. In 
fact, on the basis of his recent investigation of the localization of prospective 
neural plate in the early chick blastoderm, in which the results obtained 
by transplantation and explantation are shown to be in close agreement with 
those obtained by marking methods, Spratt’ has concluded that both 
types of methods reveal primarily the prospective fates of different regions 
f the blastoderm. In the case of the heart-forming areas, however, as 
Rudnick has emphasized,’ the areas are amazingly large, exceeding in 
size the actual myocardial anlagen. Thus it appears that the heart-forming 
areas are areas of prospective potency and not of prospective fate. 

Implicit in our definition of an organ-forming area, e.g., the heart- 
forming area, has been the assumption that the cells of that area of the 
embryo are destined to become specialized for the synthesis of the specific 
proteins characteristic of the definitive organ, and that the progressive 
localization of the heart-forming area is a localization of prospective myosin 
(and other cardiac protein) synthesizing cells. Yet the present findings 
demonstrate the progressive localization of the specific reactive groups of 
cardiac myosin, and not merely of prospective cardiac cells. Cardiac myo- 
sin is synthesized in the early embryo and can be detected in stage 3; in the 
embryo at stage 4, the protein is detected in all sections of the embryo 
analyzed (presumably in the mesoderm, although the present experiments 
do not prove this), being more concentrated in the mid and posterior, than 
in the anteriormost, regions. This finding is not in complete agreement with 
the results of explantation experiments. Spratt” studied the development of 
pulsating heart masses from anterior and posterior pieces of early blasto- 
derms cultivated im vitro, and concluded that although heart-forming po- 
tency is wide-spread in pre-streak embryos (see also ref, 2S), it later becomes 
localized in the median portion of the posterior half of the streak at stage 4. 
Rudnick™ ® likewise found that prior to its definitive localization at 
stage 5, the center of the heart-forming region lies posterior in the streak. 
The present findings do suggest that the synthesis of myosin is greater in 
the posterior than in the anterior part of the blastoderm; nevertheless some 
cardiac myosin is synthesized in the anterior part of the embryo. It 1s 
interesting to note that both Spratt and Rudnick used the classical plasma- 
embryo extract medium in their explantation studies, a medium which 
according to the reports cited above** ** might suppress cardiac differen- 
tiation, since in the absence of a statement to the contrary, the embryo 
extract is assumed to be an extract of the whole embryo, including the heart. 

In the embryo at stage 5, and even more clearly at stages 6-7, the lo- 
calization of cardiac myosin is in good agreement with the data of Rawles. 
Cardiac myosin is progressively restricted in distribution; at stage 5 the 


protein is limited at least to the region bounded by the anterior and pos 
. ” 
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terior borders of the heart-forming areas defined by Rawles, viz., 0.6 mm. 
forward of, and 0.45 mim. behind, the node. The localization of cardiac 
myosin in the medio-lateral axis is not completed until stages 6-7, although 
the protein is largely confined to the lateral heart-formuing areas at stage 5. 

Itis probable that the bulk synthesis of cardiac myosin is initiated only 
following the invagination of the prospective cardiac cells through the 
patent streak. ‘The first condition being fulfilled, myosin 1s synthesized in 
all areas of the embryo tested. 

Emphasis has already been placed on the fact that it is those cells al 
ready activately synthesizing cardiac myosin, as defined by its combining 
groups which are localized, and not merely those cells destined to synthesize 
the protein. Obviously the present findings do not fall readily into any cate 
gory with respect to current theories of differentiation. In fact, no single 
interpretation is offered here. If it 1s postulated either that myosin syn 
thesis 1s a self-limiting reaction or that other factors are limiting (e.g., 
availability of substrate, to cite one of the arguments frequently advanced), 
it then becomes necessary to explain the disappearance of myosin from the 
noneardiae areas of the blastoderm. At first thought it appears unlikely 
that the specificity of the protein might be altered and the protein redi 
rected into other channels. — Yet Daly and Mirsky“! have offered evidence of 
the synthesis of a protein “pool” in the exocrine pancreas, a pool from which 
highly specific proteins are rapidly synthesized as the secretory activity 
of the gland demands that the supply of specific enzymes be replenished. 
The degree of change of specificity in these transformations is not known. 
It is possible that the localization of the proteim is merely a sorting-out or 
segregation phenomenon, the result of cellular movements within the meso 
derm. Our present knowledge of morphogenetic movements within the 
mesoderm does not permit a eritical evaluation of this possibility, 

The question is posed as to whether it will be possible to characterize 
other organ-forming areas by immunochemical or biochemical methods. 
Obviously an organ or tissue must fulfill two criteria for successful experi 
mentation in this direction, both of which are admurably fulfilled in the 
heart-forming areas: it must be an organ which attains early morpho 
logical or functional status and for which a sound chemical basis of funetion 
has been established. Several of the organ-forming areas which fulfill 
the second criterion, e.g., the thyroid-forming area, fail when the first 
criterion 1s considered. — In our laboratory at present one possibility is being 
investigated (in addition to the analysis of the synthesis of cardiac actin as 
related to the present study), viz., an immunochemical study of the “pig 
ment-forming area’’ of the blastoderm. Further discussion and inter 
pretation of the findings and their bearing on the significance of the con 


cept of “organ-forming areas” 1s reserved pending the outcome of these 


experiments. 
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THE STRUCTURE OF THE MACRONUCLEUS OF PARAMECIUM 
AlURELIA* 

By R. F. KIimBayi 


BrotoGy Diviston, OAK RIpGE NATIONAL LABORATORY, OAK RipGe, TENNESSE! 


Communicated by T. M. Sonneborn, February 2, 19538 


Two types of bodies can be seen in the macronucleus of living Parame 


cium aurelia by phase-microscope observations.' One consists of spherical 


granules, about | uw in diameter, scattered more or less at random. The 
other type is seen only when the macronucleus is broken by compression. 
It consists of what appear to be many fine filaments, almost completely 
filling the nucleus. However, the filaments are so near the limits of resolu 
tion that their form cannot be determined with certainty. 

Paramecia exposed to lethal doses of ultraviolet radiation show a typical 
series of changes in the larger granules over the course of the several hours 
before death.” The granules clump together and fuse into one or a few 
large vacuolated masses. All stages in this fusion process can be followed 
in living material. The filaments are not visibly altered and remain uni 
formly distributed even after the granules have clumped into one large 
body. On the basis of the similarity of this behavior of the granules to the 
known reaction of nucleoli of other organisms to radiation, Kimball! 
tentatively concluded that they were nucleoli and the filaments were chro 
mosomes, 

Confirmation has been obtained for this interpretation by staining reac- 
tions. Animals of stock 90 of variety | were exposed to ultra-violet from 
a germicidal lamp and fixed in one per cent osmic acid after observations 
on living animals showed that the granules had fused into large masses 
Controls were also fixed at this time. Osmiuc acid was used as the fixative 
since observations under the phase microscope showed that, of all the 
fixatives tried, only osmic acid did not seriously alter the vacuolated fusion 
masses. The fixed material was dehydrated and embedded in paraffin. 
Two-micron-thick sections were cut and stained either with the Feulgen 
light green procedure following Rafaleo* and Conn and Darrow! or the pyro 
nin-methyl green procedure according to Lee.’ 

Examination of the stained sections of the irradiated animals showed that 
the fusion masses were stained with light green or pyronin, while the rest 
of the macronucleus was uniformly stained with Feulgen or methyl green. 
In control animals, light-green- or pyronin-stained material appeared to be 
scattered throughout the macronucleus, giving the Feulgen or methyl 
green staining a spongy appearance quite different from the uniform tex 
ture in the irradiated macronucleus. However, it was not possible to 
recognize with certainty distinct granules stained with the acid dye. Thus 
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the granules, which were so clearly visible in the living macronucleus, were 
no longer clear in the fixed nucleus, and the only evidence for them was 
scattered distribution of the acid staining. However, this scattered dis 
tribution and the results from the observations on living animals and on 
irradiated animals combine to make it clear that there are many small 
nucleoli seattered throughout the macronucleus. Apparently, the rest 
of its volume is largely occupied by many small chromosomes. 

It is of interest that an essentially identical interpretation of the struc 
ture of the macronucleus was reached by Breitschneider® on the basis of 
electron-microscope observations of sections of Paramecium caudatum. 
He saw many irregularly spherical bodies, connected by thin strands, and 
fewer larger bodies approximately | « in diameter. He interpreted the 
former as chromosomes, and the latter as nucleol. No proof was offered 
for these identifications. 

The macronucleus has been considered polyploid (see Sonneborn’ for 
review) or a multiple structure consisting of a number of subunits each of 
which is, perhaps, diploid.* The latter interpretation was based upon evti- 
dence from macronuclear regeneration taken in conjunction with genetic 
considerations. It 1s supported by Moses’ evidence® that there is only 
about 40 as much DNA in the macronucleus as in the micronucleus of 
Paramecium caudatum and by Nanney's work on mating types as given by 
Sonneborn.'’ In the present work, no evidence was found for subunits, 
and all the observations were consistent with the hypothesis of polyploidy 
with no subunits. It is not clear whether this means that the connections 
between the chromosomes of a subunit are too fine and diffuse to be ob- 
served with the techniques which have been used or whether the evidence 
which suggests subunits must be reinterpreted. 

The published observations on the ciliate macronucleus suggest at least 
two kinds of structure. Ina number of cases, scattered acid-stained granules 
are reported more or less as in the present paper. However, another 


group of observations indicate that large Feulgen-positive bodies are scat- 


tered throughout the nucleus, e.g., Turner'! for /Kuplotes. These large 
Feulgen-positive bodies are not readily homologized with the structures 
seen in the vegetative macronucleus of Paramecium. Commonly used 
fixatives, such as Schaudinn’s, produce considerable alteration in the struc 
ture of the macronucleus. Thus final judgment about Feulgen-positive 
bodies must be reserved until careful comparisons have been made between 
living and fixed material. Nevertheless, it seems probable that such bodies 
exist and must be taken into account in any generalized scheme of the 
structure and function of the macronucleus. The writer is unable to offer 
a satisfactory hypothesis to account for them. 


* Work performed under Contract No. W-7405-eng-26 for the Atomic Energy Comunis 


sion 





Vor. 39, 19538 ZOOLOGY: R. Fo KIMBALI 


' Kimball, R. F., Anat. Rec., 105, 543 (1949 
> Kimball, R. F., and Gaither, N., J. Cell. Comp. Physiol. 37, 211-234 (1951 
Rafalco, J. S., Stain Technol., 21, 91-93 (1946 
‘Conn, H. J., and Darrow, M. A., Staining Procedures Used by the Biological Stain 
Commission, Inc., Biotech Publications, Geneva, N. Y., ILB-141-1B-15, (1946) 
5 Lee, A. B., The Microtomist's Vade-Mecum, 10th ed., The Blakiston Co., Philadel 
phia, 165, (1937) 
® Breitschneider, L. H., Mikroskopie, 5, 257-269 (1950) 
7 Sonneborn, T. M., Ann. Rev. Microbiology, 3, 55-80 (1949 
* Sonneborn, T. M., Adv. in Genetics, 1, 263-358 (1947 
* Moses, M. J., J. Morphol., 87, 493-536 (1950 
‘8 Sonneborn, T. M., Cold Spring Harbor Symposia, 16, 481-503 (1951) 
' Turner, J. P., Univ. Catif. Publ. Zool., 33, 193-258 (1930), 











INFORMATION TO CONTRIBUTORS 


THE PROCEEDINGS is the official organ of the NATIONAL ACADEMY OF 
Scrences and of the NationaL Resgarcu Councit for the publication of 
brief accounts of important current researches of members of the ACADEMY 
and of the Councit and of other American investigators. The ProckEp- 
INGS will aim especially to secure prompt publication of original announce- 
ments of discoveries and wide circulation of the results of American re- 
search among investigators in other countries and in all branches of science. 


ARTICLES should be brief. The viewpoint should be comprehensive in 
giving the relation of the paper to previous publications of the author or 
of others and in exhibiting, where practicable, the significance of the work 
for other branches of science, Elaborate technical details of the work and 
long tables of data should be avoided, but authors should be precise in 
making clear the new results and should give some record of the methods 
and data upon which they are based.’ 


Manuscripts should be prepared with a current number of the Pro- 
CEEDINGS as a model in matters of form, and should be typewritten in 
duplicate with double spacing, the author retaining one copy. Illustrations 
should be confined to text figures of simple character, though more elabo- 
rate illustration may be allowed in special instances to authors willing to 
pay for their preparation and insertion. Particular attention should be 
given to arranging tabular matter in a simple and concise manner. 


REFERENCES to literature, numbered consecutively, will be placed at the 
end of the article and short footnotes should be avoided. It is suggested 
that references to periodicals be furnished in some detail and in general in 
accordance with the standard adopted for the Subject Catalogue of the 
International Catalog of Scientific Literature, viz., name of author, with 
initials following (ordinarily omitting title of paper), abbreviated name 
of Journal, with place of publication, series (if any), volume, inclusive 
pages, year. For example: Montgomery, T. H., J. Morph., Boston, 22, 
781-815 (1911); or, Wheeler, W. M., Konmigsberg, Schr. physik. Ges., 25, 
1-142 (1914). : 

Manuscripts will be accepted only from members of the Academy or 
from chairmen of the divisions of the National Research Council who will 
assume responsibility for the propriety and scientific standards of the paper 
and should be addressed to Proceedings of the NATIONAL ACADEMY OF 
Scrences, 2101 Constitution Ave., Washington, 25,D.C. Articles of which a 
member of the Academy 1s author or a co-author may extend to twelve prinied 
pages, those by non-members only to six printed pages; space above these 

imits will be billed to the author at $17.50 per page. 


Proor will ordinarily be sent,—with the understanding that charges for 
an author’s correction shall be billed to him. Authors are therefore re- 
quested to make final revisions on the typewritten manuscripts. 


Reprints should be ordered at the time of submission of manuscript 
They will be furnished to authors at cost, approximately as follows: 








No, Copies 50 100 


4 pp. ‘ 5. 
8 pp. t 11.20 
Covers 





4.80 








CONTENTS 


BACTERIOLOGY.—OXIDATION OF INTERMEDIATES OF THE TRICARBOXYLIC AcID 
Cyc.e py Exrracrs or Asotobacter Agile. . 
Z By Roy Repaske ond P, WwW. Wilson 
Brocumursrey. he Brocunmrcazy, Speciric MetHop ror Enzyme ISoLATION 
. By Leonard S. Lerman 
Cumersray, Tae Soueentry OF THE ‘Deoeceaitnitane ick or ALUMINUM Bro- 
MIDE IN MBTHYL BROMIDE IN THE PRESENCE OF SATURATED HYDROCARBONS 
, By Luther M. Foster and Charles A. Kra raus 
Cumusrny. —Tar Porenrtat UseruLness OF NATURAL TRITIUM... 
‘ . ByW. F. L pr 
Cutnnrarey. —Two Piedras Samer Commctmations OF POLYPEPTIDE CHAINS 
Invotvine Bors Cis anp TRANS Amipe Groups... 
By Linus Pauling and Robert B. Covey 
Cumaivrny: Two Rivpi.po-Ginent CONFIGURATIONS OF POLYPEPTIDE CHAINS, 
AND A NOTE ABOUT THE PLEATED SHEETS.  .. 
ae By Linus Pauling ond Robert B. Corey 
ENaneeunano. —On THE Encuniarson or Prastic Stip BANDS IN Fiat Bars IN 
TENSION TEsts . . -. . . By T. Y. Thomas 
ENGINEERING.—THE Evvect¢ or Conraxsemtuiry ON THE INCLINATION OF PLASTIC 
Sure Banpstn Frat Bars... ... . . . ByT.Y. Thomas 
ENGINBERING.—NONLINEAR MULTIPOLES .. . . . . By L. A. Zadeh 
Genetics.—A Srupy oF MUTABILITY IN Sraxpromycrn-Dapenvewt STRAINS OF 
Escherichia Cok. . . . . . «+ By Edgar L. Labrum 
Genetics.—INDUCED POLLEN Laraaus. From SEEDs or Datura Stramonium 
TREATED WITH THERMAL NEUTRONS ._. 
By J. L. Spencer, W. R. Sissies, aed A. P. Blakeslee 


Page 


225 © 


232 


236 


tf 


245 


Ganetics.~Tas Errect of THERMAL NEUTRON RADIATION ON THE CHROMO- 


somes OF Daiura . . 
a By Henry r. Yost, Je, W. R. Singleton, and.A. Fr. Blakeslee 
Marusizarics.—On Some VARIATIONAL PROBLEMS OCCURRING IN THE THEORY 
or D¥NamMIC PROGRAMMING ee 
See » By Richard Bellman, feving Glichsberg, end Oliver Gresi 
MATHEMATICS. — -FourigR TRANSFORMS OF TIME Series . .° By S. Bochner 
MatTHEMATICS,—ON THE BASIS THEOREM FOR FINITE ABELIAN Groups: IV ; 
By Jesse Douglas 
Marumatics. —Cauns WITH | Bounassrom Pomwrs oF 
‘ By Richard Otter and John Dunne 
MarHEmatics. —On THE Nosamn or Husmt Trees, I. . 
ae By Frank Harary end Sesrni E. Uhlenbeck 
Marnemarics.—Srupies 1 IN THE CONFORMAL MAPPING OF RIEMANN SURFACES, I. 
Eee ts By Maurice Heins 
Maruemartcs. —Tas Pasanaransarion AND ELeMent OF VOLUME OF THE UNI- 
TARY SympLectic GROUP... . . , By F. D. Murnaghan 
MATHEMATICS.—HEAT CONDUCTION ON “Amprre any RIEMANNIAN MANIFOLDS 
Seite, Cae eS aie Rye whit SN gas By D.\C. Spencer 
Matuematics.—ABsOLUTE VALUED ALGEBRAS . . . . By Fred B, Wright 
ZooLoGy.—AN ANALYSIS OF THE SYNTHESIS AND DisTRIBUTION OF THE Con- 
TRACTILE PROTEIN, MYOSIN, IN THE DEVELOPMENT OF THE HEART... 
hae By James D. Ebert 
Zoovocy. —Tas Srayerver o OF THE Matzonvcievs or Paramecium Aurelia. . 
a RS Se CMe eee catiieg, Botercs Gn. BRR Rie 





